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PREFACE 

This  is  the  second  in  a  series  of  Technical  Notes  (910-  )  entitled  "Self- 
Study  Manual  on  Optical  Radiation  Measurements."   It  contains  the  fourth  and 
fifth  chapters  of  that  Manual.   Along  with  Chapter  6,  which  has  already  been 
published  in  Technical  Note  910-3  (June  1977) ,  this  completes  the  first  six 
chapters  of  Part  I — Concepts. 

There  may  be  confusion  about  a  gap  in  the  numbering  of  the  appendices. 
There  were  two  appendices  in  TN  910-1  and,  when  TN  910-3  went  to  press  in 
June  1977,  we  were  planning  three  more  for  this  issue,  so  we  designated  those  in 
TN  910-3  as  numbers  6  and  7.   However,  the  third  appendix  for  this  issue  was  to 
have  been  on  "Apertures,  Stops,  Pupils,  Windows,  and  Baffles  in  Radiometry,"  and 
it  developed  that  this  topic  had  more  ramifications,  not  directly  pertaining  to 
chapters  4  and  5,  that  should  be  developed  in  greater  detail.   Accordingly,  in 
order  not  to  delay  further  the  publication  of  these  chapters  and  to  provide 
still  broader  and  more  adequate  coverage  of  the  theory  of  stops  in  relation  to 
radiometry,  it  was  decided  to  leave  out  the  appendix  and  revise  and  expand  the 
material  for  a  separate  chapter  in  the  near  future.   So  there  is  no  Appendix  5. 

For  background  information  on  the  whole  project  and  on  the  plans  for  the 
"Self-Study  Manual"  (SSM) ,  we  reproduce  here  (immediately  following)  the  Preface 
to  the  first  issue,  Technical  Note  910-1,  of  March  1976.   Although  the  exact 
details  of  the  outline  and  plans  for  the  individual  chapters  are  still  develop- 
ing, our  principal  aims  and  overall  plans  are  still  as  set  forth  in  that  Preface. 

Here,  in  Chapter  4,  we  are  primarily  concerned  with  developing  the  rest  of 
the  important  concepts  relating  to  the  distribution  of  optical  radiation  with 
respect  to  the  spatial  parameters  of  position  and  direction  and  the  spectral 
parameter  (wavelength,  frequency   v,   or  wave  number).   All  of  the  important 
radiometric  quantities  relating  to  these  parameters  are  defined  and  discussed, 
and  are  related  to  the  radiance  or  spectral  radiance  and  the  flux  or  spectral 
flux.   Important  interrelationships  are  also  explored  to  clarify  the  signifi- 
cance of  the  different  quantities  and  the  distinctions  between  them. 

In  Chapter  5  we  treat  the  measurement  equation,  the  heart  of  our  approach 
to  the  problems  of  radiometric  measurements.   The  basic  ideas  are  first  in- 
troduced by  deriving  the  measurement  equation  for  each  of  three  simple  illus- 
trative examples.   This  is  followed  by  a  more  general  discussion  based  on  the 
points  developed  in  the  sample  problems.   A  systematic  step-by-step  approach  is 
presented  for  setting  up  and  solving  any  measurement  equation  with  respect  to  the 
spatial  and  spectral  parameters.   Considerations  involving  the  other  radiation 
parameters  (polarization  and  time  or  frequency   f)   and  the  instrumental  and 
environmental  parameters  will  be  developed  in  future  chapters .   The  measurement 
equation  provides  the  basis  for  a  systematic  approach  to  the  difficult  multi- 
dimensional problems  of  radiometric  measurements.   Such  an  approach  is  an  aid 
for  planning  and  carrying  out  the  measurements  and  the  analysis  of  the  resulting 
data  and  for  preventing  possible  neglect  of  significant  parameters  that  affect 
those  results. 

We  are  grateful  to  many  for  valuable  comments  and  criticisms,  received  both 
informally  as  well  as  in  formal  reviews,  concerning  all  phases  of  the  SSM.   In 
connection  with  this  Technical  Note,  we  are  particularly  indebted  to  W.  L.  Wolfe 
and  F.  0.  Bartell  for  informal  reviews  of  portions  of  the  text  of  Chapter  4,  to 
L.  V.  Spencer  for  his  helpful  criticisms  and  discussion  of  Chapter  4,  and  to 
A.  T.  Hattenburg  and  J.  B.  Shumaker  for  valuable  discussions  and  comments  on  all 
phases  of  the  project.   Mrs.  Betty  Castle  has  done  an  outstanding  job  of  typing 
another  difficult  text,  including  the  effective  layout  of  text  and  figures.   And 
again,  we  acknowledge  H.  J.  Zoranski's  able  assistance  with  the  preparation  of 
the  figures. 

We  continue  to  earnestly  solicit  your  "feedback" — positive  or  negative.   We 
realize  that  most  of  you  are  probably  waiting  for  the  applications  chapters  and 


we'll  try  to  complete  some  of  them  as  soon  as  possible.   In  the  meantime,  how- 
ever, are  you  finding  these  first  chapters  on  basic  concepts  useful  and  helpful 
in  your  work?  Are  they  clear  and  complete  enough?   Have  you  found  any  inconsist- 
encies or  ambiguities  in  them,  or  between  them  and  other  material?  Have  they 
been  useful  for  indoctrination  of  new  employees  or  for  solving  particular  prob- 
lems?  Have  you  suggestions  for  improving  them?  Please  let  us  hear  from  you. 

Fred  E.  Nicodemus,  Editor 

Henry  J.  Kostkowski,  Chief, 
Optical  Radiation  Section 

PREFACE  to  NBS  TN  910-1 

This  is  the  initial  publication  of  a  new  series  of  Technical  Notes  (910)  entitled 
"Self-Study  Manual  on  Optical  Radiation  Measurements."   It  contains  the  first  three  chapters 
of  this  Manual.   Additional  chapters  will  be  published,  similarly,  as  they  are  completed. 
The  Manual  is  being  written  by  the  Optical  Radiation  Section  of  NBS.   In  addition  to  writing 
some  of  the  chapters,  themselves,  Fred  E.  Nicodemus  is  the  Editor  of  the  Manual  and  Henry  J. 
Kostkowski,  Chief  of  the  Section,  heads  the  overall  project. 

In  recent  years,  the  economic  and  social  impact  of  radiometric  measurements  (including 
photometric  measurements)  has  increased  significantly.   Such  measurements  are  required  in  the 
manufacture  of  cameras,  color  TV's,  copying  machines,  and  solid-state  lamps  (LED's).   Ultra- 
violet radiation  is  being  used  extensively  for  the  polymerization  of  industrial  coatings,  and 
regulatory  agencies  are  concerned  with  its  effects  on  the  eyes  and  skin  of  workers.   On  the 
other  hand,  phototherapy  is  usually  the  preferred  method  for  the  treatment  of  jaundice  in  the 
newborn.   Considerable  attention  is  being  given  to  the  widespread  utilization  of  solar  energy. 
These  are  just  a  few  examples  of  present  day  applications  of  optical  radiation.  Most  of 
these  applications  would  benefit  from  simple  measurements  of  one  to  a  few  per  cent  uncertain- 
ty and,  in  some  cases,  such  accuracies  are  almost  essential.   But  this  is  rarely  possible. 
Measurements  by  different  instruments  or  techniques  commonly  disagree  by  10%   to   50%,   and 
resolving  these  discrepancies  is  time-consuming  and  costly. 

There  are  two  major  reasons  for  the  large  discrepancies  that  occur.   One  is  that  optical 
radiation  is  one  of  the  most  difficult  physical  quantities  to  measure  accurately.   Radiant 
power  varies  with  the  radiation  parameters  of  position,  direction,  wavelength,  time,  and 
polarization.   The  responsivity  of  most  radiometers  also  varies  with  these  same  radiation 
parameters  and  with  a  number  of  environmental  and  instrumental  parameters,  as  well.   Thus, 
the  accurate  measurement  of  optical  radiation  is  a  difficult  multi-dimensional  problem.   The 
second  reason  is  that,  in  addition  to  this  inherent  difficulty,  there  are  few  measurement 
experts  available.   Most  of  the  people  wanting  to  make  optical  radiation  measurements  have 
not  been  trained  to  do  so.   Few  schools  have  had  programs  in  this  area  and  tutorial  and 
reference  material  that  can  be  used  for  self-study  is  only  partially  available,  is  scattered 
throughout  the  literature,  and  is  generally  inadequate.   Our  purpose  in  preparing  this  Self- 
Study  Manual  is  to  make  that  information  readily  accessible  in  one  place  and  in  systematic, 
understandable  form. 

The  idea  of  producing  such  a  manual  at  NBS  was  developed  by  one  of  us  (HJK)  in  the 
latter  part  of  1973.   Detailed  planning  got  under  way  in  the  summer  of  1974  when  a  full-time 
editor  (FEN)  was  appointed.   The  two  of  us  worked  together  for  about  one  year  developing  an 
approach  and  format  while  writing  and  rewriting  several  drafts  of  the  first  few  chapters. 
These  are  particularly  important  because  they  will  serve  as  a  model  for  the  rest  of  the 
Manual.   During  this  period,  a  draft  text  for  the  first  four  chapters  was  distributed,  along 
with  a  questionnaire,  for  comment  and  criticism  to  some  200  individuals  representing  virtu- 
ally every  technical  area  interested  in  the  Manual.   About  50  replies  were  received,  varying 
widely  in  the  reactions  and  suggestions  expressed.   Detailed  discussions  were  also  held  with 
key  individuals,  including  most  of  the  Section  staff,  particularly  those  that  will  be  writing 
some  of  the  later  chapters.   In  spite  of  the  very  wide  range  of  opinions  encountered,  all  of 
this  feedback  has  provided  valuable  guidance  for  the  final  decisions  about  objectives,  con- 
tent, style,  level  of  presentation,  etc. 
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In  particular,  we  have  been  able  to  arrive  at  a  clear  solution  to  difficult  questions 
about  the  level  of  presentation.   Both  of  us  started  out  with  the  firm  conviction  that,  with 
enough  time  and  effort,  we  should  be  able  to  present  the  subject  so  that  readers  with  the 
equivalent  of  just  elementary  college  mathematics  and  science  could  easily  follow  it.   That 
conviction  was  based  on  our  experience  of  success  in  explaining  the  subtleties  of  radiometric 
measurements  to  technicians  at  that  level.   What  we  failed  to  consider,  however,  was  that,  in 
making  such  explanations  to  individuals  we  always  were  able  to  relate  what  we  said  to  the 
particular  background  and  immediate  problem  of  the  individual.   That's  just  not  possible  in  a 
text  intended  for  broad  use  by  workers  in  astronomy,  mechanical  heat-transfer  engineering, 
illumination  engineering,  photometry,  meteorology,  photo-biology  and  photo-chemistry,  optical 
pyrometry,  remote  sensing,  military  infrared  applications,  etc.   To  deal  directly  and  explic- 
itly with  each  individual's  problems  in  a  cook-book  approach  would  require  an  impossibly 
large  and  unwieldy  text.   So  we  must  fall  back  on  general  principles  which  immediately  and 
unavoidably  require  more  knowledge  and  familiarity  with  science  and  mathematics,  at  the  level 
of  a  bachelor's  degree  in  some  branch  of  science  or  engineering,  or  the  equivalent  in  other 
training  and  experience. 

In  its  present  form,  the  Manual  is  a  definitive  tutorial  treatment  of  the  subject  that 
is  complete  enough  for  self  instruction.   This  is  what  is  meant  by  the  phrase  "self-study"  in 
the  title.   The  Manual  does  not  contain  explicitly  programmed  learning  steps  as  that  phrase 
sometimes  denotes.   In  addition,  through  detailed,  yet  concise,  chapter  summaries,  the  Manual 
is  designed  to  serve  also  as  a  convenient  and  authoritative  reference  source.   Those  already 
familiar  with  a  topic  should  turn  immediately  to  the  summary  at  the  end  of  the  appropriate 
chapter.   They  can  determine  from  that  summary  what,  if  any,  of  the  body  of  the  chapter  they 
want  to  read  for  more  details. 

The  basic  approach  and  focal  point  of  the  treatment  in  this  Manual  is  the  measurement 
equation.   We  believe  that  every  measurement  problem  should  be  addressed  with  an  equation 
relating  the  quantity  desired  to  the  data  obtained  through  a  detailed  characterization  of  the 
instruments  used  and  the  radiation  field  observed,  in  terms  of  all  of  the  relevant  param- 
eters.  The  latter  always  include  the  radiation  parameters,  as  well  as  environmental  and 
instrumental  parameters,  as  previously  pointed  out.   The  objective  of  the  Manual  is  to 
develop  the  basic  concepts  and  characteristics  required  so  that  the  reader  will  be  able  to 
use  this  measurement-equation  approach.   It  is  our  belief  that  this  is  the  only  way  that 
uncertainties  in  the  measurement  of  optical  radiation  can  generally  be  limited  to  one,  or  at 
most  a  few,  per  cent. 

Currently,  the  Manual  deals  only  with  the  classical  radiometry  of  incoherent  radiation. 
The  basic  quantitative  relations  for  the  propagation  of  energy  by  coherent  radiation  (e.g., 
laser  beams)  are  just  being  worked  out  [1,2,3,4].*   Without  that  basic  theory,  a  completely 
satisfactory  general  treatment  of  the  measurement  of  coherent  (including  partially  coherent) 
optical  radiation  is  not  possible.   Accordingly,  in  spite  of  the  urgent  need  for  improved 
measurements  of  laser  radiation,  we  won't  attempt  to  deal  with  it  now.   Possibly  this  sit- 
uation will  be  changed  before  the  current  effort  has  been  completed  and  a  supplement  on  laser 
measurements  can  be  added. 

As  stated  above,  we  first  hoped  to  prepare  this  Manual  on  a  more  elementary  level  but 
found  that  it  was  impossible  to  avoid  making  use  of  both  differential  and  integral  calculus 
of  more  than  one  variable.   However,  to  help  those  that  might  be  a  bit  "rusty"  with  such 
mathematics,  we  go  back  to  first  principles  each  time  a  mathematical  concept  or  procedure 
beyond  those  of  simple  algebra  or  trigonometry  is  introduced.   This  should  also  throw  addi- 
tional light  on  the  physical  and  geometrical  relationships  involved.   Where  it  seems  inap- 
propriate to  do  this  in  the  text,  we  cover  such  mathematical  considerations  in  appendices. 
It  is  also  assumed  that  the  reader  has  had  an  introductory  college  course  in  physics,  or  the 
equivalent. 


^Figures  in  brackets  indicate  literature  references  listed  at  the  end  of  this  Technical  Note. 


The  Manual  is  being  organized  into  three  Parts,  as  follows: 

Part  I.    Concepts 

Step  by  step  build  up  of  the  measurement  equation  in  terms  of  the  radiation  parameters, 
the  properties  and  characteristics  of  sources,  optical  paths,  and  receivers,  and  the  environ- 
mental and  instrumental  parameters.   Useful  quantities  are  defined  and  discussed  and  their 
relevance  to  various  applications  in  many  different  fields  (photometry,  heat-transfer  engi- 
neering, astronomy,  photo-biology,   etc.)  is  indicated.   However,  discussions  of  actual 
devices  and  measurement  situations  in  this  Part  are  mainly  for  purposes  of  illustrating 
concepts  and  basic  principles. 

Part  II.    Instrumentation 

Descriptions,  properties,  and  other  pertinent  data  concerning  typical  instruments, 
devices,  and  components  involved  in  common  measurement  situations.   Included  is  material 
dealing  with  sources,  detectors,  filters,  atmospheric  paths,  choppers  (and  other  types  of 
optical  modulators),  prisms,  gratings,  polarizers,  radiometers,  photometers,  spectroradiom- 
eters,  spectrophotometers,  etc. 

Part  III.   Applications 

Measurement  techniques  for  achieving  a  desired  level  of,  or  improving,  the  accuracy  of 
a  measurement.   Included  will  be  a  very  wide  variety  of  examples  of  environmental  and  instru- 
mental parameters  with  discussion  of  their  effects  and  how  to  deal  with  them.   This  is  where 
we  deal  with  real  measurements  in  the  real  world.   The  examples  will  also  be  drawn  from  the 
widest  possible  variety  of  areas  of  application  in  illumination  engineering,  radiative  heat 
transfer,  military  infrared  devices,  remote  sensing,  meteorology,  astronomy,  photo-chemistry 
and  photo-biology,  etc. 

Individual  chapter  headings  have  been  assigned  only  to  the  first  five  chapters: 

Chapter  1.    Introduction 

Chapter  2.    Distribution  of  Optical  Radiation  with  respect  to  Position  and 
Direction  —  Radiance 

Chapter  3.    Spectral  Distribution  of  Optical  Radiation 

Chapter  4.    Optical  Radiation  Measurements  —  a  Measurement  Equation 

Chapter  5.   More  on  the  Distribution  of  Optical  Radiation  with  respect  to 
Position  and  Direction 

Other  subjects  definitely  planned  for  Part  I  are  thermal  radiation,  photometry,  distribution 
with  respect  to  time,  polarization,  diffraction,  and  detector  concepts.   It  is  not  our 
intention,  however,  to  try  to  complete  all  of  Part  I  before  going  on  to  Parts  II  and  III.   In 
fact,  because  we  realize  that  a  great  many  readers  are  probably  most  interested  in  the 
material  on  applications  to  appear  in  Part  III,  we  will  try  to  complete  and  publish  some 
chapters  in  Parts  II  and  III  just  as  soon  as  adequate  preparation  has  been  made  in  the 
earlier  chapters  of  Part  I.   However,  because  our  approach  to  radiometry  differs  so  much  from 
the  traditional  treatment,  we  feel  that  unnecessary  confusion  and  misunderstanding  can  be 
avoided  if  at  least  the  first  nine  chapters  of  Part  I  are  published  first  and  so  are  avail- 
able to  readers  of  later  chapters. 

Finally,  we  invite  the  reader  to  submit  comments,  criticisms,  and  suggestions  for 
improving  future  chapters  in  this  Manual.   In  particular,  we  welcome  illustrative  examples 
and  problems  from  as  widely  different  areas  of  application  as  possible. 

As  previously  stated,  we  are  indebted  to  a  great  many  individuals  for  invaluable  "feed- 
back" that  has  helped  us  to  put  this  text  together  more  effectively.  Notable  are  the  inputs 
and  encouragement  from  the  Council  on  Optical  Radiation  Measurements  (CORM) ,  especially  the 
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CORM  Coordinators,  Richard  J.  Becherer,  John  Eby,  Franc  Grum,  Alton  R.  Karoli,  Edward  S. 
Steeb,  and  Robert  B.  Watson,  and  the  Editor  of  Electro-Optical  Systems  Design,  Robert  D. 
Compton.  In  addition,  for  editorial  assistance,  we  are  grateful  to  Donald  A.  McSparron, 
Joseph  C.  Richmond,  and  John  B.  Shumaker,  and  particularly  to  Albert  T.  Hattenburg. 

We  are  especially  grateful  to  Mrs.  Betty  Castle  for  the  skillful  and  conscientious 
effort  that  produced  the  excellent  typing  of  this  difficult  text.   We  also  want  to  thank 
Henry  J.  Zoranski  for  his  capable  help  with  the  figures. 

Fred  E.  Nicodemus ,  Editor 

Henry  J.  Kostkowski,  Chief, 
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SELF- STUDY  MANUAL  on  OPTICAL  RADIATION  MEASUREMENTS 
Part  I.   Concepts 

This  is  the  second  in  a  series  of  Technical  Notes  (910-)  entitled 
"Self-Study  Manual  on  Optical  Radiation  Measurements."   It  contains  the 
fourth  and  fifth  chapters  of  this  Manual.   Additional  chapters  will 
continue  to  be  published,  similarly,  as  they  are  completed.   The  Manual 
is  a  comprehensive  tutorial  treatment  of  the  measurement  of  incoherent 
radiation  that  is  complete  enough  for  self  instruction.   Detailed  chapter 
summaries  make  it  also  a  convenient  authoritative  reference  source. 

The  following  radiometric  quantities  are  defined  and  discussed  in 
Chapter  4:   radiant  energy,  radiant  exposure,  radiant  fluence,  radiant 
density,  radiant  intensity,  radiant  flux  (surface)  density,  irradiance, 
radiant  exitance,  radiant  fluence  rate,  radiant  sterisent,  as  well  as 
spectral  radiant  energy,  and  the  other  corresponding  spectral  quantities. 
In  particular,  each  quantity  is  related  to  the  quantities  previously 
introduced  in  Chapters  1-3:   radiant  flux  or  power,  radiance,  and  the 
corresponding  spectral  quantities.   Important  interrelationships  between 
the  different  quantities  are  also  presented,  particularly  where  they 
help  to  clarify  their  significance  and  the  distinctions  between  them. 
Treatment  of  one  of  the  most  important  interrelationships,  the  inverse- 
square  law,  also  covers  commonly  used  approximations  involving  actual 
(extended)  sources  and  receivers. 

The  measurement  equation,  central  to  our  approach  to  all  of  radiom- 
etry,  is  introduced  in  Chapter  5  through  three  illustrative  measurement 
problems:   a  spectral-radiance-comparison  measurement;  a  spectral-irra- 
diance  measurement  near  a  large  source;  and  an  irradiance  measurement 
with  a  wide-band  radiometer.  Normalization  of  spectrally  broad-band 
measurements  is  briefly  discussed.   A  general  discussion  of  the  meas- 
urement equation  summarizes  and  enlarges  on  the  points  brought  out  by 
the  examples.   A  set  of  orderly  steps  for  solving  the  measurement  equa- 
tion is  presented  and  the  limitations  of  the  measurement  equation  de- 
veloped in  this  chapter  are  summarized. 

Key  Words:      Measurement  equation;  optical  radiation  measurement; 
photometry;  radiometry;  spectroradiometry. 

Chapter  4.   More  on  the  Distribution  of  Optical  Radiation  with  Respect  to  Position 
and  Direction 

by  Fred  E.  Nicodemus 

In  this  CHAPTER.    We  introduce  all  of  the  remaining  radiometric  quantities  listed  in 
table  4-1  (a  slightly  revised  version  of  table  Al-3  in  Appendix  1  [5]).   In  Chapter  2  [5], 
we  concentrated  only  on  the  concept  of  radiance  as  the  basic  spatial  distribution  of 
radiant  flux  with  respect  to  both  position  and  direction  and  in  Chapter  3  [5]  we  similarly 
considered  only  spectral  radiance.   Now  we  treat  all  of  the  other  spatial  distributions  of 
flux,  and  of  energy  (time- integrated  flux).   Included  is  radiant  sterisent,  a  quantity  for 
characterizing  distributed  volume  sources.   In  every  case  the  quantity  is  defined  in  terms 
of  radiance  or  is  explicitly  related  to  the  radiance.   Interrelationships  between  the  dif- 
ferent quantities  are  also  given  to  clarify  their  meaning  and  significance  and  to  help  to 
distinguish  between  them.   A  particularly  important  interrelationship  is  the  inverse- 
square  law  which  gives  the  irradiance  element  at  a  given  distance  from  a  source  element 


(point)  of  known  radiant  intensity.   The  problems  involved  in  the  use  of  common  approxi- 
mations to  the  inverse-square  law,  applied  to  actual  (extended)  sources  and  receivers,  are 
discussed. 

In  Chapter  1  (on  p.  7  of  [5]),  we  said,  "It  cannot  be  too  strongly  emphasized  -  -  - 
that  everything  said  here  about  the  fundamentals  of  radiometry,  even  though  stated  in 
terms  of  watts,  applies  equally  to  all   forms  of  optical  radiation  measurements.   For 
example,  the  measurement  of  illumination,  for  application  to  vision  needs,  involves  all  of 
the  fundamentals,  not  just  those  discussed  in  the  chapter  on  photometry."  We  can  now  be 
more  explicit.   Compare  tables  4-1,  4-2,  and  4-3  (previously  tables  Al-3,  -4,  and  -6  in 
Appendix  1  [5]).   Note  that  in  both  table  4-2  and  4-3  there  is  an  entry  corresponding  to 
each  entry  in  table  4-1.   The  main  difference  is  apparent  in  the  units  or  unit-dimensions. 
Wherever  watts  [W]  appear  in  table  4-1,  they  are  replaced  by  lumens  [lm]  in  table  4-2  and 
by  quanta  per  second  [q«s   ]  in  table  4-3.   In  addition,  table  4-2  gives  some  equivalent 
units  that  are  also  commonly  used  for  the  photometric  quantities.   For  lack  of  space  and 
because   they  are  exactly  the  same,    the  defining  relations  in  the  third  column  of  table  4-1 
are  not  repeated  in  the  other  two  tables.   The  symbols,  too,  are  the  same  except  that, 
when  it  is  necessary  to  distinguish  between  them,  the  subscript   e  may  be  used  for  the 
radiometric  quantities,   v  for  the  photometric  quantities,  and  p   for  the  photon-flux 
quantities,  as  indicated  by  the  alternate  symbols  in  the  second  column  of  each  table. 
Ordinarily,  when  the  context  makes  it  clear  which  quantities  are  being  used,  the  subscripts 
may  be  omitted.   Also,  when,  as  in  this  chapter,  we  want  to  emphasize  the  common  geometri- 
cal relationships  that  apply  equally  to  all,  the  quantity  symbols  are  also  used  without 
subscripts.   We  will  present  the  material  in  this  chapter  in  terms  of  the  radiometric 
quantities,  with  unit-dimensions  in  terms  of  watts,  but  at!  of  the  geometrical  relation- 
ships are  equally  applicable  to  photometric  and  photon-flux  quantities .   It  is  only 
necessary  to  replace  each  term,  quantity  symbol,  and  unit-dimension  with  the  corresponding 
one  from  the  appropriate  table  to  show  this  explicitly  in  each  case. 

Many  readers  may  use  only  a  few  of  the  quantities  in  these  tables  and  those  of 
interest  to  them  may  change  from  time  to  time.   Accordingly,  they  may  prefer  to  turn 
directly  to  Chapter  5,  where  we  first  introduce  the  measurement  equation  and  deal  with  the 
concepts  relating  directly  to  the  measurement  of  some  radiometric  quantities.   Then,  as 
they  encounter  new  quantities  here  or  in  their  reading  elsewhere,  they  can  return  to  the 
pertinent  parts  of  this  chapter  where  each  one  is  defined  and  discussed  in  some  detail. 
To  find  the  chapter  where  each  quantity  is  defined,  see  the  numbers  in  parentheses  in  the 
first  column  of  table  4-1. 

"SIMPLE"  SPATIAL  DISTRIBUTIONS  (with  RESPECT  to  POSITION  or  DIRECTION) . 

Until  now,  we  have  analyzed  all  spatial  relations  connected  with  the  propagation  of  optical 
radiation  in  terms  of  the  basic  spatial-distribution  function,  namely  radiance  (as  a  func- 
tion of  ray  position  and  direction).   However,  a  glance  at  table  4-1  shows  that  there  are 
some  other  simpler  spatial  distributions  in  common  use  with  special  names.   These  are  dis- 
tributions with  respect  to  position  only  (flux  per  unit  area)  and  with  respect  to  direction 
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only  (flux  per  unit  solid  angle),  rather  than  with  respect  to  both  at  once.   Customarily, 
these  "simpler"  quantities  are  presented  and  discussed  first,  before  taking  up  the  more 
general  treatment  of  radiance.   However,  we  have  deliberately  reversed  the  usual  order 
because  it  has  been  our  experience  that  the  apparent  simplicity  can  be  misleading.   Lack 
of  attention  to  the  simultaneous  variation  of  ray-radiance  with  respect  to  both  position 
and  direction  in  actual  beams  of  radiation  can  often  lead  to  error.   For  consistently 
accurate  measurement  results,  the  sound  approach  is  to  look  at  every  situation  first  in 
terms  of  ray-radiances.  A  simpler  distribution  function,  with  respect  to  position  or 
direction  alone,  should  be  used  only  when  it  is  then  clear,  from  the  analysis  in  terms  of 
ray- radiances,  that  its  use  is  appropriate  and  will  not  introduce  any  unacceptable  approx- 
imations or  errors.   That  analysis  can  often  be  very  brief,  but  evevy   situation  should  be 
considered  first  in  this  way. 

In  subsequent  sections  we'll  take  up  each  of  these  "simple"  distributions,  in  turn, 
but,  here,  we'll  first  review  those  presentations  briefly  to  put  them  into  perspective. 
This  means  that  we'll  be  using  some  of  the  terms  before  they  are  fully  defined  and  dis- 
cussed in  detail.   For  that  reason,  it  may  be  worthwhile  to  reread  this  section  after 
going  through  the  rest  of  the  chapter,  to  more  fully  appreciate  the  significance  of  the 
intercomparisons  between  the  different  quantities. 

We  take  up  first  the  distribution  of  flux  with  respect  to  position  only,  in  the  form 
of  flux  per  unit  area  as  a  function  of  position.   Here  we  find  two  possibilities.   To 
distinguish  them,  we  refer  to  one  as  a  "directed-surface  distribution"  and  the  other  as  an 
"omni-directional-surface  distribution."  An  example  of  the  first  would  be  the  quantity 
measured  by  a  plane  detector  responding  to  incident  radiation  from  the  hemisphere  above 
it;  the  second  would  be  that  measured  by  a  spherical  detector  responding  equally  to  inci- 
dent radiation  from  all  directions  over  a  full  sphere  surrounding  it. 

The  first,  the  directed-surface  distribution,  is  called  by  the  CIE-IEC  [6]  the 
"radiant  flux  (surface)  density  (at  a  point  of  a  surface)."  They  also  use  the  terms 
"irradiance"  [for  incident  radiant  flux  (surface)  density]  and  "radiant  exitance"  [for 
exitent — leaving  the  surface — radiant  flux  (surface)  density].   However,  they  have  no 
nomenclature  at  all  for  the  second,  the  omni-directional-surface  distribution.   It  is  used 
primarily  by  photobiologists  and  photochemists,  nuclear  physicists  and  engineers,  and 
others,1  many  of  whom  refer  to  it  as  the  "fluence  rate."  Both  kinds  of  surface-area 


^Nuclear  scientists  generally  use  the  term  "energy  current  density"  for  the  directed- 
surface  distribution  that  we  call  "radiant  flux  (surface)  density"  and  they  use  "flux 
density"  or  "energy  flux  density"  for  the  omni-directional-surface  distribution  that  we 
call  "radiant  fluence  rate"  [7].   However,  "fluence  rate"  is  also  given  as  an  acceptable 
alternate  for  the  latter  [8,9],  so  our  choice  of  that  term  does  not  create  any  direct 
conflict.   With  respect  to  the  term  "flux  density",  however,  there  is  an  unfortunate,  but 
unavoidable,  ambiguity  in  the  existing  literature  that  even  the  unit-dimensions  won't 

(Footnote  continued  on  page  8) 


resolve  because  both  the  CIE-IEC  [6]  "flux  (surface)  density"  and  the  nuclear  [7,8,9] 

"flux  density"  (or  "fluence  rate")  are  measured  in  units  of  propagated  quantity  per  unit 

-1  -2         -2 
time  per  unit  area,  e.g.   [J*s   *m  ]   or  [W*m  ].   All  we  can  do  is  to  warn  that  when  a 

nuclear  scientist  designates  a  quantity  as  a  "flux  density"  he  probably  means  an  omni- 
directional-surface distribution  that  we  call  "fluence  rate"  rather  than  the  directed- 
surface  distribution  that  we  call  "flux  (surface)  density"  and  which  he  will  more  likely 
call  "current  density".   The  only  way  we  know  to  verify  this,  in  the  absence  of  explicit 
definitions,  is  through  relationships  that  may  be  given  between  each  quantity  and  other 
radiometric  quantities.   That's  one  of  the  important  reasons  for  presenting  such  inter- 
relationships between  the  different  radiometric  quantities  in  considerable  detail  in  this 
chapter. 


_2 
distribution  are  measured  in  units  of  watts  per  square  meter  [W*m  ].  Up  to  this  point, 

and  in  tables  4-1  and  -2,  we've  used  the  descriptive  phrases  "directed-surface"  and  "omni- 
directional-surface" to  emphasize  the  distinction  between  these  two  forms  of  flux  distribu- 
tion with  respect  to  position.   Hereafter,  and  in  table  4-3,  we  use  the  less  cumbersome 
terms  "flux  (surface)  density"  and  "fluence  rate",  without  the  descriptive  phrases,  except 
where  they  may  be  introduced  for  emphasis.   The  parenthetical  term  "surface"  is  retained, 
however,  because  the  term  "density"  alone  more  often  refers  to  a  volume  concentration. 

The  radiant  flux  (surface)  density  is  the  flux  per  unit  area  through  a  surface-area 
element  of  fixed  orientation,  from  (or  to)  all  (or  any)  directions  within  the  hemisphere 
on  one  side  of  it.   The  radiant  fluence  rate  is  the  flux  per  unit  cross-sectional  area 
incident  on  a  spherical  volume  element  (which  has  the  same  cross-sectional  surface  area 
in  all  directions  lying  within  a  complete  sphere  centered  on  the  volume  element) .   The 
flux  (surface)  density  is  the  appropriate  quantity  to  relate  to  a  reference  surface  that 
intersects  a  radiation  beam,  such  as  the  (imaginary)  plane  surface  across  each  one  of  a 
pair  of  beam-defining  apertures.   Situations  involving  the  fluence  rate  are  those  where 
radiation  is  incident  throughout  a  large  solid  angle  (up  to  a  complete  sphere)  without  a 
beam-defining  aperture  at  the  receiver.  An  example  would  be  a  small  amount  of  matter 
suspended  in  space  where  it  is  exposed  to  irradiation  from  all  sides,  such  as  a  microbe  in 
an  irradiated  transparent  culture  medium  or  a  small  volume  element  at  a  point  in  the 
lasing  gas  of  an  "optically  pumped"  laser. 

Note  that  flux  (surface)  density  may  refer  either  to  incident  radiation  (irradiance) 
or  exitent  radiation  (exitance) — either  receiver  or  source.   Fluence  rate,  on  the  other 
hand,  involves  only  incident  radiation — a  receiver  quantity.   [Exitent  (emitted  or  scat- 
cattered)  radiation  from  a  volume  source,  not  referred  to  a  reference  surface,  will  be 
treated  later  when  we  take  up  radiant  sterisent.] 

The  last  of  these  "simple"  distributions,  the  distribution  of  radiant  flux  with 
respect  to  direction  only,  in  the  form  of  exitent  flux  per  unit  solid  angle  as  a  function 
of  direction,  in  contrast  to  fluence  rate,  involves  only  exitent  radiation;  it  is  a  source 
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quantity.  While  it  is  defined  for  any  source,  including  extended  sources  as  we'll  see 
shortly,  it  isn't  very  useful  except  when  the  source  can  be  treated  as  a  "point"  source; 
i.e.,  when  the  source  is  far  enough  away  so  that  the  irradiance  it  produces  on  an  element 
of  receiving  surface  is  inversely  proportional  to  the  square  of  its  distance  from  that 
receiver  element.   This  directional  distribution  of  radiant  flux  is  called  the  "radiant 
intensity  (of  a  source  in  a  given  direction)"  by  the  CIE-IEC  [6],  and  they  use  that  term 
consistently  for  just  the  one  quantity.   Unfortunately,  however,  there  are  many  others  who 
use  the  term  "intensity"  in  a  variety  of  ways,  resulting  in  widespread  ambiguity  and 
confusion.   For  example,  many  physicists  use  the  term  "intensity"  for  the  flux  (surface) 
density.   Others,  particularly  astrophysicists  and  meteorologists,  use  the  term  as  a 
synonym  of  radiance,  although  many  of  them  add  the  modifier  "specific"  so  that  they  speak 
of  the  "specific  intensity"  of  flux  per  unit  area  and  solid  angle.   There's  just  no  way  to 
avoid  this  unfortunate  situation  in  the  existing  literature,  so  here  is  one  place  where 
it's  particularly  important  to  pay  close  attention  to  the  units  (unit-dimensions)  given — 
for  any  quantity  called  "intensity."  In  this  Manual,  we  consistently  follow  the  CIE-IEC 
[6]  and  speak  of  "intensity"  only  when  we  mean  flux  per  unit  solid  angle  in,  e.g.,  watts 
per  steradian   [W«sr   ]. 

We  also  discuss  in  some  detail  the  considerations  governing  the  valid  application  of 
the  quantity  called  radiant  intensity,  the  circumstances  under  which  a  real  source  may  be 
treated  usefully  as  a  "point  source";  i.e.,  one  where  the  irradiance  produced  by  it  on  a 
distant  element  of  a  receiver  obeys  the  inverse-square  law.   This  analysis  is  in  terms  of 
radiance  and  throughput1  and  it  also  introduces  the  engineering  "configuration  factors" 
that  are  closely  related  to  the  concepts  of  throughput  and  projected  solid  angle2.   The 
relationships  between  configuration  factors  and  throughputs  and  projected  solid  angles  are 
developed  in  Appendix  3..  With  these  relationships,  the  extensive  tables  of  configuration 
factors  that  are  available  can  be  used  to  conveniently  evaluate  the  throughputs  and/or 
projected  solid  angles  for  a  great  many  beam  configurations,  between  sources  and  receivers, 
covered  by  those  published  tables. 

DIRECTED- SURFACE  DISTRIBUTIONS:   RADIANT  FLUX  (SURFACE)  DENSITY;  IRRADIANCE;  RADIANT 
EXITANCE .   There  are  a  number  of  practical  situations  where  the  quantity  of  significance 
is  just  the  amount  of  flux  per  unit  area  flowing  through  (into  or  out  of)  a  surface, 
regardless  of  the  orientation  of  the  rays  along  which  the  energy  is  propagating.   For 
example,  in  determining  the  required  size  for  a  black  surface  to  dissipate  a  certain 
amount  of  power  (flux)  by  radiation,  an  important  quantity  is  the  amount  of  flux  per  unit 
area,  the  number  of  watts  per  square  meter   [W*m  ],   that  the  surface  emits.   For  another 
example,  the  flux  from  a  very  distant  radiation  source  is  almost  always  uniformly  distrib- 
uted in  rays  of  equal  radiance  over  any  small  plane  receiving-surface  area  that  intersects 

^ee  eq.  (2.27)  [5]. 
2See  Appendix  2  [5]. 


a  narrow  portion  of  the  beam  from  that  source.   Accordingly,  the  amount  of  received  power 
or  flux  is  directly  proportional  to  the  receiving  area;  it  is  the  product  of  that  plane 
area  and  the  incident  flux  (surface)  density  or  irradiance,  which  we  will  now  define. 

In  order  to  examine  the  distribution  of  radiation  with  respect  to  position,  consider 
a  photocell  lying  on  a  bench  or  desk  below  two  long  tubular  fluorescent  ceiling  lamps,  as 
illustrated  in  figure  4.1.   The  radiation  (light)  from  the  lamps  reaches  the  photocell 
through  a  circular  hole  in  a  flat  opaque  baffle,  placed  on  top  of  the  photocell  so  that 

its  center  point  P  is  just  above  the  center  of  the  photocell.   Baffles  with  holes  of 

2 
different  areas   AA  [m  ]   are  used,  always  centered  at  P   and  all  small  enough  so  that 

all  of  the  flux  through  them  reaches  the  sensitive  surface  of  the  photocell.   The  output 

from  the  photocell  is  then  a  measure  of  the  flux  A$  [W]   through  the  hole  from  both  lamps. 


LAMP 


LAMP 


PHOTOCELL 


Figure  4.1.    Experiment  to  develop  the  concept  of  irradiance. 
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In  fact,  we'll  assume  that  the  photocell  has  an  ideal  response  that  is  the  same  for  a 
given  amount  of  incident  flux  A$  reaching  its  sensitive  surface  at  any  point  and  from 
any  direction.   As   AA  is  reduced,  by  using  smaller  and  smaller  holes  in  the  baffle,  the 
flux  A$,   as  measured  by  the  photocell  output,  also  decreases,  while  the  quotient  A$/AA 
approaches  a  constant.   That  limiting  value,  toward  which  the  quotient  converges  as  the 
area  AA  becomes  vanishingly  small,  can  also  be  described  as  the  quotient  of  the  element 
of  flux  d$(x,y)   through  a  surface  element   dA(x,y)   containing  the  point  x,y  by  the 
area  dA  of  the  element,  and  is  called  the  flux  (surface)  density  (at  a  point  of  a 
surface)  * 

W(x  y)   =   L±m  M(x»y}    =   d$(x,y)   [Wm"2]  (4.1) 

The  orientation  of  the  rays  passing  through  dA  need  not  be  specified.  They  may  be 
concentrated  in  a  small  solid  angle,  as  part  of  a  well-collimated  beam  (e.g.,  that  from  a 
very  distant  source,  such  as  a  star)  or,  at  the  other  extreme,  they  may  be  spread  more  or 
less  evenly  over  the  complete  hemisphere  of  directions  above  dA  in  a  widely  diffuse 
field  of  radiation.  In  the  more  general  case,  if  dA  is  an  element  of  a  curved  surface, 
the  hemisphere  is  that  above  the  tangent  plane  that  contains   dA. 

In  order  to  cover  all  possibilities  when  ray  directions  are  confined  to  the  hemisphere 
on  one  side,  only,  of  a  surface  element,  we  also  take  account  of  the  sense  in  which  propa- 
gation occurs  along  a  ray.   It  is  incident  when  it  is  toward  the  surface;  exitent  when  it 
is  away  from  the  surface.   As  we'll  see  when  we  take  up  thermal  radiation,  we  are  some- 
times concerned  with  simultaneous  propagation  in  both  senses,  incident  and  exitent,  along 
each  ray  direction,  as  well  as  with  the  net  energy  transfer  given  by  the  difference  between 
them.   It  is  this  net  value  that  is  measured  by  most  detectors.   This  can  be  a  particularly 
important  consideration  in  heat-transfer  configurations  near  thermal  equilibrium  where 
each  of  two  surfaces  at  nearly  the  same  temperature  may  emit  substantially  toward  the 
other  but  where  there  is  little  or  no  net  exchange  of  energy  between  them. 

If  the  sense  of  propagation  along  all  rays  is  toward  a  surface  element  dA  (from  one 
side  of  dA) ,   the  incident  radiant  flux  (surface)  density  W  is  usually  designated  as 
the  irradiance,  with  the  symbol  E,(E  )   [6].   Similarly,  for  propagation  away  from  dA, 
the  exitent  radiant  flux  (surface)  density  W  leaving  the  surface  (on  one  side)  is  usually 
designated  as  the  radiant  exitance,  with  the  symbol  M, (M  )   [6]. 

In  previous  chapters  [5]  we  pointed  out  that  radiance  and  spectral  radiance,  being 
derivatives,  are  quotients  of  vanishingly  small  quantities  that,  due  to  limiting  noise 
levels  in  all  measuring  equipment,  can  never,  even  in  principle,  be  measured  exactly. 
More  seriously,  reduction  of  the  quantities  involved  reaches  a  point  where  the  basic 


^he  CIE-IEC  [6]  has  no  symbol  for  this  quantity  other  than  the  more  restricted  M  or  E. 
We  have  adopted  W,   formerly  widely  used  for  radiant  exitance  [now  M,(M  )]. 
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assumptions  of  the  geometrical-optics  model  of  ray  propagation  are  no  longer  valid.   The 
same  is  true  with  respect  to  flux  (surface)  density,  including  exitance  and  irradiance, 
and  with  respect  to  all  of  the  other  distributions  treated  in  this  chapter.   However,  our 
use  of  calculus,  based  on  the  assumption  of  an  underlying  continuous  distribution  of 
infinite  resolution,  provides  useful  results  that  agree  well  with  experimental  measure- 
ments as  long  as  we  are  careful  not  to  try  to  apply  those  results  outside  the  domain  of 
geometrical  optics. 

In  Chapter  2  [5],  we  defined  the  radiance  at  the  point  x,y  on  a  reference  surface 
in  the  direction  6,.<|>  of  an  intersecting  ray  through  that  point  (through  the  surface 
element  dA  at  that  point — see  figure  4.2)  as 

2 
T  e    Q  j.\      d  <Kx,y,6,(j>)   i-tt  ~2   -In  m   i  /  \ 

L(x,y,e,<f>)   =    dA.c;^.^    [W«m   -sr   ]  (2.14) 

From  this  defining  equation  for  radiance,  we  were  also  able  to  write  the  expression  for 
the  element  of  flux  associated  with  a  ray  of  radiance  L(x,y,0,<(>)   as 

d$(x,y,e,<j>)   =  L(x,y,9,<j))'COse«d(WA  (2.24) 

=  L(x,y,e,<)>)'dft«dA  [W] .  (4.2) 

Accordingly,  the  total  flux  through  the  surface  element  dA  at  the  point  x,y  is 

d$(x,y)   =  dA'/  L(x,y,e,<J>)«dfi  =  dA*/  L(x,y,0  ,<j>)  -cos0  'dw 

=  dA*/  /  L(x,y,6,(f))'cos0«sin0'de'd(j)  [W]  ,  (4.3) 

<(,  0 

where  the  limits  of  integration  for  0   (in  the  range  0  -  tt/2  [rad])   and  <f>   (in  the 
range  0  -  2tt  [rad])   include  all  rays  of  the  beam  of  interest  that  pass  through  dA  at 
x,y.   The  flux  (surface)  density  (irradiance  or  radiant  exitance,  depending  on  whether 
propagation  is  toward  or  away  from  dA  on  the  side  in  question)  through  the  surface 
element  dA  at  the  point  x,y   is  then 

W(x,y)  =      d$(x,y)/dA  (4.1) 

=  /  L(x,y,0  ,<f>)  «dfi  =  /  L(x,y,0  ,<}>) 'COsO'dco 

=  /  /  L(x,y,0,cf>)«cos0«sin0'd0'd(|>  [W-m   ]  (4.4) 

$  0 

where,  again,  the  limits  of  integration  include  all  of  those  rays  of  the  beam  that  pass 
through  dA  at  the  point  x,y.  In  particular,  if  the  beam  fills  the  entire  hemisphere 
above  the  element  dA   (above  the  tangent  plane  containing   dA) ,   this  becomes 
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Figure  4.2.   Geometry  of  ray-surface  intersection  (for  the  definition 
of  radiance) .   (This  figure  previously  appeared  as 
figure  2.7  [5].) 
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W(x,y)   =  /  L(x,y,6,(f>)'dfi  =  /  L(x,y,6  ,$)  -cose'dw 

h  2tt 

2tt  it/2  _„ 

=  /   /   L(x,y,6,(())«cose«sine«de«d<()  [W-m   ].  (4.4a) 

o  o 

Here,   J     denotes  integration  over  a  full  hemisphere  of   2tt  [sr]   with  respect  to  solid 

h 
angle.   Of  course,   E(x,y)   (irradiance)   or  M(x,y)   (exitance)   can  always  be  substituted 

for  W(x,y)   {flux  (surface)  density},   in  eqs.  (4.4)  and  (4.4a)  or  elsewhere,  whenever  it 

is  appropriate  to  be  more  explicit  about  whether  propagation  is  respectively  toward  or 

away  from  the  surface  element   dA. 

Although,  as  we  have  pointed  out,  flux  (surface)  density  (including  irradiance  and 
exitance)  is  not  a  directional  quantity,  it  is  useful  to  recognize  the  directionality  of 
the  element   of  flux  (surface)  density  by  relating  it  to  the  flux  element  (radiance  times 
throughput  element)  associated  with  a  single  ray  in  the  direction  0,<j>   through  the  surface- 
area  element  dA  at  the  point  x,y   (see  figure  4.2): 

dW(x,y,e,(f>)   5   d$(x,y,e,<j>)/dA 

=  L(x,y,e,<j>)-d9/dA  [W'irf2],  (4.5) 

where  the  throughput  element  {eq.  (2.28)  [5]}   is 

d0   =   dA«dft  =   dA«cos6'dw 

2 
=  dA* cos6 • sine 'dO'd^  [m  «sr]  (4.6) 

so  that 

dW(x,y,e,c(>)   =  L(x,y,e,(j)) -dfi  [W-m~  ].  (4.7) 

Again,  of  course,   dE(x,y,6  ,<j>)   or  dM(x,y,6  ,<)>)   may  replace  dW(x,y  ,6  ,<)>)   to  designate, 
explicitly,  the  sense  of  propagation  (toward  or  away  from  dA)   along  the  ray  in  the 
direction  0,<f>   through  dA  at  the  point  x,y.   This  also  provides  an  alternate  deriva- 
tion of  eq.  (4.4),  the  last  two  lines  of  which  follow  directly  from  eq.  (4.7),  since 

W(x,y)   =   /  dW(x,y,6,<j))  [W-m"2]  .  (4.8) 
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FLUX  (SURFACE)  DENSITY  in  ISO-RADIANCE  BEAMS.    The  relationship  between  radiance  and  flux 
(surface)  density  across  any  intersecting  reference  surface,  is  greatly  simplified  for  an 
iso-radiance1  beam  in  which  the  value  of  radiance  is  the  same  for  all  rays .   Then  L   is  a 
constant  that  may  be  brought  outside  the  integral  in  eq.  (4.4)  and  the  flux  (surface) 
density  through  a  surface  element   dA  at  any  point  x,y  becomes 

W  =  L-n  [W-m~2].  (4.9) 

It  is  the  product  of  the  constant  radiance  L  and  the  projected   solid  angle  fi  =  J  cosO'dco 

to 
{eq.  (2.31)  [5]},   where  w  is  the  solid  angle  filled  by  the  rays  of  the  beam  that  inter- 
sect at  the  point  x,y  where  the  quantity  is  evaluated.   In  particular,  when  all  rays 
into,  or  from,  the  full  hemisphere  above  an  iso-radiance  surface  element  or  plane  surface 
are  included  in  the  beam, 

W  =  tt'L  [W-m~2],  (4.10) 

since  the  projected  solid  angle  of  a  full  hemisphere  is  it  [sr]   {see  eq.  (A2-15)  in 
Appendix  2  [5]}. 

Iso-radiance  beams  and  their  properties  are  very  important  in  radiometry  (including 
photometry) .   It  is  obvious  from  the  foregoing  that  analysis  is  greatly  simplified  when 
iso-radiance  conditions  exist.   Even  when  all  rays  are  not  of  exactly  the  same  radiance 
but  are  only  approximately  so,  the  use  of  an  average  value  as  though  it  were  everywhere 
the  same  provides  approximate  results  that  are  useful  for  many  purposes.   However,  this 
has  led  to  such  widespread  treatment  of  almost  all  problems  in  terms  of  source-surface 
radiant  exitance  and  receiver  irradiance,  particularly  in  the  fields  of  heat-transfer  and 
illumination  engineering,  that  the  directional  distribution  of  radiance  is  too  easily 
overlooked  in  situations  where  it  is,  in  fact,  significant.   Accordingly,  we  reemphasize 
that  the  sound  approach  is  always  to  analyze  a  situation  first  in  terms  of  the  radiance  of 
individual  rays  and  its  spatial  distribution  and  to  employ  one  of  the  "simpler"  distribu- 
tions, in  terms  of  position  or  direction  alone,  only  after  that  analysis  has  shown  clearly 
that  this  can  be  done  without  unacceptable  approximation  or  error. 


1 Iso-radiance,  rather  than  lambertian,  because,  strictly,  a  lambertian  surface  is  one  that 
is  isotropic  (the  same  radiance  in  all  directions — at  least  within  the  beam  of  interest) 
at  each  point,  although  it  might  be  non-uniform  (varying  from  point  to  point).   If  it  is 
non-uniform,  a  second  surface  cutting  the  beam  at  a  distance  from  the  first  will  no  longer 
be  lambertian.   The  rays  crossing  at  any  one  point  of  the  second  surface  will  come  from 
different  points  of  the  original  surface,  so  the  radiance  of  a  point  on  the  second  surface 
is  no  longer  isotropic.   An  iso-radiance  surface  or  beam  is  both  uniform  and  labertian;  it 
has  the  same  radiance  everywhere  and  in  all  directions  so  that  every  surface  cutting  the 
beam  anywhere  is  lambertian  (within  the  beam) . 
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(0,0,z) 


(x,y,z)  or 
(p>0,<t>) 


(x,y,0) 


Figure  4.3.    A  small  sphere  of  radius   Ar   and  cross  section 

2 
Aa  =  Tr-(Ar)  ,   centered  at  the  point  x,y,z 

(or   p,6,4>). 
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OMNI-DIRECTIONAL- SURFACE  DISTRIBUTION:   FLUENCE  RATE.    The  directed-surface  distributions, 
defined  by  eq.  (4.1),  are  seen,  by  eq.  (4.4),  to  depend  on  the  orientation  of  the  reference- 
surface  element  dA.   Accordingly,  they  are  not  scalar  field  quantities,  since  they  do  not 
have  a  unique  value  at  each  point.   The  omni-directional-surface  distribution,  on  the  other 
hand,  involves  no  reference  surface;  it  does  have  a  unique  value  at  each  point  in  space, 
making  it  a  scalar  field  quantity.   The  fluence  rate,  as  we  call  it,  is  defined  with  ref- 
erence to  the  cross-sectional  area  da  of  a  spherical  volume  element  (which  is  the  same  in 
all  directions)  rather  than  the  area  dA  of  a  surface  element  of  fixed  orientation. 
Sometimes  it  is  called,  instead,  the  "scalar  irradiance." 

In  order  to  analyze  the  concept  of  omni-directional-surface  distribution  or  fluence 
rate  by  means  of  a  "thought  experiment,"  consider  a  small  sphere  of  radius  Ar  with  its 
center  at  the  point  x,y,z  in  space,  as  shown  in  figure  4.3.   The  area  of  the  circular 
plane  surface  through  x,y,z  perpendicular  to  any  incident  ray,  i.e.,  the  area  of  the 
circular  cross-section  of  the  sphere,  is  then 

a  =  ir-(Ar)2  [m2].  (4.11) 

We  also  postulate  that,  in  some  way,  we  can  determine  the  total  flux  A$  [W]   incident  on 
the  sphere  from  all  directions.   (This  might  be  approximated  by  replacing  the  sphere  with 
an  ideal  black  spherical  receiver  and  making  calorimetric  measurements  based  on  its  temper- 
ature history,  etc.)   Then,  if  we  keep  the  center  of  the  sphere  fixed  at  x,y,z  and 
reduce  the  cross-sectional  area  Aa  by  reducing  the  radius  Ar,   the  fluence  rate  at  the 
point  x,y,z   is  defined  as  the  limit,  as   Aa  approaches  zero,  of  the  quotient   A$/Aa: 


FCx.y.z)   E   L^  A*(x,y,z)   E   d*(x,y,z)  {       -2 
t   '  -"      .  Aa-*0     Aa  da  ' 


(4.12) 


where 


_2 
F  (x,y,z)  [W«m  ]   is  the  fluence  rate  at  the  point  Xjyjz;1 


and 


d$(x,y,z)  [W]   is  the  flux  or  power  incident  on  a  spherical  volume 

2 
element,  of  cross-section  da  [m  ] ,   centered  at  the  point  x,y,z. 

No  limitation  is  placed  on  the  direction  of  arrival  of  the  incident  radiant  flux  d$, 


1Neither  the  CIE-IEC  [6]  nor  the  ANSI  [10,11]  has  either  a  term  or  a  symbol  for  this 
quantity,  or  for  the  related  quantity,  fluence.   After  choosing  the  symbol  F   for  fluence 
[12]  which  we  discuss  later  in  this  chapter,  we  preferred  to  avoid  an  additional  symbol,  as 
used  elsewhere  [8],  for  fluence  rate  by  adopting  F   [=  dF/dt]   since,  as  we  will  see 
presently,  the  fluence  rate  is  the  time  derivative  of  the  fluence  at  the  same  point.   In 
this  Manual,  we  employ  the  notation  X  =  dX/dp  for  any  radiometric  quantity  X  and  the 
following  radiation  parameters  p:   position,  direction,  any  spectral  parameter  (wavelength, 
wave  number,  frequency  v,  etc.),  or  time  or  frequency  f  <<  v. 
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It  may  be  incident  from  any  direction  which,  unless  explicitly  limited  in  a  particular 
case,  includes  the  entire  sphere  of  directions  surrounding  the  point  x,y,z,   a  full  4ir-[sr] 
solid  angle. 

Next,  let's  consider  the  contribution  to  the  fluence  rate,  defined  in  eq.  (4.12),  by 
the  incident  element  of  flux  d$(x,y,z,0  ,<j>)   associated  with  a  single  ray,  of  radiance 
L(x,y,z,6  ,<f>)   through  the  point  x,y,z   from  the  direction  6,<j>.   That  element  of  flux  is 
given  by  {see  eq.  (2.24)  [5]} 

d<S>(x,y,z,6,<f>)   =  L(x,y,z,6,<J>)«cos0  «d<Wa  [W]  ,  (4.13) 

3. 

where  0   is  the  angle  between  the  ray  and  the  normal  to  da.   But  da  is  the  circular 
intersection  of  the  spherical  volume  element  with  a  plane  through  x,y,z  perpendicular  to 
the  ray,  so  that  6  =0  and  cos6  =  1  for  every  incident  ray  and  eq.  (4.13)  becomes 

3.  3 

d$(x,y,z,e,c(>)   =  L(x,y,z,9,(|))'diD'da  [W] .  (4.14) 

The  corresponding  element  of  fluence  rate  is  then 

dFt(x,y,z,6,<f>)   =   d$(x,y,z,e,<}>)/da 

=  L(x,y,z,0,<j))-dw  [W-m~2],  (4.15) 

and  the  total  fluence  rate  at  the  point  x,y,z  is  the  integral  over  all  such  elements  of 
incident  fluence  rate: 

Ft(x,y,z)   =   /  dFt(x,y,z,6,4) 

=  /  L(x,y,z,e,<J>)-du) 

0) 

=  /  /  L(x,y,z,e,<f>)«sin6«de'd<|>  [W-m  ],  (4.16) 

<j>  e 

where  the  limits  of  integration  for  6  and  <}>  include  all  rays  incident  on  the  elementary 
spherical  volume  at  x,y,z.   In  particular,  if  the  incident  radiation  comes  from  all 
surrounding  directions,  the  solid  angle  w  is  a  full  sphere  of  4tt  [sr]   and  this  becomes 

Ft(x,y,z)   =  /      dFt(x,y,z,6,4>) 
sphere 

=  /  L(x,y,z,0,<}>)  «d(o 
4tt 

=  /2lT/7rL(x,y,z,e,<j))'sine«de«d<f)  [W«m~2].  (4.16a) 

o  o 
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As  was  stated  earlier,  this  is  a  scalar  field  quantity,  since'  it  has  a  unique  value  at 
each  point  x,y,z   in  space.   By  contrast,  the  incident  radiant  flux  (surface)  density  or 
irradiance  is,  from  eq.  (A. 4a), 

W(x,y)   =  /  L(x,y,0,(|>)'dfi  =  /  L(x,y,6,<f>)  'cose-da)  [W«m  ].  (4.17) 

h  2tt 

The  latter,  as  we  have  said,  does  not  have  a  unique  value  at  a  point  x,y,z  in  space. 
Instead,  the  value  of  W(x,y)   in  eq.  (4.17)  depends  on  the  orientation  of  the  reference 
surface  at  the  point  x,y;   i.e.,  on  the  direction  of  the  normal  to  the  element  dA  at  x,y 
which  determines  the  value  of  6,   and  hence  of  cos6,   for  each  incident  ray.1 

The  integrals  in  eqs.  (4.16a)  and  (4.17)  differ  in  two  important  respects:   (1)  the 
first  does  not  contain  the  obliquity  factor  cosO,  which  appears  as  a  weighting  factor 
that  varies  with  surface  orientation  in  the  second;   (2)  the  first  integral  is  taken  over  a 
solid  angle  that  may  extend  to  the  full  sphere  of  4tt  [sr]  while  the  second  may  extend,  at 
maximum,  only  to  the  hemisphere  of  2ir  [sr]   on  just  one  side  of  the  reference-surface 
element  (of  the  tangent  plane  containing  the  element  dA) . 

Conversely,  we  can  write  expressions  for  radiance,  in  terms  of  radiant  flux  (surface) 
density  and  in  terms  of  fluence  rate,  from  eqs.  (4.7)  and  (4.15),  respectively.   They  are 

L(x,y,0,<j>)   =   <*W(*»7,M>    =   dW(x,y,M)   [w.B-2   ^  (    } 

att  cos6»dw 

and 

dF  (x,y,z,6,<j>)      __    . 
L(x,y,z,e,(f>)   =  - ^ [W-m^-sr].  (4.19) 

Note,  particularly,  that  eq.  (4.19)  is  an  expression  for  the  radiance  in  a  given  direction 
at  a  point  in  space  that  does  not  involve  any  reference  surface.  While  this  is  appealing 
to  the  theoretician,  we  don't  feel  that  it  is  of  strong  practical  significance  in  connec- 
tion with  optical  radiation  measurements  since  such  measurements  almost  invariably  involve 
an  instrument  with  a  defined  beam,  or  throughput,  involving  a  reference  surface  formed  by 
the  (imaginary)  plane  across  the  receiving  aperture  or  entrance  pupil  (and  the  solid  angle 
subtended  at  each  point  across  that  aperture  or  pupil  by  the  entrance  window,  the  field 


^ote  that,  in  the  rare  cases  where  we  directly  compare  the  incident  flux  (surface)  density 
W(x,y)   or  irradiance  E(x,y)   and  the  fluence  rate  F  (x,y,z)   at  the  same  point,  the 
coordinates  x,y  in  eq.  (4.17)  are  position  coordinates  on  the  reference  surface  while 
x,y,z  in  eq.  (4.16)  are  coordinates  in  space.   The  x,y  coordinates  of  the  common  point 
can  be  the  same  in  both  equations  only  if  the  reference  surface  is  a  plane  surface.   Then, 
by  appropriate  transformation  (translation  and  rotation)  the  origins  and  axes  of  both 
coordinate  systems  can  be  brought  into  coincidence. 
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solid  angle).1  Accordingly,  the  appropriate  definition  and  treatment  of  radiance,  for  our 
purposes,  is  in  terms  of  a  reference  surface,  as  in  eq.  (2.14)  [5]   {repeated  just  before 
eq.  (4.2),  above} . 

On  the  other  hand,  there  certainly  are  occasions  when  a  distribution  of  incident 
radiant  flux  with  respect  to  position  is  more  appropriately  or  meaningfully  given  in  terms 
of  fluence  rate  rather  than  irradiance  [incident  radiant  flux  (surface)  density] .  Accord- 
ingly, we  need  to  clearly  establish  the  relationship  between  these  two  quantities,  both 
given  in  [W'm   ].   We  substitute  E  for  W  in  eq.  (4.17),  since  we  are  concerned  only 
with  incident  radiation,  and  compare  it  directly  with  eq.  (4.16).   In  making  this  compari- 
son, we  recognize  that  the  ray-direction  coordinates  6,<|>  may  be  referred  to  arbitrarily 
oriented  fixed  axes  but  that  cosO   in  eq.  (4.17)  is  the  cosine  of  the  angle  between  the 
ray  direction  9,<f>   and  the  normal  to  the  receiving  surface  element  dA   upon  which  it  is 
incident.   Accordingly,  we'll  designate  the  polar  angle  from  the  normal  to  dA   as   6  , 
to  distinguish  it  from  the  direction  coordinate  9   referred  to  the  fixed  axes.   Then  it  is 
clear  that  the  element  of  fluence  rate  at  a  point  x,y  on  some  reference  surface  is 

dFt(x,y,e,<i>)  =     dE(x,y,e,<j>)/coser  =     dE^x.y.G,*)  [W-m~2]  ,  (4.20) 

where  dE   is  the  normal  irradiance  and  is  equal  to  the  irradiance  element  dE  only  when 
cos6  =  1,   i.e.,  when  9  =0   (when  the  ray  direction  0,<J>  coincides  with  the  normal  to 
dA  ) .   Thus  the  element  of  fluence  rate  at  a  point  is  equal  to  the  normal  irradiance  at 
that  point,  defined  as  the  irradiance  on  a  surface  element  at  that  point  divided  by  the 
cosine  of  the  angle  between  the  incident  ray  and  the  normal  to  the  surface  element.   Only 
the  elements   of  fluence  rate  and  irradiance,  associated  with  individual  rays  (elements  of 
throughput)  can  be  directly  related  in  this  way.   The  integrated  total  fluence  rate  and 
irradiance  at  a  point  cannot  be  directly  converted  from  one  to  the  other.   If  values  for 
both  are  needed,  they  are  most  easily  obtained  from  the  incident  radiance  distribution  and 
eqs.  (4.16a)  and  (4.17).   But  this  seems  unlikely;  ordinarily  only  one  or  the  other  will  be 
useful  or  significant. 

Discussion  of  illustrative  examples  involving  radiant  flux  (surfrce)  density  and 
radiant  fluence  rate  is  deferred  until  the  concept  of  radiant  intensity  has  been  presented. 
Then  the  examples  will  help  to  clarify  all  three  of  these  concepts  and  their  interrelation- 
ships. 

DIRECTIONAL  DISTRIBUTION:   RADIANT  INTENSITY— TRADITIONAL  APPROACH. 

Consider  the  beam  of  radiation  flowing  from  all  points  on  the  surface  of  an  extended  source 
S   to  a  small  receiving  surface  of  area  AA,   as  shown  in  figure  4.4.   Designate  the  mag- 
nitude of  the  total  flux,  reaching  the  receiver  from  S,   as  A$.   Let  the  solid  angle 
subtended  at  any  point  P  on  the  source  surface  by  the  receiver  be  designated  as  Au). 


^The  theory  of  stops  (pupils  and  windows)  will  be  found  in  any  standard  text  on  geometrical 
optics  [13,14]. 
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Figure  4.4.     Experiment  to  develop  the  concept  of  (radiant)  intensity 
—  traditional  approach. 


If  the  receiver  is  moved  arbitrarily  far  from  the  source,  both  A$  and  Aw  decrease  to 
the  vanishing  point,  but  their  quotient  A$/Aw  approaches  a  constant  value  called  the 
radiant  intensity  I  of  the  source  in  the  direction  of  the  receiver: 


I   E 


d3>   rTT   -1, 

"ST  [w*sr  ]' 


(4.21) 


Also,  as  the  receiver  moves  away,  the  small  solid  angle  that  it  subtends  at  every  other 
point  on  the  source  surface  becomes  more  nearly  equal  to  and  more  nearly  parallel  to  Aco. 
Ultimately,  the  rays  from  all  points  of  S  to  AA  "at  infinity"  can  be  thought  of  as 
forming  a  set  of  parallel  solid-angle  elements  dcu,   each  containing  a  ray  from  that  source 
point  in  the  direction  of  AA,   as  depicted  in  figure  4.5. 


Figure  4.5.    Beam  from  the  source  of  figure  4.4  to  a  receiver 
"at  infinity." 
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If  L  is  the  radiance,  in  the  direction  of  the  receiver  AA,   of  an  element  of  source 

surface  dA  ,   then,  by  eq.  (2.15  [5],   the  element  of  flux  propagated  from  dA   to  AA 
s  s 

is  L'dA  •cos9*do),  where  0  is  the  angle  between  the  normal  to  dA   and  the  ray  from  dA 
s  s  s 

to  AA.   Then  the  flux  in  the  solid-angle  element  dw,   propagated  from  all  of   S   to  AA, 
is  given  by  the  integral 

d$  =  /  L»dA  'Cos6*dw  =  dw/  L*cos9«dA  ,  (4.22) 

S    s  S       s 

so  that 

I  e  -P-  =   f  L.cos6«dA  [W-sr-1].  (4.23) 

dw      s        s 

This  gives  us  an  alternative  expression  for  the  radiant  intensity  of  an  extended  source  in 
a  given  direction  in  terms  of  the  radiance  of  each  point  (surface  element)  of  the  source  in 
the  given  direction. 

Measurements  of  the  incident  flux,  per  unit  solid  angle  subtended  at  the  source, 
reaching  a  detector  from  a  source  at  a  finite  distance  are  only  approximations  to  the 
defined  quantity  given  by  eq.  (4.21)  or  (4.23) 

I  m    -£*-  [W-sr-1].  (4.24) 

We'll  examine  the  factors  that  contribute  to  the  degree  of  approximation  after  we've 
defined  radiant  intensity  a  bit  more  carefully  and  have  established  its  relationships  to 
other  radiometric  quantities  in  a  somewhat  different  manner  and  in  greater  detail. 

The  definition  or  description  of  the  concept  of  radiant  intensity  seems  to  give  rise 
to  more  problems  than  that  of  any  other  radiometric  quantity.   Some  readers  are  puzzled  by 
the  fact  that  we  have  related  the  flux  A$,   coming  from  all  parts  of  the  source  in  figure 
4.4  (not  just  that  flowing  in  the  solid  angle  Aw  from  the  particular  point  P  shown  in 
the  illustration) ,  to  the  more  limited  solid  angle  Aw  in  the  quotient  that  converges  to 
the  limiting  value  of  eq.  (4.21)  and  that  appears  in  eq.  (4.24).   The  point  P   in  figure 
4.4  is  any   point  on  the  emitting  surface  of  the  source.   There  is  (not  shown)  a  similar 
solid  angle,  approximately  equal  to  Aw,   subtended  at  every  other  source  point  from  which 
source  radiation  reaches  the  receiver.  We  use  the  value  Aw,   then,  as  a  representative 
approximation  to  all  of  these  solid  angles,  subtended  at  all  points  of  the  source.  As  the 
distance  between  source  and  receiver  increases,  the  approximation  improves  until,  in  the 
limit,  as  suggested  by  figure  4.5,  all  of  the  elementary  solid  angles  become  parallel  and 
equal  to  dw.   In  a  very  large  number  of  cases,  the  approximation  becomes  practically 
satisfactory  (within  useful  limits  of  accuracy)  within  distances  available  in  the  labora- 
tory.  As  we'll  see,  we  soon  reach  a  point,  with  most  sources,  where  the  flux  A$  and  the 
solid  angle  Aw  both  become  inversely  proportional  to  the  square  of  the  distance  between 
convenient  points  on  both  source  and  receiver  (the  size  of  which  may  also  be  appreciable) , 
so  the  quotient  becomes  and  remains  constant  as  that  distance  is  further  increased.   Then 
we  have  a  useful  value  of  source  radiant  intensity  that  characterizes  the  source  emission 
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in  the  given  direction.   It  is  often  overlooked,  however,  that  this  depends  on  the  fact 
that  the  radiance  of  a  point  on  the  surface  of  a  source  usually  doesn't  change  abruptly 
with  a  small  change  in  direction.  We'll  examine  this  point  in  more  detail  later,  in  the 
text  accompanying  figure  4.12,  below. 

DIRECTIONAL  DISTRIBUTION:   RADIANT  INTENSITY.    For  a  more  detailed  treatment  of  radiant 
intensity,  consider  again  the  experiment,  first  introduced  in  Chapter  2  [5],  with  an  iso- 
radiance  visible  source,  two  beam-forming  apertures  in  opaque  screens,  and  a  photocell, 
arranged  as  illustrated  in  figures  4.6  and  4.7(a)  (figures  2.2  and  2.3  of  Chapter  2  [5]). 
A  uniform,  clear  propagation  path  is  assumed,  with  no  variation  in  refractive  index  and  no 
attenuation  (by  absorption  or  scattering).   Previously,  we  found  that,  as  we  varied  the 
areas  AA..   and  AA„,   and  the  respective  tilts  6..   and  6„,   of  the  apertures  as  well  as 
the  distance  D  between  them,  the  radiance  L,  was  given  by  eq.  (2.4).   That  equation  is 
repeated  here,  but  as  an  approximation 


L  » 


A$>D 


AA..  •cos01  *AA  'COS0- 


(4.25) 


because  we  also  found,  subsequently,  that  the  exact  expressions  for  L  were  given  in  eqs. 
(2.5)  and  (2.13)  as 


Source 


Screen 


Screen 


Figure  4.6.   Experiment  to  develop  the  concept  of  radiance. 

(This  figure  appeared  previously  as  figure  2.2  [5].) 
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Lim  A$ 

L  =  AA„->0 


AaJ-X)   AA^cose^CAA^cosG^D2) 

=  aT+O  — ^—A =    ..   d2?   . [W-m-^sr"1].  (4.26) 

1     AA1  •cos61 'Aa)1?       dA- •cos01  •dw1  „ 

Aa)12"*° 

These  equations  define  the  radiance  L  that  is  geometrically  invariant  along  any  ray  in  a 

2 
lossless,  isotropic  (iso-index-of-refraction)  medium.   Here,  the  quantity  duo..  „  =  dA„#cos6„/D 

is  the  element  of  solid  angle  subtended  by  dA„  at  dA.   {see  eq.  (A2-13) ,  Appendix  2  [5]}. 

The  first  aperture  screen  and  the  source,  taken  together,  are  equivalent  to  an  iso- 

radiance  source,  of  area  AAn   filling  the  aperture  in  the  screen,  because  every  ray  of 

the  beam,  through  every  point  across  that  aperture,  has  the  same  radiance  L.  Whenever  it 

will  simplify  our  discussion,  then,  we  may  replace  the  source-screen  combination  by  a  bare 

iso-radiance  plane  source  of  radiance  L,  of  area  AA  =  AA..   tilted  at  an  angle  6  =  6,  , 

si  si 

at  that  same  location.   Similarly,  the  combination  of  the  second  aperture  screen  and  the 
photocell,  taken  together,  is  equivalent  to,  and  may  be  replaced  by,  a  bare  plane  receiver- 
detector  of  area  AA  =  AA_,   tilted  at  an  angle  G  =  6„,   centered  at  a  distance  D  from 
r     2  r    2 

the  center  of  the  source,  and  with  uniform  isotropic  response  (equal  response  to  flux 
incident  at  any  point  and  from  any  direction  within  the  beam) .   These  equivalent  configura- 
tions, that  may  be  interchanged  at  will  without  altering  the  propagation  of  radiant  energy 
in  the  beam,  are  illustrated  in  figures  4.7(a)  and  (b) . 

In  addition  to  rewriting  eq.  (4.26)  in  terms  of  the  quantities  in  figure  4.7(b),  it 
will  be  convenient  to  rearrange  it  somewhat  to  examine  what  happens  as  we  go  to  the  limits 
in  different  ways. 

Llm  A$ 

L  =  AA  -K) 


s     AA  «cos6  •  Aaj 
a    n         s     s   sr 
Aw  -K) 

sr 


Llm            A$  -2   -1 

=  AA  -*0  — = —  [W-m   -sr   ],   or  (4.27) 

_  S    AA  -cose  -(AA  «cos6  /D  ) 
D  ->-  °°    s     sr     r 

Lim 

=  AAg->0  — —  [W-m"  -sr   ].  (4.28) 

AAS  _   AA  'CosG  «(AA  -cosO  /D  ) 
AA  -K)    s     s    r     r 
r 

Note  that  the  solid  angle  Aw    (subtended  at  AA   by  AA  )   may  be  reduced  either  by 
increasing  D,   as  in  eq.  (4.27),  or  by  reducing  the  receiver  area  AA  ,   as  in  eq.  (4.28), 
or,  of  course,  by  doing  both.  We  examine  the  results  of  doing  one  or  the  other  alone, 
first  looking  at  the  effects  of  increasing  the  separation  distance  D  as  in  eq.  (4.27) 

24 


Figure  4.7(a).   Tilted  Apertures. 

(This  figure  appeared  previously  as  figure  2.3  [5].) 


ftr-if— 


Figure  4.7(b).   Bare  source  and  detector  equivalent  to  figure  4.7(a) 
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while  holding  the  source  area  AA   fixed.   Then,  holding  D  constant,  we'll  examine  the 
effects  of  independently  varying  AA  .   A  similar  analysis  of  eq.  (4.28)  will  follow  in  the 
next  section. 

As  the  separation  D  is  increased,  both  the  flux  A$   and  the  solid  angle  Aw    are 
reduced.   However,  if  the  small  source  AA   is  strong  enough  and  the  small  detector  AA 
is  sufficiently  sensitive,  a  point  is  soon  reached  where  the  quotient  A$/Aoj    becomes  and 
remains  constant  (within  achievable  measurement  accuracy)  with  further  increase  of  D 
until  there  is  no  longer  enough  flux  A$  to  produce  a  measurable  response.   The  limiting 
value  of  the  quotient   A$/Aw    is  called  the  radiant  intensity  AI  of  the  source  AA   in 
the  direction  of  the  detector  (in  a  direction  making  an  angle  6   with  the  normal  to  the 
plane  source  surface) : 

AT     Lim   AJ> d$    r.T   -1,  ,.    OQ. 

AI  =    ^  a  y~     =  [W-sr   ].  (4.29) 

AA  •cos0  /D        dw 
r     r  sr 

This  can  be  repeated  for  different  values  of  source  area  AA   with  similar  results  except 
that,  as  AA   is  reduced,   AI   is  correspondingly  reduced.   The  maximum  distance  D  where 
there  is  still  enough  flux  A$   to  produce  a  measurable  response  is  also  reduced  with 
decreasing  AA  .   Accordingly,  we  can  substitute  from  eq.  (4.29)  in  eq.  (4.27)  to  obtain 

T     Lim       AI dl rTT  -2   -1,  ,.  oriN 

L     a  a  ^n  ~~n 5 —  =  ~n s —  [w*m  *sr  ]»  (4.30) 

AA  -K)   AA  -cose        dA  «cos6 
s      s     s         s     s 

from  which  we  have 

dl  =  L«cos0  -dA  [W«sr_1].  (4.31) 

Up  to  this  point,  for  simplicity,  we  have  postulated  an  iso-radiance  source.   More 
generally,  however,  the  limiting  values  of  radiant  intensity  can  be  obtained  in  exactly  the 
same  way  even  when  the  source  radiance  varies  with  position  and/or  direction.   It  is  only 
necessary  to  go  to  smaller  areas  and  solid  angles  to  obtain  the  constant  values,  corre- 
sponding to  the  limiting  values  that  are  the  radiant  intensities.  Accordingly,  we  will 
restate  eq.  (4.31)  more  completely  and  explicitly:   The  element  of  radiant  intensity 

dl(x,y,9,<f>) ,   in  the  direction  0,<f>,   of  a  source-surface  element  dA   of  radiance 

s 

L(x,y,e,(j))   at  the  point  x,y  and  in  that  direction,  is  given  by 

dl(x,y,e,<|>)  -  L(x,y,e,<J>)»cose  «dA  [W«sr-1].  (4.32) 

Then,  by  integrating  both  sides  of  eq.  (4.32),  we  obtain  an  alternative  defining  equation 
for  the  radiant  intensity  of  an  extended  source  of  radiating  (emitting  and/or  reflecting  or 
scattering)  area  A  : 
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1(9,*)   =  /  dl(x,y,e,<f>) 


s 


-1, 


=  /  L(x,y,e,<t>)'cos6  «dA   [W«sr~  ]. 

A  S   S 

s 


(A. 33) 


Note  that,  in  these  equations,  the  angle 

be  equal  to  the  ray-direction-coordinate  angle  0  only  when  the  element  dA 


(between  a  ray  and  the  normal  to  dA  )  will 


is  oriented 


with  its  normal  parallel  to  the  polar  axis  of  the  direction  coordinates.   Generally,  for  a 
non-planar  reference  surface,  the  angle  6   for  parallel  rays  will  vary  from  point  to 
point  across  that  reference  surface,  as  illustrated  in  figure  4.8. 


^r  L(x,,y,,0,<£) 


L(x2,y2  ,0,<£) 


Ux3,y3,0,</>) 


Figure  4.8.    Illustration  to  clarify  the  designations 


and   6 


Given  a  fixed  direction   6,<f>   and  a  curved  reference 

surface,  the  angle   6   between  the  direction  6,<j> 

and  the  normal  to  a  surface  element   dA   varies  from 

s 

point  to  point  across  the  curved  surface. 
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POINT  SOURCES  and  RECEIVERS  and  the  INVERSE- SQUARE  LAW.   Although  eq.  (4.33)  explicitly 
defines  a  quantity  called  the  "radiant  intensity"  for  any  source  of  extended  area  A    it 
turns  out  that  this  quantity  is  useful,  in  general,  only  when  the  source  can  be  regarded  or 
treated  as  a  "point  source."  The  classic  example  of  what  is  meant  by  "treated  as  a  'point 
source'"  is  a  star,  even  the  nearest  of  which,  in  spite  of  being,  in  many  cases,  incredibly 
large  in  comparison  with  the  earth,  or  even  the  sun,  is  so  far  distant  that  our  telescopes 
can't  distinguish  or  "resolve"  different  points  on  its  surface.  Accordingly,  for  all 
practical  purposes,  the  rays  reaching  the  earth  from  all  parts  of  the  emitting  surface  of  a 
star  are  those  that  travel  from  it  in  parallel  lines  in  a  single  direction  {the  direction 
6,<)>  in  eq.  (4.33)}.   In  a  great  many  cases,  when  shorter  distances  are  involved,  a  criterion 
that  determines  whether  or  not  we  can  treat  an  extended  source  as  a  point  source,  with  a 
useful  value  of  radiant  intensity,  is  whether  or  not  the  irradiance  or  radiant  fluence  rate 
at  a  point  of  interest,  due  to  the  source  in  question,  obeys  the  inverse-square  law  within 
a  given  tolerance.  Accordingly,  we  next  want  to  look  more  closely  at  the  relationship 
between  radiant  intensity  and  irradiance  or  radiant  fluence  rate  known  as  the  "inverse- 
square  law." 

We  can  obtain  the  inverse-square  law  by  going  back  to 

Lim   A$    ,..  -2   -1,  ,.  OON 

L  =  AA  -»0  2 —  [W'm   'sr   ]  (4.28) 

s    AA  -cose  '(AA  «cos0  /D  ) 
a  a   r.    s     s    r     r 
AA  -*0 

r 

to  see  what  happens  as  we  first  reduce  AA   while  holding  AA   constant.  As  AA 
becomes  smaller,  the  flux  A$  reaching  the  detector  from  the  source  is  correspondingly 
reduced,  but  the  quotient  A$/AA   becomes  and  remains  constant,  and  we  recognize  [from  eq. 
(4.1)]  that,  in  the  limit,  this  quotient  is  the  incident  radiant  flux  (surface)  density  or 
irradiance  AE.   We  designate  it  as  AE,   rather  than  just  E,  because  it  is  also  depend- 
ent on  AA  .   Accordingly,  eq.  (4.28)  becomes 
s 

T     Lim        AE'D2  dE'D2       r   -2   -11 

AA  -K)  AA  -cose  -cos0        dA  -cose  «cos6    L     *Sr   J'         IA.J4J 
s      s     s     r        s     s     r 

When  this  is  rearranged  and  compared  with  eq.  (4.31),  we  have 

2 

L-dA  -cosO   =  dl  =   dE'°   [W-sr-1]  (4.35) 

s     s  cos6 

r 

from  which  we  obtain  the  inverse-square  law,  giving  the  element  of  normal  irradiance  dE 

n 

or  radiant  fluence  rate  dF   [see  eq.  (4.20)]  at  a  distance  D  from  the  source  element 
dA   in  a  direction  in  which  its  intensity  is  dl: 

dER  =  dF   =  dE/cosO   =  dl/D2  =  L«cos9  «dA  /D2  [W«nf2],  (4.36) 

or,  more  explicitly, 
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dEn(xr,yr,zr,e,<t>)      =     dFt(xr,yr,zr,e,(t»)      =     dE(xr,yr,zr,6  ,<j>)  /cos6r 

=     dI(xs,ys,zs,e,<t>)/D2 

2  -2 

=     L(x    ,y    ,z    ,e,<t>)*cos6s-dAs/D      [W-m      ], 


(4.36a) 


where   (see  figure  4.9) 


^8>y»  »^s 


^r»)fr  t  Zr 


Figure  4.9.    A  section  through  source  and  receiver  surfaces  of  arbitrary  shape. 

A  source-surface  element   dA   at  x  ,y  ,z   and  a  receiver- 

s       s   s   s 

surface  element  dA   at  x  ,y  ,z   are  separated  by  a  distance  D. 
For  purposes  of  illustration,  the  normals  to  the  two  elements  are 
shown  as  lying  in  the  same  plane  although  they  are  not  actually  so 
constrained. 
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D2  =   (xr-  xs)2  +  (yr-  ys)2  +  (z^  zg)2  [m2] 

is  the  square  of  the  (slant)  distance  between  the  source  element  at  x  ,y  ,z   and  the 
receiver  element  at  x  ,y  ,z  , 

6  =  tan"1  {[(xr-  xs)2  +  (yr-  y,,)2^/^-  zg)  }  [rad]  , 

<j>  =  tan"  {(yr-  y  )/(x  -  xg)  }  [rad], 

6   [rad]   is  the  angle  between  the  direction  8,(f>  and  the 

normal  to  the  element  of  source  emitting  surface  dA 

at  x  ,y  ,z  ,   and 
s  J  s     s 

G    [rad]   is  the  angle  between  the  direction  0,<()  and  the 

normal  to  the  element  of  receiving  surface  dA   at 

x  ,y  ,z  . 
r  yr  r 

The  inverse- square  law  is  a  relationship  for  a  single  ray  or  throughput  element  through  a 
source-surface  element   dA  ,   of  intensity  dl  =  L«cos6  *dA  ,   in  the  direction  6,<J> 
from  the  point  x  ,y  ,z  .   It  gives  the  normal  irradiance  dE   [=  dE/cos6  ]   or  radiant 
fluence  rate  dF   at  a  distance  D  in  the  direction  9,<j>  from  that  source  element. 
Integrals  based  on  these  relationships  can  successfully  account  for  the  observed  values  of 
radiation  flux  in  beams  involving  extended  source  and  receiver  surfaces.  Nevertheless, 
although  a  "point  source"  is  often  specified,  it  is  not  unusual  to  find  less  precise  state- 
ments to  the  effect  that  the  flux  received  from  a  distant  source  varies  inversely  with  the 
square  of  that  distance.   Such  a  statement  can  only  be  approximately  true,  in  general,  when 
applied  to  real  (extended)  sources  and  receivers.   The  inverse-square  law  [eq.  (4.36)],  as 
we  have  seen,  applies  exactly  to  each  source-element — receiver-element  pair  and  the  direct 
or  slant  distance  between  them.   That  distance,  however,  will,  in  general,  be  different, 
perhaps  only  slightly  different  but  still  different,  for  different  pairs.   Thus,  an  overall 
inverse-square  relationship  based  on  just  a  single  (average  or  nominal)  source-receiver 
distance  can,  in  general,  only  be  an  approximation.   Next,  we'll  take  a  closer  look  at  such 
approximations  to  the  inverse- square  law  and  try  to  evaluate  the  extent  of  their  usefulness. 
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INVERSE- SQUARE- LAW  APPROXIMATIONS  for  EXTENDED  SOURCES  and  RECEIVERS.    There  are  three 
main  considerations  that,  together,  determine  when,  in  a  given  measurement  configuration, 
it  is  useful  to  treat  a  source  as  a  point  source,  characterized  by  its  radiant  intensity 
{eq.  (4.33)},  or  as  an  extended  source,  characterized  by  its  (average)  radiance: 

(1)  whether  or  not  the  source  fills  the  entrance  window1  or  angular 
field  of  view  of  the  radiometer;2 

(2)  the  size  of  the  maximum  transverse  dimensions  of  the  source,  and  of 
the  instrument  receiving  aperture  or  entrance  pupil,1  relative  to 
the  distance  D  between  them;  and 

(3)  the  directional  distribution  of  radiance  from  each  surface  element 

dA   across  the  radiating  surface  of  the  source  A  ,  and  particularly 
the  effects  of  obscuration  (when  an  opaque  portion  of  an  irregularly 
shaped  source  may  block  rays  from  other  parts  of  the  emitting 
surface) . 

So  far,  we  have  considered  only  fairly  regular  well-behaved  sources,  most  of  them  plane 
surfaces  or  equivalent  to  such  a  plane  emitting  surface,  but  when  we  start  dealing  with 
real  sources,  such  as  the  very  common  tungsten  lamps  with  coiled  filaments,  severe  problems 
can  arise  due  to  the  fact  that  parts  of  the  emitting  filament  lie  behind  other  opaque 
portions  of  the  coil.   We'll  discuss  this  more  a  bit  later.   In  addition,  it  should  be 
noted  that  the  problem  can  also  be  complicated  by  directional  variations  in  instrument 
response.   But  we'll  leave  it  with  just  that  caution  and  assume  now  that  we  have  good 
instruments  that  respond  equally  to  incident  flux  from  any  direction  within  the  angular 
field,  keeping  in  mind  that  things  may  be  quite  different  if  this  important  assumption  is 
not  valid. 

Although  the  three  considerations  are  listed  above  in  the  order  in  which  they  are 
usually  checked,  we'll  start  with  the  second,  the  transverse  dimensions  of  source  and/or 
receiver  in  relation  to  their  nominal  separation  distance  D.   The  geometrical  relations 
involved  can  be  seen  more  clearly  if  we  simplify  our  experimental  configuration  to  one 
often  employed  in  photometry,  observing  a  small  source  with  a  small  detector  on  an  optical 
bench  where  the  distance  between  them  is  easily  changed  and  easily  measured.   In  figure 
4.10  we  see  an  iso-radiance  source  consisting  of  a  tungsten  lamp  in  an  opaque  housing  with 
an  opal-glass  diffuser  across  one  side  to  produce  an  iso-radiance  radiating  surface.   An 
opaque  screen  with  a  small  circular  hole  of  area  AA  ,   perpendicular  to  the  bench  so  that 
6  =0   (cosO  =  1),   limits  the  beam  of  rays  from  source  to  receiver,  establishing  the 


1See  any  standard  text  on  geometrical  optics  [13,14],   A  future  chapter  is  planned  dealing 
with  stops,  pupils,  windows,  and  baffles  and  their  role  in  radiometry. 

In  a  scanning  system,  we  are  concerned  with  filling  the  instantaneous  entrance  window  or 
angular  field  of  view;  in  an  imaging  system  we  are  concerned  with  filling  the  instantaneous 
entrance  window  or  angular  field  of  view  of  a  resolution  element. 

31 


1 

1 
cu 

TJ 

OJ 

CU 

u 

Tl 

a) 

a 

> 

cu 

c 

u 

•H 

X 

a> 

cu 

U-l 

4-1 

o 

o 

4-1 

4-1 

a 

cu 

CO 

o 

c 

C 

o 

o 

•H 

o 

4J 

cu 

E 

x: 

■H 

j-i 

O 

•4-1 

M 

O 

a 
a. 

en 

ca 

CO 

cu 

rH 

C 

cfl 

rH 

o 

3 

•H 

4-t 

■U 

cu 

CJ 

en 

cd 

3 

M-4 

O 

<4H 

o 

CU 

on 

Tl 

d 

C 

Cfl 

cfl 

j-j 

>, 

. 

cu 

4-1 

C 

Ji 

■H 

o 

j-i 

en 

•H 

C 

4-1 

a. 

cu 

cfl 

o 

ij 

•H 

H 

a 

T3 

cu 

•H 

ctt 

> 

M 

cu 

y^ 

!-i 

T3 

4-1 

c 

■H 

O 

cfl 

14-1 

u 

•H 
Tl 

o 

4-1 

Cfl 

13 

c 

S-i 

cfl 

cu 

*■ — ' 

r- 

B 

•H 

cu 

cu 

^ 

CJ 

u 

CU 

1-1 

cfl 

Pu 

3 

3 

X 

o 

V 

w 

CO 

cn 

cu 

3 
00 

•H 


32 


position,  orientation,  and  area  of  the  effective  source  AA  ,   as  before  in  figure  4.7(b). 
All  of  these  source  components  are  carried,  and  move  together  as  a  unit,  on  a  singl  optical- 
bench  carriage  with  an  index  or  cursor  mark  aligned  with  the  vertical  plane  of  the  screen 
(the  plane  of  the  source  AA  ) .   A  second  carriage  carries  the  receiver  components,  which 
consist  of  another  vertical  screen  with  a  small  aperture  of  area  AA   (6  =0;  cos6  =  1), 
followed  by  a  photocell  with  a  sensitive  area  large  enough  to  be  sure  to  accept  all  of  the 
rays  of  the  beam  through  this  second  aperture.   Again,  the  cursor  mark  or  index  on  the 
receiver  carriage  is  aligned  with  the  vertical  plane  of  the  aperture  screen  so  the  distance 

D  between  the  source  AA   and  the  receiver  AA   is  easily  measured  by  the  difference 

s  r 

between  the  cursor-index  readings  for  the  two  carriages  on  the  horizontal  scale  of  the 
optical  bench. 

The  value  of  the  flux  6$  in  the  beam  is  measured  by  the  photocell.   Its  approximate 
value  in  terms  of  beam  parameters  is  given  by 

6$  w   (L-AA  /D2)-AA   «   (AI  /D2) «AA   ~  AE   «AA   [W] ,  (A. 37) 

s       r        s       r       rs    r 

since,  from  eq.  (4.31),  a  plane  iso-radiance   (L  =  a  constant)   source  of  area  AA   has  a 

radiant  intensity  AI  «=>  L«cos6  «AA   and  since  cos0  =  cos6  =  1*.   For  a  still  larger 
s        s   s  s       r  ° 

source  and  receiver,  of  area  A   and  A  ,   respectively,  this  same  approximation  exists 
and  we  want,  now,  to  determine  the  degree  of  approximation  involved  in 

$  m  E-A   »  (I/D2)«A   «  L«(A  «A  /D2)  [W] ,  (4.38) 

or,  putting  it  another  way,  we  want  to  evaluate  the  discrepancy 

A$  E  $  -  (L'A  «A  /D2)  [W]  (4.39) 

in  order  to  assess  the  percentage  or  relative  inaccuracy  A$/$  as  a  function  of  the  trans- 
verse dimensions  of  A   and  A   relative  to  D. 

s        r 

In  order  to  evaluate  eq.  (4.39),  we  need  an  exact  expression  for  the  flux  $  in  the 
beam  between  the  extended  source  and  the  extended  receiver.  But  we  already  have  such  an 
expression  in  eq.  (2.12)  [5]: 


$ 


-2 


=  L«J  /  D  -cose  -cos6  «dA  «dA   =  L»0  [W] ,  (4.40) 


A  A 
r  s 


where  D   is  the  slant  distance  between  each  pair  of  surface  elements  dA   and  dA  ,  not 
s  r  s        r 

the  nominal  separation  or  the  perpendicular  distance  D  between  the  vertical  planes  of  A 
and  A  ,   as  read  on  the  optical  bench.   The  throughput  0  between  two  parallel  coaxial 
circular  plane  areas  A   and  A   separated  by  a  distance  D  along  the  common  axis  through 
their  centers  is  available  in  Appendix  3.   If  the  radius  of  A   is  r  ,   and  of  A   is 
r  ,  we  can  write  [see  eq.  (A3-15)] 

0  =  (Tr^/2).{(r^  +  r'  +  D")  -  [(r*  +  r'  +  dV  -  4r^r^i  [m^-sr]. 
s    r  s    r  s  r 
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We  can  also  combine  eqs.  (4.39)  and  (4.40)  to  simplify  the  expression  for  the  relative 
discrepancy  A$/$. 

A$  =  L'0  -  L»(A  »A  /D2)   =  L-A6  [W] ,  (4.42) 

where 

A6  =     G  -  (A  »A  /D2)  [m2«sr]  (4.43) 

from  which  it  is  clear  that 

A$/$  =  A0/0  [dimensionless]  (4.44) 

is  the  relative  discrepancy. 

There  is  a  common  "rule  of  thumb"  in  radiometry  (photometry)  that  application  of  the 
approximate  inverse- square  relationship 

0  «  A  -A  /D2  [m2«sr]  (4.45) 

s  r 

will  be  adequate  for  accuracies  of  about  99  per  cent  when  the  separation  D  is  at  least  10 

times  the  maximum  transverse  dimension  of  A   or  A  .   For  example,  in  the  present  case, 

2  2  s       r 

where  A  =  ir-r   and  A  =  ttt  ,   those  transverse  dimensions  are  2r   and  2r  ,  respec- 
s      s        r      r  s        r 

tively,  so  we  want  to  evaluate  the  inaccuracy  A0  when  D  s  10' (2r  )   and  D  S  10* (2r  ). 
To  make  it  simple,  let   r  =  r   and  D  =  20  r  =  20  r  ,   and  solve  eqs.  (4.43)  and  (4.41) 
to  obtain  the  relative  inaccuracy  A0/0.   The  result  is  just  under  0.005  or  0.5%.  But, 
to  get  this,  we  have  assumed  ideal  conditions,  an  iso-radiance  source  and  a  receiver- 
detector  that  responds  uniformly  and  isotropically  to  flux  incident  at  any  point  across 

A   from  any  direction  within  the  beam  from  A  . 
r  s 

For  a  somewhat  better  idea  of  the  applicability  of  this  "rule  of  thumb"  to  real  con- 
figurations, consider  the  extreme  possibility  that  the  coaxial-parallel-disc  source  and 
receiver  are  both  so  non-uniform  that  all  of  the  radiation  emitted  isotropically  by  the 
source  disc  actually  comes  only  from  a  small  region  about  a  point  S   at  the  edge  of  the 
disc  while  all  of  the  non-uniform  isotropic  responsivity  of  the  receiver  is  similarly 
concentrated  in  a  small  sensitive  region  about  a  point  R  at  the  edge  of  that  disc, 

diametrically  opposite  to  point  S  on  the  source  disc,  as  depicted  in  figure  4.11.  The 

2     2    5s 
actual  slant  distance  separating  these  points  is  then  D  =  2{r  +  (D  /4) } 2,   where 

D  =  20r.   Then  the  relative  or  fractional  error  resulting  from  the  use  of  D  rather  than 

2    2   2 
D   in  the  inverse-square  relationship  is   (D  -  D  )/D  =  0.01  or  1%.   Thus,  the  "rule  of 
s  s 

thumb"  still  holds  for  an  estimated  accuracy  of  99%,   even  with  such  extreme  non-uniform- 
ity as  long  as  both  source  and  receiver  are  still  isotropic  at  each  point  of  A   and  A  . 

This  brings  us  to  the  third  consideration,  the  directional  distribution  of  radiation 
[and  of  receiver-detector  response] .   Although  the  isotropic  assumption  is  not  realistic 
over  a  full  hemisphere  for  most  surfaces,  it  is  often  a  good  approximation  over  a  fairly 
wide  range  of  directions  about  the  normal.   Accordingly,  our  "rule  of  thumb"  approach  is  a 
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Figure  4.11.    Extremely  non-uniform  isotropic  source  and  receiver  to  assess  a 
"worst-case"  approximation  to  the  inverse-square  law. 

Source  and  receiver  are  equal,  parallel,  co-axial,  circular 
discs  of  radius   r   separated  by  distance  D.   The  slant 
distance  Dg  between  extreme  points   S   and   R   is  seen  to  be 

0  9  J- 

D   =   (4r  +  D  )  2. 

s  ' 

2  0 

Evaluated  for  D  =  20r,    D   =   1.01-D 

s 


reasonable  one  for  most  sources  and  receivers  with  reasonably  flat  emitting  and  receiving 
surfaces  that  are  oriented  perpendicular  to  the  axis  between  their  centers,  as  in  the 
illustrative  experimental  configurations  we've  been  discussing.  We've  also  seen,  in  figure 
2.10  of  Chapter  2  [5],  that,  with  truly  iso-radiance  emitting  surfaces,  the  shape  is 
immaterial  as  long  as  it  subtends  the  same  solid  angle  at  the  receiver.   However,  with  real 
emitting  surfaces  the  radiance  may  change  when  the  direction  to  a  receiver  from  an  off-axis 
point  changes  slightly  as  the  receiver  moves  farther  away.   For  example,  if  the  source  in 
figure  4.11  is  not  isotropic,  we  can't  base  our  conclusion  on  just  the  geometrical  rela- 
tions but  must  take  into  account  that  the  radiance  L  of  the  rays  from  S  to  R  is  also 
a  function  of  the  separation  distance  D.   The  situation  becomes  even  worse  with  sources  of 
such  irregular  shape  that  rays  from  some  parts  of  the  emitting  surface  can  be  blocked  in 
some  directions  by  other  parts  of  the  source  itself  while  other  rays  in  directions  very 
close  to  them  are  not  blocked,  thus  producing  violent  variations  in  the  directional  dis- 
tribution of  source  radiance  (see  figure  4.12).  Not  only  does  this  alter  the  directional 
distribution  of  radiance;  it  may  also  change  the  effective  emitting  area  A   as  a  function 
of  direction.  With  the  regularly  periodic  structure  of  a  helical  lamp  filament  (see  figure 
4.12),  one  would  expect  that  from  some  distant  observation  points,  most  of  the  inner  sides 
of  the  far  loops  of  the  coil  would  be  hidden  behind  the  near  loops  while  a  small  rotation 
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of  the  lamp  would  bring  them  all  into  view  in  the  spaces  between  the  near  loops;  and  that 
this  might  be  repeated  as  the  lamp  is  rotated  further.   Just  such  a  periodic  variation  (of 
about  one  per  cent)  in  received  flux  with  lamp  rotation  has  been  observed  (see  figure  18 
and  accompanying  text  in  [15]).   This  situation  is  too  complex  for  a  precise  analysis 
except  to  point  out  that,  for  such  sources,  there  may  not  be  a  sufficiently  accurate 
inverse-square  relationship  such  as  eq.  (4.38),  even  though  the  inverse-square  law  of  eq. 
(4.36a)  is  still  perfectly  valid  between  every  pair  of  source-  and  receiver-surface  elements. 
The  overall  radiant  intensity  of  the  entire  source s    in  such  cases,  may  be  too  variable  to 
be  a  useful  quantity  by  which  to  characterize  the  lamp  to  better  than  a  few  per  cent. 

Finally,  to  return  to  the  first  criterion,  we  look  briefly  at  the  relationship  between 
a  source  and  the  angular  field  of  view  of  an  instrument  (the  solid  angle  u  subtended  at 
the  entrance  pupil  by  the  entrance  window — see  figure  4.13).   Our  discussion  of  the 
experimental  situation  depicted  in  figure  4.6  (figure  2.2  of  [5]),  has  established  that  the 
flux  A$  through  the  two  beam-defining  apertures  is  the  same,  regardless  of  the  position 
and  orientation  of  the  source,  as  long  as  all  rays  of  the  beam  originate  on  an  iso-radiance 
surface.   Since  the  entrance  window,  of  area  A^,,   and  entrance  pupil,  of  area  A  ,   of  an 
optical  system  (see  figure  4.13)  are  equivalent  to  just  such  a  pair- of  beam-defining 
apertures  for  the  rays  that  will  pass  through  that  optical  system,  an  iso-radiance  source 
that  fills  the  entrance  window  or  field  of  view   (as  well  as  the  entrance  pupil)  with  rays 
of  the  same  radiance     L    will  produce  the  same  flux  through  the  optics  regardless  of  its 
position  and  orientation.      In  fact,  we've  already  used  the  exact  expression  for  the  flux  in 
such  a  beam  in  eq.  (4.40): 

$  =  L'/  /  D   «cose  -cosG  -dA  «dA^  =  L«0  [W] .  (4.40a) 

h   Ar  S 

The  throughput  0  between  entrance  window  and  entrance  pupil  is  a  constant  of  the  instru- 
ment.  It  is  the  same  as  the  throughput  between  the  aperture  stop  (same  as  entrance  pupil 
in  this  oversimplified  case)  and  the  field  stop.   Accordingly,  if  the  flux  responsivity 
R  (x,y,6,<f>)  of  the  instrument  is  not  a  variable  function  of  position  and  direction   but 
has  the  same  constant  value  R   for  every  ray  of  the  incident  beam,2  the  instrument  output 
can  be  treated  as  a  measurement  of  the  incident  radiance  L  =  $/0,   as  well  as  of  the  flux 
$  in  that  beam.   If  the  source  is  not  an  iso-radiance  source  so  that  the  radiance  cannot 
be  taken  outside  the  integral  in  eq.  (4.40a),  that  equation  can  be  written  with  the  variable 
radiance  inside  the  integral  as  $  =  jL«d0  and  the  quantity  $/6  then  gives  us  an 
average  value  of  the  incident  radiance,  thus 


1See  any  standard  text  on  geometrical  optics  [13,14].  A  future  chapter  is  planned  dealing 
with  stops,  pupils,  windows,  and  baffles  and  their  role  in  radiometry. 

2The  effects  of  a  variable  R  (x,y,6,<f>),  when  this  important  assumption  is  not  satisfied, 
will  be  covered  in  Chapter  5  on  the  measurement  equation. 
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APERTURE    STOP    8 
ENTRANCE    PUPIL 
(RECEIVING    AREA) 


ENTRANCE  WINDOW 
(IMAGE   OF  FIELD    STOP) 


Figure  4.13.    A  very  simple  radiometer  with  a  single  lens,  showing  the 
entrance  window  (image  of  field  stop  in  object  space)  and 
the  entrance  pupil  (image  of  aperture  stop  in  object  space- 
here  coincident  with  the  aperture  stop) . 

S  —  source 

S '  —  source  image 
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*/0  =   (l/0)»/L«d0  =  L  [W»m  2-sr  1] .  (4.46) 

The  foregoing  applies  only  when  radiation  from  an  extended  source  completely  fills  the 
entrance  window  or  field  of  view  (and  the  entrance  pupil)  and,  in  the  case  of  the  non-iso- 
radiance  source,  the  average  radiance  of  eq.  (4.46)  is  averaged  only  over  that  portion  of 
the  source  that  lies  within  the  entrance  window,  only  that  portion  that  contributes  rays  to 
the  beam  accepted  by  the  instrument  throughput.   On  the  other  hand,  when  a  source  fills 
only  a  portion  of  the  entrance  window,  as  illustrated  in  figure  4.13,  so  that  rays  from 
each  point  of  the  entire  exposed  portion  of  the  emitting  surface  reach  the  full  entrance 
pupil  and  are  included  in  the  measured  beam,  the  instrument  output,  as  a  measure  of  the 
flux  $  in  that  beam  {again  assuming  constant  flux  responsivity     R   that  is  the  same  for 
every  ray  of  the  incident  beam) , 1  when  divided  by  the  area  A   of  the  entrance  pupil 
(also  a  constant  of  the  instrument) ,  is  a  measure  of  the  (average)  irradiance 
E  =  $/A   [W*m   ]   at  the  entrance  pupil.   Then  if,  in  addition,  the  size,  distance,  and 
configuration  of  the  source  are  such  that,  as  discussed  earlier  with  respect  to  the  second 
and  third  considerations,  it  has  a  meaningful  radiant  intensity  I  at  that  distance  (from 
source  to  entrance  pupil) ,  producing  an  irradiance  E  that  is  inversely  proportional  to 
the  square  of  that  distance,  the  instrument  output  (again  assuming  constant     R*)1  may  also 
be  treated  as  a  measure  of  the  apparent  source  intensity  at  the  given  distance  D.   In 
general,  this  apparent  radiant  intensity  will  be  equal  to  the  actual  source  intensity  only 
if  there  is  no  attenuation  (loss  by  absorption  or  scattering)  along  the  optical  path  between 
source  and  measuring  instrument.   The  apparent  radiant  intensity  is  given  by 

I1  =     t*-I  =  D2«$/A   [W'sr-1],  (4.47) 

where 

I'   [W°sr  ]   is  the  apparent  radiant  intensity  of  the  source  as 
viewed  from  the  instrument  at  a  distance  of  D  [m] , 

I   [W*sr   ]   is  the  radiant  intensity  of  the  source  as  given  by 
eq.  (4.33),  and 

t*  is  the  propagance,  as  defined  by  eq.  (2.38)  [5],  of  the  optical 

path,  of  length  D  [m] ,   between  source  and  instrument  (entrance 
pupil  or  receiving  aperture)  where  the  measured  flux  is  $  [W] 

(R.  =  const,    is  the  same  for  all  rays^)    through  an  aperture  of 

a   r  2n 
area  A  [ml. 
r 

For  example,  the  source  quantity  (radiance  or  radiant  intensity)  measured  by  the  flux 
reaching  an  ideal  photocell,  with  uniform  and  isotropic  responsivity     R,  for  all  incident 


1The  effects  of  a  variable  R  (x,y,0,<}>)  ,  when  this  important  assumption  is  not  satisfied, 
will  be  covered  in  Chapter  5  on  the  measurement  equation. 
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rays j 1  positioned  to  receive  the  entire  beam  through  a  telescope  that  is  focused  on  the 
moon,  will  depend  on  the  magnification  used.   At  low  magnification,  where  the  entire  moon 
is  seen  through  the  telescope,  the  output  is  a  measure  of  the  irradiance  at  the  telescope 
objective  and,  in  combination  with  the  distance  to  the  moon,  a  measure  of  the  moon's 
apparent  radiant  intensity.   On  the  other  hand,  at  higher  magnification,  where  only  a 
portion  of  the  moon's  surface  is  seen  in  the  field,  the  output  is  a  measure  of  the  average 
apparent  radiance  of  just  that  portion  of  the  moon's  surface.   (Again  we  use  the  word 
"apparent"  to  include  the  effects  of  attenuation;  the  apparent  radiance  is  L'  E  t*«L.) 

One  further  consideration  that  may  occasionally  arise  is  the  effect  of  the  longitu- 
dinal extent  of  a  source  toward,  or  away  from,  the  receiver  (we  have  so  far  considered  only 
the  transverse  dimensions).   It  is  clear  from  the  discussion  accompanying  figure  2.10  [5] 
that  the  irradiance  produced  at  a  distant  point  by  an  iso-radiance  source  depends  only  on 
the  solid  angle  subtended  at  that  point  by  its  periphery  or  its  extremities  (as  "seen"  from 
the  receiver)  and  is  otherwise  independent  of  its  size,  position,  and  configuration  within 
that  solid  angle  (strictly,  it  is  proportional  to  the  projected  solid  angle  SI   =  Jcos6*doj). 
Hence,  longitudinal  effects  do  not  occur  for  iso-radiance  sources.   For  non-uniform  sources, 
we  can  obtain  an  indication  of  the  effects  of  longitudinal  extent  on  approximations  to  the 
inverse-square  law  by  examining  the  configuration  in  figure  4.14.   Two  "point-source" 
elements  of  an  extended  volume  source  (see  next  section)  each  of  radiant  intensity 
dl  [W*sr  ],   are  at  distances  of   (1+6)  «D  and   (1-6) *D  from  an  observation  point  0, 
where  the  resulting  element  of  normal  irradiance  or  radiant  fluence  rate  is 


dE   =  dF 
n       t 


dl 


(D  +  6-D)' 


dl 


(D  -  6-D)' 


2dl   riT  -2, 
7T-   [W«m  ]. 


(4.48) 


Figure  4.14.    Source-point  configuration  for  assessing  the 

effects  of  longitudinal  extent  on  approximations 
to  the  inverse-square  law. 


^he  effects  of  a  variable  R  (x,y,0  ,((>)  ,  when  this  important  assumption  is  not  satisfied, 
will  be  covered  in  Chapter  5  on  the  measurement  equation. 
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The  approximation  involved  is,  clearly, 

=-  + x-    «  2.  (4.48a) 

(1  +  &y  (1  -  6)Z 

It  is  easily  shown  that  the  approximation  will  be  within  1%  if  6  s  0.06.   The  approxi- 
mation improves  if  more  source  points  are  distributed  over  the  intervening  space. 

VOLUME-SOURCE  DISTRIBUTION:   RADIANT  STERISENT.   Up  to  this  point,  we  have  considered  only 
sources  that  have  emitting  surfaces  or,  at  least,  only  those  where  it  is  feasible  to  refer 
the  exitent  radiance  to  some  convenient  reference  surface.   However,  there  are  also  dis- 
tributed sources  where  radiation  is  emitted,  or  incident  radiation  is  scattered,  from 
points  throughout  a  volume  where  it  may  not  be  feasible  to  establish  the  source-radiance 
distribution  directly  with  respect  to  any  reference  surface.   For  example,  in  studying  the 
night-sky  airglow,1  the  radiance,  or  spectral  radiance,  due  to  the  emitting  layer  is 
observed  from  the  ground  where  the  radiation  coming  from  altitudes  of  about  50  [km]   or 
more  seems  to  come  from  an  emitting  surface  having  the  observed  radiance.   However,  the 
radiation  is  actually  coming  from  excited  atmospheric  molecules  distributed  throughout  an 
emitting  region  that  may  extend  vertically  over  an  altitude  range  measured  in  kilometers. 
Theoreticians  concerned  with  the  details  of  the  emission  processes  need  to  relate  the 
observed  sky  radiance  to  the  volume  distribution  of  source  emission,  so  we  need  a  radio- 
metric quantity  to  characterize  that  distribution. 

A  more  familiar  example  is  the  daytime  sky  where  the  radiance  (actually,  the  corre- 
ponding  luminance — see  table  4-2)  that  we  see  consists  of  scattered  sunlight  reaching  us 
from  points  distributed  throughout  a  huge  volume  of  air.  When  the  sky  is  clear  and  blue, 
the  scattering  points  are  primarily  the  air  molecules  themselves,  as  well  as  some  dust  and 
aerosol  particles;  when  the  sky  is  hazy,  misty,  or  "smoggy",  there  are  many  more  and  larger 
scattering  particles  of  various  aerosols  and  atmospheric  contaminants.  We're  not  con- 
cerned, however,  with  the  detailed  mechanisms,  only  with  the  fact  that  the  source  of  radia- 
tion is  widely  distributed  throughout  an  extended  volume. 

The  International  Lighting  Vocabulary  (CIE-IEC)  [6]  lists  no  radiometric  quantity  that 
is  a  volume-source  distribution.   We  have  adopted  the  term  radiant  sterisent  and  the  symbol 
L*,   proposed  by  Jones  [16,17]  and  by  a  Nomenclature  Committee  of  the  Optical  Society  of 
America  [18],  for  the  "generated"  (emitted  and/or  scattered)  radiance  per  unit  path  length 
along  a  ray,  or  the  equivalent  "generated"  radiant  intensity  per  unit  volume,  by  which  such 
a  distributed  source  can  be  characterized.   This  concept  will  now  be  developed. 


^Nocturnal  airglow  is  the  emission  of  light  from  extremely  high-altitude  atmospheric  regions 
where  gases  at  very  low  pressures  and  densities  have  been  excited  by  incident  solar  radia- 
tion during  the  daytime  so  that  they  luminesce  faintly  over  the  whole  sky  at  night  at 
certain  characteristic  wavelengths. 
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dL, 


Figure  4.15.   A  ray  path  through  an  extended  volume  of  radiating 
(emitting  and/or  scattering)  medium. 


Consider  a  region  in  which  radiation  is  being  generated  (emitted  and/or  scattered)  at 
all  points  throughout  an  extended  volume.   Let  the  point  0  in  figure  4.15  be  a  point  of 
observation  or  measurement,  and  consider  the  incident  radiance  at  0  from  a  single  direc- 
tion (along  a  single  ray) .   Let  the  distance  along  that  ray  from  0  to  any  point  P  on 
the  ray  be  s  [m] .   If  contributions  to  the  incident  radiance  at  0  are  generated  (emitted 

and/or  scattered)  along  the  ray  path  at  the  rate  of  L*(s)   units  of  radiance  per  meter  of 

-2   -1  -1 
path  length   [W*m  «sr   *m  ],   the  radiance  of  an  element  of  path  length  ds  at  P  in 

-2   -1 
the  direction  of  the  ray  to  0  is  dL  =  L*(s)«ds  [W»m  *sr   ].   If  there  is  no  attenua- 
tion (and  if  the  intervening  medium  is  uniform  and  isotropic — iso-ref ractive- index) ,  this 

element  of  path  ds  will  then  contribute  an  amount  dL  =  dL   to  the  observed  or  measured 

g 
radiance  incident  from  that  direction  (along  the  given  ray  from  P)   at  0,   because  of  the 

geometrical  invariance  of  radiance — see  Chapter  2  [5] .   However,  if  there  is  attenuation 

over  the  intervening  path  of  s  [m]   from  P  to  0,   so  that  the  propagance1  is  t*(s), 

the  incident  radiance  at  point  0  from  just  the  element  ds  at  P  is 


dL  =   x*(s)-dL   =  T*(s)«L*(s)«ds  [W-m2«sr1]. 


(4.49) 


The  observed  or  measured  incident  radiance  at  0,   from  the  direction  of  the  ray  through 
P,   is  then  the  integral  of  all  such  contributions  from  the  entire  path,  or 


L  =  /  L*(s)«T*(s)-ds  [W«m  2-sr  1] , 


(4.50) 


Accordingly,  referring  to  any  point  P  as  the  point  x,y,z  and  designating  the  ray 
direction  at  P(x,y,z)   by  the  direction  coordinates  0,<t>,   for  greater  generality,  the 
volume  distribution  of  generated  radiance,  as  a  function  of  position  and  direction  through- 
out any  such  distributed  source,  is  characterized  by  the  radiant  sterisent  (in  a  given 
direction  at  a  point) 


L*(x,y,z,9,cf>) 


dL  (x,y,z,6,<j)) 
ds 


[W«m   «sr   ], 


(4.51) 


^ee  Chapter  2  [5]  for  these  relations  involving  propagance  t*. 
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where 


-3   -1 

L*(x,y,z,G,<j>)  [W'm  *sr  ]   is  the  radiant  sterisent  at  the  point 

x,y,z  in  the  direction  8,cf>; 

-2   -1 
dL  (x,y,z,0,<j>)  [W*m  »sr   ]   is  the  generated  (emitted  and/or 

scattered)  radiance  of  the  medium  along  the  path  element  ds, 

in  the  direction  6,<j>,   at  the  point  x,y,z;   and 

2     2     2  ^ 
ds  =  (dx  +  dy  +  dz  ) 2  [m]   is  an  element  of  path  length  along 


the  given  ray  in  the  direction 


at  the  point   x,y,z   so 


that  8  =  cos   (dz/ds)   and  <J>  =  tan  (dy/dx).1 

-3   -1 
From  the  unit-dimensions,   [W«m  «sr   J,   it  appears  that  radiant  sterisent  is  radiant 

-1  -3 

intensity  [W*sr   ]   per  unit  volume   [m  ],   and  this  is,  indeed  the  case,  as  is  easily 

demonstrated.   Consider  a  cylindrical  element  of  volume  dV,   along  the  path  element  ds 

of  figure  4.15,  with  parallel  end  faces,  each  of  area  dA  and  tilted  at  an  angle  0 

between  the  normal  to  dA  and  the  ray  path  along  the  element  ds,   as  illustrated  in 

figure  4.16.   The  volume  element  dV  =  ds*cos6  *dA  is  indefinitely  small   (As-K)  and  AA-K)) 


Figure  4.16.   An  elementary  volume  along  the  ray  path  at  P   in 
figure  4.15. 


so  that  the  generation  (emission  and/or  scattering)  of  radiation  is  uniform  throughout  dV 

and  the  radiance  dL   is  the  same  over  all  parts  of  the  end-face  element  dA.   Then,  from 

g 
eq.  (4.32),  the  element  of  radiant  intensity  of  the  end  face  dA  in  the  direction  9,<J> 

is 


dl   =  dL  'COS0  -dA  [W«sr   ], 
g       g     s 


(4.52) 


so  that 


dL   =  dl  /(cos0  «dA). 
g       g      s 


(4.53) 


^The  relationships  between  spherical  and  rectangular  coordinates  are  summarized  in  Appen- 
dix 2  [5]  for  those  who  may  wish  to  refresh  their  memories  on  that  topic. 
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Then,  from  eqs.  (4.51)  and  (4.53),  we  can  write 

dL  dl  dl        _ 

L*  ~=    -£-    E  ds-cose  »dA  E  -df"  »-•«].  <4-54> 


s 


or,  more  explicitly 


dl  (x,y,z,e,<(>) 


L*(x,y,z,e,(())   e  8 — LI    [W.m3-sr1],  (4.55) 


where 


dl  (x,y,z,9  ,<(>)  [W«sr   ]   is  the  generated  radiant  intensity  of  the 

g  -3 

volume  element   dV  =  ds'cosG  -dA  [m   ]   at  the  point  x,y,z 

in  the  direction  9,<(>,   where 

6,<j>  [rad]  ,   as  before,  designates  the  direction  of  the  ray  (and  of 
its  path-length  element   ds)   at  the  point  x,y,z. 

In  developing  the  concept  of  radiant  sterisent  and  its  use  in  the  radiative-transfer 
relation  of  eq.  (4.50),  we  have  integrated  along  the  entire  path,  from  the  point  0  to 
infinity,  without  considering  that  there  might  be  something  other  than  the  distributed- 
source  medium  along  that  path.   Now,  we  want  to  look  at  the  case  where  the  distributed- 
source  medium  is  found  only  along  a  finite  path  extending  between  the  observation  or 
measurement  point  0   and  an  opaque  source.   Let  the  opaque  surface  have  a  radiance  L..   in 
the  direction  of  point  0  at  a  distance  s.   from  0  along  the  ray  path.   The  propagance, 
again,  is  given  by  x*(s)   for  0s  s  <.   s..  ,   but  now,  because  the  path  is  blocked  by  the 
opaque  surface  at   3  =  3^   x*(s)  =  0   for  all   s  >  s..  .   Then,  if  L*(s)   is  still  the 
sterisent  along  the  ray  path  in  the  intervening  medium,  the  observed  radiance  from  the 
direction  of  the  given  ray  incident  at  point  0  is  now  given  by 

g 
L  =  L1-t*(s1)  +  /  1L*(s)-x*(s)-ds  [W'm"2-sr_1].  (4.56) 

o 

The  integral  vanishes  for  all  values  of  s  greater  than  s..   because  the  propagance  x*(s) 
goes  to  zero,  so  it  is  immaterial  whether  the  upper  limit  is  shown  as  s..   or  as  infinity. 

While  equations  like  eq.  (4.56)  and  the  approach  used  in  setting  it  up  are  most  likely 
to  be  used  in  any  applications  to  evaluate  the  radiometric  quantities  involved,  some 
readers  may  be  interested  in  the  possibility  for  making  the  sterisent  a  more  general  dis- 
tribution function  that  can  account  for  discrete  as  well  as  continuously  distributed 
sources.   This  is  conceptually  appealing  because  eq.  (4.50)  then  becomes  a  very  general 
equation  of  radiative  transfer.   It  is  accomplished  by  introducing  the  discontinuous  dis- 
tribution function  known  as  the  Dirac  delta-function  [19].   For  example,  in  the  situation 
covered  by  eq.  (4.56),  let  the  sterisent  be,  instead, 
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L*(s)   =  L..-6(s  -  sn)  +  L  *(s)  [Win  3-sr  X]  (4.57) 

1        1     m 

where 

-3   -1 
L*(s)  [W«m  «sr     is  the  sterisent  of  the  intervening  medium;  and 
m 

6(s  -  s.)   is  a  Dirac  delta-function  that  satisfies  the  following 
defining  relations  [19] : 

6(u)  =  0  for  u  i   0, 

J 6(u) *du  =  1,   and 

Jf(u)-6(u-a)-du  =  f(a), 

when  the  integration  is  carried  out,  in  each  case,  over  any 
range  of  the  variable  that  includes  the  zero  of  the  argument 
of  the  6-function. 

It  is  easily  verified  that,  following  the  above  rules  for  the  6-function,  substitution  of 
eq.  (4.57)  into  eq.  (4.50)  leads  to  the  same  result  as  eq.  (4.56).  While  it  is  appealing 
and  conceptually  helpful  to  thus  treat  eq.  (4.50)  as  a  very  general  equation  of  radiative 
transfer,  any  attempt  to  apply  it  to  an  actual  situation  can  rapidly  develop  considerable 
complexity. 

Since,  as  we  have  shown,  the  radiant  sterisent  is  the  radiant  intensity  per  unit 
volume,  the  radiant  intensity  of  an  extended  volume  source  in  which  there  is  no  significant 
attenuation  is  just  the  volume  integral  of  the  sterisent 

l(0,(j>)   =  /  L*(x,y,z,6,(|))-dV  [W-sr-1],  (4.58) 

V 

3 
where  the  element  of  volume  dV  =  dx«dydz  [m  ] .   In  the  presence  of  appreciable  attenua- 
tion, on  the  other  hand,  it  is  necessary  to  first  establish  ray-radiances  at  an  appropriate 
reference  surface  by  eq.  (4.50)  and  then  to  compute  the  source  intensity  for  the  radiance 
distribution  at  that  reference  surface  by  eq.  (4.33).   Whether  this  is  useful  or  not  in 
either  case  depends,  of  course,  on  the  size-distance  relationships  already  discussed  with 
respect  to  the  intensity  of  extended  sources  and  the  approximations  to  the  inverse-square 
law. 

ENERGY  (TIME- INTEGRATED  FLUX)  DISTRIBUTIONS.   As  stated  in  Chapter  1  [5],  we're  mainly 
interested,  in  this  Manual,  in  the  measurement  of  radiant  flux  or  power  in  watts   [W] . 
However,  there  are  also  many  occasions  when  the  effects  of  interest  produced  by  a  radiation 
beam  are  related,  instead,  to  the  time- integrated  flux,  i.e.,  the  energy,  in  joules   [J]. 
A  very  common  example  is  that  of  photographic  exposure,  including  that  from  the  pulsed 
output  of  various  photoflash  illumination  devices.   The  incident  flux  reaching  a  particular 
area  of  a  photographic  film  may  vary  greatly  with  time  due  to  the  action  of  a  camera 
shutter  and/or  the  pulsed-light  output  of  a  flash  lamp  but,  within  rather  wide  limits  of 
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so-called  "reciprocity"  (between  time  and  irradiance) ,  it  is  the  time  integral  of  the 
incident  flux,  i.e.,  the  incident  energy,  that  determines  the  density  of  the  resulting 
photographic  image — along  with  other  factors,  such  as  temperature  and  "development"  (chem- 
ical processing).   The  radiation  "dose"  in  various  photobiological  or  medical  applications, 
such  as  photosynthesis  in  plants  exposed  to  sunlight  or  the  erythemal  or  "sunburn"  effect 
of  ultraviolet  radiation  on  the  human  skin,  also  depends  primarily  on  the  total  energy 
involved.   Accordingly,  we  take  a  brief  look  at  radiometric  quantities  that  are  spatial 

distributions  of  radiant  energy,  rather  than  radiant  flux.   Such  quantities  include 

_2 
radiant  (directed-surface)  exposure  H  [J*m  ] ,   radiant  (omni-directional)  fluence 

_2  -3 

F  [J«m   ],   and  radiant  (volume)  density  w  [J«m   ]. 

DIRECTED-SURFACE  DISTRIBUTION  of  INCIDENT  RADIANT  ENERGY;   RADIANT  EXPOSURE. 

Consider  a  time-varying  beam  of  radiation,  incident  on  a  given  surface.   If  the  total 

2 
energy,  incident  on  a  small  area  of  that  surface  AA  [m  ]   about  the  point  x,y  in  a  time 

interval  from  t,   to  t„  is  AQ[J],  we  define  the  radiant  exposure  (at  that  point  and  in 

that  time  interval)  as 

AA-+0   AA       dA 

Since  the  beam  is  varying  with  time,  the  radiant  flux  element  d$(t)   incident  on  each  area 
element  dA  is  a  function  of  time,  making  the  irradiance  at  any  point  x,y  also  a  func- 
tion of  time 

d§(t)      -2 
E(t)   =   Q^  ;   [W-m  Z].  (4.60) 

In  terms  of  the  time-varying  incident  flux  d$(t),   the  total  incident  energy  in  the  inter- 
val from  t..   to  t~   is  then 

fc2 
dQ     =     /     d$(t)«dt    [J],  (4.61) 

h 

so  that 

t2  t2 

H  =-  "dA-  ■  /  ^dT^--dt  =  /  E(tWt  [J]-  (4-62) 

tl  tl 

The  radiant  exposure  is,  thus,  the  time  integral  of  the  irradiance  [the  incident  radiant 
flux  (surface)  density].   In  fact,  each  of  the  relations  previously  given  between  irradi- 
ance [incident  radiant  flux  (surface)  density]  and  another  radiometric  quantity  transforms 
directly  to  that  for  radiant  exposure  by  expressing  the  other  quantity  as  a  function  of 
time  and  integrating  it  over  the  appropriate  time  interval.   For  example,  eq.  (4.7)  becomes 
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fc2 
dH(x,y,6,<j>)   =  dfi  •  /  L(x,y,0,<j>,t)'dt  [J«m   ].  (4.63) 

h 

Note,  however,  that,  in  unusual  situations  where  there  is  interaction  between  the  time 
dependence  and  the  effects  of  other  radiation  parameters,  as  in  some  scanning  instruments, 
such  transformations  cannot  be  made  so  simply. 

OMNI-DIRECTIONAL- SURFACE  DISTRIBUTION  of  INCIDENT  RADIANT  ENERGY:   RADIANT  FLUENCE. 

Radiant  fluence,  like  radiant  fluence  rate,  is  defined  in  terms  of  the  radiation  incident 

2   2 
on  a  small  sphere  of  radius  Ar  and  cross-sectional  area  Aa  =  TT-(Ar)   [m  ]   [eq.  (4.11)]. 

The  total  energy  incident  on  the  sphere  from  all  directions  in  a  given  time  period  of 

interest,  from  t..   to  t„,   is  AQ  [J].   Then,  if  we  keep  the  center  of  the  sphere  fixed 

at  the  point  x,y,z   (figure  4.3)  while  we  reduce  the  cross-sectional  area  Aa  by  reducing 

the  radius  Ar,   the  radiant  fluence  at  the  point  x,y,z  is  defined  as  the  limit,  as  Aa 

approaches  zero,  of  the  quotient  AQ/Aa   (since  the  CIE-IEC  [6]  have  no  symbol,  we  choose 

F,   used  by  photobiologists  [12]): 

'«*».■>  =  I£o-!2-  a  -£-i*--2i-  <*•«*> 


If  the  element  of  flux  incident  from  all  directions  at  the  point  x,y,z  is  varying  with 
time  d$(t" 
eq.  (4.61) 


time  d$(t),   the  total  incident  energy  in  the  interval  from  t..   to  t„   is  then  given  by 


h 


dQ     =     /     d$(t)-dt    [J].  (4.61) 

Then,    from  eqs.    (4.64),    (4.61),    and   (4.12), 


fc2  t2 

F     =~     "la"     *     /     -J7T5--    dt     =     /     Ft(t)'dt    [J'm_2]-  (4'65> 

Thus,  as  might  be  expected,  radiant  fluence  is  the  time  integral  of  the  radiant  fluence 
rate  at  the  same  point.   And,  again,  each  of  the  relations  previously  given  between 
radiant  fluence  rate  and  another  radiometric  quantity  transforms  directly  to  that  for 
radiant  fluence  by  expressing  the  other  quantity  as  a  function  of  time  and  integrating  it 
over  the  appropriate  time  interval  (except  when  there  is  interaction  between  the  time 
dependence  and  the  effects  of  other  radiation  parameters).   In  particular,  from  eq.  (4.20), 
we  can  write  the  relationship  between  elements  of  radiant  exposure  and  radiant  fluence  at 
a  point  x,y  on  some  reference  surface  as 

dF(x,y,0,<j>)   =  dH(x,y,e,<f>)/cose   =  dH  (x,y,0,<f>)  [J«m~2],  (4.66) 

where  dH   is  an  element  of  normal  exposure  and  is  equal  to  the  exposure  element  dH 
only  when  cos0   =  1,   i.e.,  when  0   =0   (when  the  ray  direction  6,cj>   coincides  with  the 
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normal  to  the  surface  element  dA  ) .   Thus  the  element  of  radiant  f luence  at  a  point  is 
equal  to  the  normal  exposure  at  that  point,  defined  as  the  exposure  on  a  surface  element 
that  is  normal  to  the  ray  at  that  point.   Only  elements  of  f luence  and  exposure,  associated 
with  individual  rays  (elements  of  throughput)  can  be  directly  related  in  this  way.   The 
integrated  total  fluence  and  exposure  at  a  point  cannot  be  directly  converted  from  one  to 
the  other.   If  values  for  both  are  needed,  they  are  most  easily  obtained  from  the  incident 
radiance  distribution  and  the  counterparts  of  eqs.  (4.16a)  and  (4.17).   But  this  seems 
unlikely;  ordinarily  only  one  or  the  other  will  be  useful  or  significant. 

Note  that  eq.  (4.65)  is  consistent  with 

Ft(x,y,z)   =   dF<y>z>   [Wm~2],  (4.67) 

justifying  our  use  of  the  notation  F   (which  means  =  dF/dt)1   for  fluence  rate. 

VOLUME  DISTRIBUTION  of  RADIANT  ENERGY:   RADIANT  (VOLUME)  DENSITY.  The  volume  distribution 

of  radiant  energy  is  not  included  in  the  CIE-IEC  nomenclature  [6].  However,  there  are 

some  occasions  when  it  is  used  in  radiometry,  and  it  is  included  in  the  ANSI  nomenclature 
[11],  so  we  include  it  here  for  completeness. 

The  radiant  (volume)  density  at  a  point  x,y,z  is  defined  as  the  quotient  of  the 
radiant  energy  contained  in  a  small  volume  about  the  given  point  by  the  magnitude  of  the 
volume,  as  that  volume  is  made  indefinitely  small: 

So  stated,  the  concept  is  clear  and  appears  quite  simple.  However,  when  we  try  to  relate 
it  to  other  radiometric  quantities,  we  must  take  into  account  the  fact  that  radiant  energy 
is  in  constant  motion,  propagating  at  high  speeds  that,  in  a  vacuum,  reach  the  well-known 
value  c  «  3  x  108  [m«s_1]  {=  (2.997  924  58  ±  0.000.000.012)  x  108  [m-s"1]}.   The  problem 
is  to  evaluate  the  instantaneous  radiant  energy  density  or  the  average  density  over  a  very 
short  time  period. 

In  order  to  do  this,  consider  a  ray  incident  at  the  point  x,y,z  from  the  direction 
6,4>   (see  figure  4.17).   If  the  radiance  is  L(x,y,z,6  ,<j>)  ,   the  element  of  fluence  rate 
and  the  flux  element  incident  on  an  elementary  sphere  of  cross  section  da  centered  at 
x,y,z  are  given  by  eqs.  (4.15)  and  (4.14),  respectively: 

dFt(x,y,z,e,cf>)   =  L(x,y,z,6,<j>)-dw  [W-m~2]  (4.15) 

and 

d$(x,y,z,e,<f>)   =  L(x,y,z,e,<|))«dwda  [W] .  (4.14) 


*In  this  Manual  X  =  dX/dp  for  any  radiometric  quantity  X  and  the  following  radiation 
parameters  p:   position,  direction,  a  spectral  parameter  (wavelength,  wave  number,  fre- 
quency v,   etc.),  or  time  or  frequency  f  <<  v. 
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ELEMENTARY    RIGHT-CIRCULAR   CYLINDER 
OF    LENGTH    ds    AND    CROSS-SECTION    do  - 


ELEMENTARY    SPHERE    OF 
CROSS-SECTION    do 


e, 


* 


Figure  4.17.    A  sectional  view  showing  a  ray  through  the  point  x,y,z   and 
along  the  axis  of  a  right-circular  cylinder  volume  element 
tangent  to  an  inscribed  spherical  volume  element  (see,  also 
figure  4.3) . 


Figure  4.17  also  shows  a  right  circular  cylinder,  of  length  ds  along  the  ray  as  its 
axis,  centered  at  x,y,z,   and  with  cross-sectional  area  da,   so  that  it  is  tangent  to 
the  inscribed  sphere,  also  of  cross  section  da  and  centered  at  x,y,z.   The  volume  of 
this  cylindrical  element  is 


dV  =  da-ds  [m  ] , 


(4.69) 


and  we  denote  the  amount  of  radiant  energy  contained  in  that  volume  at  any  instant  as 

dQ  [J].   In  terms  of  the  flux  element  d$,   the  energy  passing  through  a  transverse  element 

da  at  any  point  on  the  ray  in  an  elementary  time  interval  dt  is  dQ  =  d$«dt  [J]. 

But  the  time  required  for  all  of  the  energy  in  the  cylindrical  element  at  a  given  instant 

to  pass  out  through  one  end  is  just  dt  =  ds/v  =  n«ds/c.   This  is  the  time  for  the  energy 

at  one  end,  flowing  at  the  velocity  v  =  c/n  [m*s   ],   to  flow  just  the  distance  ds  to 

the  opposite  end.   Here,   n  S  c/v  is  the  index  of  refraction  of  the  medium.   Accordingly, 

the  energy  in  the  element  at  any  instant  is  given  by 
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dQ  =  d$-n-ds/c  [J].  (4.70) 

Now,  we  can  combine  eqs.  (4.68),  (4.14),  (4.69),  and  (4.70)  to  obtain  the  element  of  radi- 
ant (volume)  density  due  to  this  elementary  beam  along  the  single  ray  as1 

At  ,,  a\   -   dQ       n-d$ 

dw(x,y,z,e,cj,)   =  -gj-  =  -T75r- 

=   (n/c)-L(x,y,z,e,<f>)-du  [J-m~  ].  (4.71) 

But,  from  eq.  (4.15),  this  is  equivalent  to1 

dw(x,y,z,6,c|>)   =   (n/c)-dFt(x,y,z,6,<t.)  [J-m~3].  (4.72) 

The  same  is  true  for  the  elements  associated  with  all  other  rays  through  x,y,z   (at  a 
given  time)   from  all  directions  6,<j>,   so,  integrating  both  sides  of  eq.  (4.72)  with 
respect  to  solid  angle  (direction)   du  =  sin0-d6-dc)>,  we  can  write1 


/  dw(x,y,z,e,(f))   =   (n/c)-J  dFt(x,y,z,6  ,<j>) 


w(x,y,z)   =   (n/c).Ft(x,y,z)  [J-m3];  (4.73) 

the  radiant  (volume)  density  at  a  point  (at  a  given  time)  is  equal  to  n/c  times  the 
fluence  rate  at  that  point  (at  the  given  time).   Also,  from  eq.  (4.16a),  we  can  relate 
this  to  the  incident  radiance  by 

w(x,y,z)   =   (n/c)-/   (L(x,y  ,z,9  ,<)>)  -dco  [J-m"3].  (4.74) 

4ir 

Note,  particularly,  that  all  quantities  in  this  equation  are  evaluated  at  the  point 
x,y,z,   including  the  radiance  L.   If  the  radiance  has  a  value  L..   at  some  other  point 
x1  ,y1  ,z..   along  the  ray,  where  the  ray  direction  is  0..  ,<J>,   and  the  refractive  index  is 
nx,   then,  by  eq.  (2.23)  [5], 

2  -2   -1 

L(x,y,z,6,(f>)   =   (n/n1)  'I^Cx^y^z^e^il^)  [W-m   -sr   ] 


and,  furthermore,  this  value  must  be  multiplied  or  divided  (depending  on  the  direction  of 
propagation)  by  the  prop; 
eq.  (2.40)  or  (2.39)  [5] 


propagation)  by  the  propagance  t*   of  the  ray  path  between  x,y,z   and  x..  ,y1  , z  ,   as  in 


*Note  that  in  a  dispersive  medium,  where  n  =  n(X) ,   these  equations  must  be  written  for 
the  spectral  quantities;  for  example, 

-3   -1 
dw  (x,y,z,0,<j>,A)   =   [n(A)/c] -L.  (x,y,z,0,(J),A) -du  [J-m   -nm   ],  (4.71a) 


and 


-3   -1 
wA(x,y,z,X)   =   [n(X)/c]«F  x(x,y,z,X)   [J-m   -nm   ].  (4.73a) 
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SPECTRAL  RADIOMETRIC  QUANTITIES.   Although  this  chapter  is  devoted  to  spatial  distribu- 
tions, it  should  be  noted  that  all  of  the  spatial  radiometric  quantities  discussed  here 
are  transformed  to  the  corresponding  spectral  quantities  in  just  the  same  way  that  radiance 
was  transformed  to  spectral  radiance  in  Chapter  3  [5]. 

In  general,  if  X  is  any  radiometric  quantity,  including  L,  L*,  F  ,  M,  E,  W,  I,  $, 
w,  F,  H,  or  Q,   and  if  y  is  a  spectral  parameter  or  variable,  including  wavelength  X, 
frequency  v,   or  wave  number  a   (see  footnote  at  beginning  of  Chapter  3  [5]),  then  the 
corresponding  spectral  radiometric  quantity  is 

X   =  dX/dy  [X-y-1].  (4.75) 

For  example:   spectral  radiant  energy,  in  terms  of  wavelength,  is  C\  =  dQ/dX  [J«nm  ]; 

-3-1    -1  -1 
spectral  radiant  sterisent,  in  terms  of  wave  number,  is  L*  =  dL*/da  [W*m  «sr   «(cm  )   ] 

-3   -1  ° 

or   [W«m  «sr   «cm]   (i.e.,  in  watts  per  cubic  meter,  steradian,  and  reciprocal  centimeter); 

-3    -1  t 
and  spectral  radiant  density  in  terms  of  frequency,  is  w  =  dw/dv  [J*m  'THz   ].   These 

examples,  like  almost  all  of  our  discussion,  are  in  terms  of  energy  units,   [J]   and   [W] . 
One  should  keep  in  mind  that  they  apply  equally  to  any  other  form  of  flux  that  is  propa- 
gated in  accordance  with  geometrical  (ray)  optics.   For  example,  spectral  photon-flux 

-1  -2   -1 
exitance,  in  terms  of  wavelength,  is  (see  table  4.3)  M  ^  =  dM  /dX  [q«s   *m  «nm  ]. 

Spectral  quantities  could  also  be  formed  in  the  same  way  from  the  luminous-flux  or  photo- 
metric quantities  of  table  4.2,  but  we're  not  aware  of  any  useful  application  for  such 
quantities. 

SUMMARY  of  CHAPTER  4.   All  of  the  radiometric  quantities  in  table  4-1  have  now  been 
introduced,  with  definitions  and  discussions  of  their  significance  and  their  interrelations 
with  other  radiometric  quantities,  particularly  radiant  flux,  radiant  energy,  radiance,  or 
spectral  radiance.   The  "simple"  spatial  distributions  of  radiant  flux  include, 

(1)  the  directed-surf ace  distributions: 

radiant  flux  (surface)  density  W,   including 
irradiance  E ,   and 
radiant  exitance  M; 

(2)  the  omni-directional-surface  distribution: 

radiant  fluence  rate  F  ;  and 

(3)  the  directional  distribution: 

radiant  intensity  I. 

Important  interrelationships  include  the  one  between  radiance  and  flux  (surface)  density 
in  iso-radiance  beams  and  the  inverse-square  law  that  relates  point-source  radiant  inten- 
sity to  the  resulting  irradiance  as  a  function  of  distance,  as  well  as  approximations  to 
the  inverse-square  law  involving  extended  sources  and  receivers.   Also  presented  is  the 


+One  teraherz  l[THz]  =  1012[Hz]. 
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volume-source  distribution  radiant  sterisent  L*  and  its  use  in  a  radiative-transfer 
equation.   In  addition,  there  are  energy  (time-integrated  flux)  distributions  corresponding 
to  all  of  the  flux  distributions.   Those  given  special  names  and  symbols,  and  listed  in 
table  4-1,  include 

(1)  the  directed-surface  distribution  of  incident  radiant  energy: 

radiant  exposure  H; 

(2)  the  omni-directional-surface  distribution  of  incident  radiant 
energy: 

radiant  fluence  F;   and 

(3)  the  volume  distribution  of  radiant  energy: 

radiant  (volume)  density  w. 

Finally,  there  is  a  spectral  radiometric  quantity  X   corresponding  to  each  radiometric 
quantity  X  and  spectral  parameter  u. 

With  respect  to  the  spatial  parameters,  a  radiometric  measurement  situation  should 
always   be  analyzed,  first,  in  terms  of  the  distribution  of  ray-radiance  with  respect  to 
both  position  (across  a  reference  surface)  and  direction  (through  each  point  of  that 
surface)  to  account  for  the  full  range  of  possibilities  (see  Chapter  2  [5]).   If  it  is 
then  found  that  the  assumption  of  uniformity,  or  the  use  of  an  average  value,  with  respect 
to  part  of  that  distribution  will  not  produce  unacceptable  error,  a  "simpler"  distribution 
of  radiant  flux  with  respect  to  position  (area)  or  direction  (solid  angle)  alone  can  be 
used. 

The  flux  per  unit  area,  as  a  function  of  position,  is  given  in  two  ways:   The  first, 
the  directed-surface  distribution,  is  called  the  radiant  flux  (surface)  density  (at  a 
point  of  a  surface) 

W(x,y,)  =      d$(x,y)/dA  [w-nf2] ,  (4.1) 

where  the  directional  distribution  of  the  flux  d$(x,y)  [W]   through  the  surface  element 

2 
dA  [m  ]   at  the  point  x,y  need  not  be  specified.   This  quantity  is  related  to  the  radi- 
ance distribution  by 

W(x,y)   -   /  L(x,y,e,<J>)«dfi  =  /  L(x,y,0  ,cf>)  •cose-dw 
w  w 

=  /  /  L(x,y,6,c|>)'Cose«sine'de'd<j)  [W-m   ],  (4.4) 

<j)  e 

where  the  intersecting  rays  at  the  point  x,y  fill  a  solid  angle  w  ^  2tt  [sr]   that  may 
include  the  entire  hemisphere  on  one  side  of  the  element   dA   (of  the  tangent  plane  con- 
taining that  element,  if  the  reference  surface  is  curved).   For  an  incident  beam, 
W(x,y)  =  E(x,y)   is  called  the  irradiance;  for  an  exitent  beam,   W(x,y)  =  M(x,y)   is  the 
radiant  exitance. 
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The  second  form  of  flux  per  unit  area,  the  omni-directional-surface  distribution,  is 
called  the  radiant  fluence  rate  (at  a  point  in  space) 

Ft(x,y,z)   e  d<Kx,y,z)/da  [W-m~2],  (4.12) 

where,  again,  the  directional  distribution  need  not  be  specified  for  the  flux  d$(x,y,z)  [W] 

2 
incident  on  a  spherical  volume  element,  of  cross  section  da  [m  ]   and  centered  at  the 

point  x,y,z.   This  quantity  is  related  to  the  incident  radiance  distribution  by 
Ft(x,y,z)   =   /  L(x,y,z,6,<(>)  -dw 

=   /  /  L(x,y,z,6,c|>)-sin0-d0'd<|>  [W«m   ] ,  (4.16) 

c(>  0 

where  the  intersecting  rays  at  the  point  x,y,z  fill  a  solid  angle  w  <,   4tt  [sr]   that  may 
include  the  entire  sphere  of  directions  surrounding  the  point  x,y,z. 

Equations  (4.4)  and  (4.16)  can  be  turned  around  to  obtain  the  following  expressions 
for  radiance  in  terms  of  radiant  flux  (surface)  density  and  of  radiant  fluence  rate, 
respectively: 

t  t          a   ^      dW(x,y,9,(j))       dW(x,y,9,(j>)   fT7  -2   -ln  ..  ... 

L(x,y,0,<j>)   =  t^       =    cos£:^'   [W-m   -sr   ]  (4.18) 

and 

dF  (w,y,z,e,<J))       „   _1 

L(x,y,z,e,(t»)   = ^ [W-m  -sr   ].  (4.19) 

The  third  of  the  "simple"  distributions,  the  purely  directional  distribution  of 
radiant  flux,  is  the  radiant  intensity  (in  a  direction) 

1(0, <j>)   e  d$(e,<f>)/do)  [W«sr-1]  (4.21a) 

=   J  L(x,y,0,<j>)«cos0  -dA   [W-sr_1].  (4.33) 

A 
s 

Although  there  is  a  value  of  radiant  intensity  for  any  source,  by  eq.  (4.33),  that  quantity 

is  not  particularly  useful  except  to  characterize  a  source  in  a  configuration  where  it  may 

be  treated  as  a  "point  source";   i.e.,  at  distances  where  the  irradiance  that  it  produces 

varies  inversely  as  the  square  of  that  distance.   The  inverse  square  law  applies,  strictly, 

only  to  elementary  or  "point"  sources  and  receivers.   The  normal  irradiance  dE 

[e  dE/cos0  ]   on  a  receiver  element  dA   at  the  point  x  ,y  ,z  ,   or  the  radiant  fluence 

rate  dF   incident  on  an  elementary  receiver  volume  of  cross  section  da  at  the  point 

x  ,y  ,z  ,   at  a  distance  D  from  a  source  element  dA   at  the  point  x  ,y  ,z   of  radiant 
r  yr  r  s  s  s  s 

intensity  dl  in  the  direction  0,<j>  of  the  receiver  at  x  ,y  ,z  ,   at  an  angle  0 
from  the  normal  to  dA  ,   is  given  by 

dE   =  dF   =  dE/cos0   =  dl/D2  =  L«cos0  «dA  /D2  [W«m~2].  (4.36) 

n       t  r  s   s 
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More  explicitly, 

dEn(xr,yr,zr,6,<l>)   =  dFt(xr,yr,zr,e,<j>)   =  dE(xr>yr  ,zr  ,9  ,<\>)  /cose^ 


=  dI(xs,ys,Zs,e,<l))/D2 

2     -2 
=  L(x  ,y  ,z  ,6,<{))«cos0  -dA  /D   [W«m   ].  (4.36a) 

S    S    S  o     o 

When  relations  of  this  form  are  written  for  sources  and  receivers  of  finite  dimensions, 

the  distance  D  can  only  be  a  nominal  or  average  distance  because  the  slant  distance  D 

varies  for  different  pairs  of  source  and  receiver  elements,  so  the  relations  are  only 

approximations  to  the  inverse-square  law.   As  an  example,  a  simple  case  of  parallel,  plane, 

circular,  co-axial,  iso-radiance  source  and  ideal  receiver-detector,  of  area  A   and  A  , 

s        r 

respectively,  separated  by  a  distance  D  along  the  axis  between  their  centers,  is  examined. 
The  incident  flux  at  the  receiver  is 

$  «=a  E«A   <=»   (I/D2)«A   «  L»(A  «A  /D2)  [W] ,  (4.38) 

where  L  is  the  source  radiance,   I  is  its  radiant  intensity,  and  E  is  the  irradiance 
at  the  center  of  the  detector-receiver.   The  discrepancy  between  the  actual  and  inverse- 
square-approximation  values  of  flux  is 

A$   E   $  -  L«(A  «A  /D2)  [W].  (4.39) 

The  relative  error  of  the  approximation  to  the  inverse-square  law  is  A$/$  or,  in  terms  of 
throughput  (dividing  through  by  L) ,   A0/0.   For  the  conditions  stated,  it  is  shown  that 
this  error  is  less  than  one  per  cent  when  the  maximum  transverse  dimension  of  source  or 
receiver  is  no  more  than  one  tenth  of  the  distance  D  between  them.   In  fact,  the  use  of 
D,   the  distance  between  centers,  rather  than  the  slant  distance  D   between  a  pair  of 
single  points  on  source  and  receiver,  does  not  introduce  an  error  exceeding  one  per  cent  in 
the  inverse- square  relationship  within  these  size  limits,  even  for  extreme  points  at 
opposite  edges  of  source  and  receiver  (see  figure  4.11.)   A  brief  examination  shows  effects 
of  about  the  same  magnitude  due  to  the  longitudinal  dimensions  of  a  source.   The  irradiance 
from  a  uniform  iso-radiance  source,  however,  depends  only  on  the  projected  solid  angle  that 
it  subtends  (see  figure  2.10  [5]). 

Altogether,  three  main  considerations  determine  when  a  source  may  be  usefully  treated 
as  a  point  source,  characterized  by  its  radiant  intensity,  in  making  measurements: 

(1)  the  size  of  the  source  relative  to  filling  the  entrance  window  or 
angular  field  of  view  of  the  measuring  instrument  (radiometer) , 

(2)  the  size  of  the  maximum  transverse  dimensions  of  source  and  receiver 
(instrument  receiving  aperture  or  entrance  pupil)  relative  to  the 
distance  separating  them,  and 
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(3)   the  directional  distribution  of  radiance  from  each  source-surface 
element  and  the  effects  of  obscuration  of  portions  of  that  surface. 

These  three  considerations  are  discussed  in  some  detail. 

The  volume-source  distribution  is  defined  as  the  radiant  sterisent  (in  a  direction  at 
a  point) 

L*  =  dL  /ds  =     dl  /dV  [W«m"3«sr_1] .  (4.54) 

g         8 

More  explicitly,  the  radiant  sterisent  at  the  point  x,y,z  in  the  direction  0,<j>  is 

L*(x,y,z,6,<j>)  =     dL  (x,y,z,0,<f>)/ds  [W«m~3«sr_1] ,  (4.51) 

-2-1 
where  dL  (x,y,z,6,<j>)  [W*m   *sr   ]   is  the  generated  (emitted  and/or  scattered)  radiance  of 

the  medium  in  the  path  element  ds  at  the  point  x,y,z  and  in  the  direction  6,<j)   (along 

a  given  ray  in  that  direction).   Alternatively, 

L*(x,y,z,e,<j>)   =   dl  (x,y,z,e,<|>)/dV  [W-m^-sr"1]  ,  (4.55) 

where  dl  (x,y,z,0  ,<|>)  [W*sr  ]   is  the  generated  radiant  intensity  of  the  volume  element 

g        3 
dV  =  ds'cos0  *dA  [m  ]   at  the  point  x,y,z  in  the  direction  0,<|).   Thus  sterisent  is  the 

generated  radiance  per  unit  path  length  or  the  generated  radiant  intensity  per  unit  volume. 

Its  significance  is  most  clearly  apparent  in  the  radiative-transfer  equation  for  the 

observed  radiance  at  a  point  from  a  given  direction 

00 

L  =  J  L*(s)-T*(s)-ds  [W-m"2-sr_1],  (4.50) 


where  L*(s)   is  the  sterisent  in  the  direction  of  the  observation  point  at  a  distance  s 
from  that  point  along  the  ray  path  in  the  given  direction  and  t*(s)   is  the  propagance 
over  the  intervening  path.   This  is  a  general  relation  of  very  wide  application  if  the 
distribution  of  generated  radiance,  the  sterisent,  is  broadened  to  include  discontinuities, 
expressed  by  using  Dirac  delta-functions,  as  well  as  continuous  distribution  functions. 

The  radiant  intensity  of  an  extended  volume  source  with  negligible  internal  attenua- 
tion is  given  by 

1(9,4))   =  /  L*(x,y,z,e,<t>)«dV  [W-sr"1].  (4.58) 

V 

However,  in  the  presence  of  substantial  attenuation  within  the  source  medium,  it  is  neces- 
sary to  first  establish  ray- radiance  at  an  appropriate  reference  surface  by  eq.  (4.50)  and 
then  to  compute  the  source  intensity  from  the  radiance  distribution  at  that  reference 
surface  by  eq.  (4.33).   Whether  this  is  useful  or  not,  in  either  case,  depends,  of  course, 
on  the  size-distance  relationships  already  discussed  with  respect  to  the  intensity  of 
extended  sources  and  approximations  to  the  inverse-square  law. 
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There  are  energy  (time- integrated  flux)  distributions,  corresponding  to  all  of  the 
spatial  distributions  of  radiant  flux  but  only  the  more  commonly  used  have  been  given 
special  names  and  symbols.   In  many  cases  it  is  the  total  energy,  the  time  integral  of 
flux,  in  a  given  time  interval  that  is  of  most  interest  and  significance. 

The  directed-surface  distribution  of  incident  radiant  energy  is  the  radiant  exposure 
(at  a  point  of  a  surface) 

H(x,y)   E  dQ(x,y)/dA  =  J   E(x,y,t)-dt  [J-m2].  (4.62a) 

h 

It  is  the  time  integral  of  the  irradiance  at  the  given  point. 

The  omni-directional-surface  distribution  of  incident  radiant  energy  is  the  radiant 
f luence  (at  a  point) 

_2 
F(x,y,z)   =  dQ(x,y,z)/da  [J-m   ]  (4.64a) 

2 
=  /   Ft(x,y,z,t)-dt  [J-m   ].  (4.65a) 

tl 
It  is  the  time  integral  of  the  radiant  f luence  rate  at  the  given  point. 

The  volume  distribution  of  radiant  energy  is  the  radiant  (volume)  density  (at  a  point) 

w(x,y,z)   E   dQ(x,y,z)/dV  [J-m"3]  (4.68a) 

=   (n/c).Ft(x,y,z)  [J-m"3]  (4.73) 

=   (n/c)-J  L(x,y,z,6,4>).du)  [J-m"3].  (4.74) 

4ir 

It  is  equal  to  n/c  times  the  f luence  rate  at  the  given  point  (and  at  a  given  time) ,  where 
n  is  the  index  of  refraction  of  the  medium  and  c  is  the  vacuum  velocity  of  electromag- 
netic radiation.   (In  a  dispersive  medium,  similar  relations  exist  but  only  between  the 
spectral  quantities  at  each  wavelength.) 

Finally,  spectral  radiometric  quantities  are  formed,  corresponding  to  each  spatial 
radiometric  quantity  X,  just  by  taking  the  derivative  with  respect  to  the  appropriate 
spectral  (variable)  parameter  y : 

X   =     dX/dy.  (4.75) 

The  reader  is  also  reminded  that  all  of  the  foregoing  quantities  in  terms  of  energy   [J] 
and  energy  flux  or  power   [W] ,   including  the  spectral  quantities ,  have  their  counterparts 
in  terms  of  any  other  form  of  flux  that  is  propagated  in  rays  that  obey  the  laws  of  geo- 
metrical optics.   The  photometric  quantities,  in  terms  of  flux  in  lumens   [lm] , 
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corresponding  to  the  radiometric  quantities  of  table  4.1  are  listed  in  table  4.2.  The 
photon-flux  quantities,  in  terms  of  flux  in  quanta  per  second  [q«s  ],  are  similarly 
listed  in  table  4.3. 
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Chapter  5.    An  Introduction  to  the  Measurement  Equation 

by  Henry  J.  Kostkowski  and  Fred  E.  Nicodemus 

In  this  CHAPTER.    The  measurement  equation,  which  is  central  to  our  approach  to  all  of 
radiometry  (optical  radiation  measurements) ,  is  introduced  by  deriving  simple  measurement 
equations  for  three  illustrative  measurement  problems.   First,  however,  there  is  a  brief 
presentation  and  discussion  of  responsivity,  the  output  signal  of  a  radiometer  per  unit 
input  of  incident  radiation,  an  essential  factor  in  the  measurement  equation.   The  measure- 
ment equation  relates  the  instrument  output  resulting  from  a  measurement  to  the  distribu- 
tion of  incident  spectral  radiance  in  terms  of  all  of  the  radiation  parameters  as  well  as 
significant  environmental  and  instrumental  parameters.   By  appropriate  substitution,  the 
output  signal  is  similarly  related  to  the  radiometric  quantity  to  be  measured,  if  it  is 
other  than  the  incident  spectral-radiance  distribution.   The  examples  here  deal  only  with 
the  radiation  parameters  of  position,  direction,  and  spectrum  (wavelength) .   The  radiation 
parameters  of  polarization  and  of  time  or  frequency  of  scintillation  or  modulation  and  the 
instrumental  and  environmental  parameters  have  not  yet  been  treated  in  this  Manual.   The 
three  measurement  problems  are:   (1)  transferring  the  spectral-radiance  calibration  of  a 
tungsten-ribbon  lamp,   (2)  determining  the  spectral  irradiance  near  a  large  source,  and 
(3)  determining  irradiance  with  a  broad-band  (broad-spectral-band)  radiometer.   The  last 
problem  also  leads  into  a  brief  discussion  of  normalization.   A  general  discussion  sum- 
marizes and  enlarges  on  the  important  features  and  properties  of  the  measurement  equation 
brought  out  by  the  illustrative  problems.   A  set  of  orderly  steps  for  approaching  the 
solution  of  the  measurement  equation  is  presented.   The  limitations  of  the  measurement 
equation  developed  in  this  chapter  are  summarized. 

RESPONSIVITY.    The  response  of  an  instrument  used  for  optical  radiation  measurements,  for 
example  a  radiometer  or  a  photometer,  is  an  output  signal  S   that  may  have  many  different 
forms.   It  may  be  a  voltage,  a  current,  a  photographic-film  density,  or  just  the  deflection 
of  a  meter  needle  or  the  displacement  of  a  line  on  a  recorder  strip  chart.   Each  of  these 
"signals"  is  measured  in  different  units  so,  for  general  statements  that  apply  to  all  of 
them,  we  just  designate  them  as  units  of  output  signal1   [S]   (we  use  the  upper-case   "S" 
to  distinguish  them  from  the  unit  of  time,  the  second   [s]). 

In  order  to  interpret  or  evaluate  this  output  signal  as  a  measure  of  the  incident 
radiation  that  produces  it,  we  need  the  relationship  between  it  and  that  incident  radiation. 
That  relationship  is  stated  mathematically  by  the  function  called  the  instrument  responsivity, 


■^his  may,  variously,  be  called  the  instrument  "response"  or  "output".   The  nomenclature 
is  not  well  standardized.   However,  since  most  such  outputs  are  electrical  "signals",  we 
have  taken  the  term  and  symbol  from  the  International  Electrotechnical  Commission  (IEC) 
[20]  that  are  readily  compatible  with  treatments  of  electrical  "signal"  processing  and 
display.   For  those  who  may  be  interested,  this  choice  of  nomenclature  is  discussed  in 
more  detail  in  Appendix  4. 
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the  output  signal  per  unit  input  of  incident  radiation  [17].  Most  instruments  respond  to 
the  time  derivative  of  the  incident  radiant  energy,  i.e.,  to  the  radiant  flux1   $  [W] . 
Accordingly,  we  will  be  concerned,  mostly,  with  the  flux  responsivity 

R$  =     dS/d$  [S-W_1],  (5.1) 

the  units  of  output  signal  per  watt  of  incident  radiant  flux.   It  follows  that  the  element 
of  output  signal  produced  by  an  element  of  incident  radiant  flux  is  given  by 

dS  =  R$'d$  [S].  (5.2) 

On  the  other  hand,  when  "reducing"  measurement  data,  some  people  prefer,  instead,  to  use 

d$  =  dS/R$  =  K$-dS  [W],  (5.3) 

where  K  =  1/R^  [W«S   ]   is  sometimes  called  the  calibration  "constant"  although,  like 
the  responsivity,  as  we'll  see,  it  is  often  quite  variable. 

With  many  radiometers  the  output  signal,  and  hence,  the  responsivity,  varies  consider- 
ably with  the  position  and  direction  of  incoming  rays  at  the  receiving  aperture.   It  is 
easy  to  demonstrate  this  fact  by  moving  the  narrow  beam  of  a  small,  intense,  well-collimated 
source  in  front  of  the  instrument.  And  we've  already  seen  in  figure  3.1  (Chapter  3  [5]) 
that  responsivity  can  also  be  a  highly  variable  function  of  wavelength.   In  general,  then, 
responsivity  [eq.  (5.1)]  and  the  element  of  output  signal  [eq.  (5.2)]  must  be  written  as 
functions  of  both  the  spatial  and  spectral  radiation  parameters : 

R$(x,y,6,(j),X)   E  dS(x,y,6,((),A)/d$(x,y,0,((),A),  (5.1a) 

and 

dS(x,y,0,<(.,X)  =  R$(x,y,e,<|),A)-d$(x,y,e,<|>,A) 

=  R.(x,y,0,<)),X)«L1  (x,y,e,<j>,A)-cos6'dwdA'dA  [S] .  (5.2a) 

The  instrument  responsivity  is,  in  general,  determined  by   (1)  the  propagance  of 
optical  elements  and  paths  internal  to  the  instrument,   (2)  the  responsivity  of  the  detec- 
tor or  transducer  element  that  transforms  the  incident  radiant  flux  into  some  other  form 
of  "signal",  and   (3)  the  signal  processing  by  components  that  modify  (amplify  or  filter) 
the  signal  from  the  detector-transducer  to  produce  the  final  instrument  output  signal, 
including  the  final  signal  display  or  recording  components .   The  effects  of  all  of  these 


10thers  (e.g.,  photographic  films  and  instruments  for  measuring  light  pulses)  respond  to 
the  energy  Q  [J],  the  time  integral  of  the  flux.  For  them,  the  energy  responsivity  is 
similarly  defined  as  R  =  dS/dQ  [S«J  ].  (The 
listed  in  table  4-1;  see  also  Appendix  1  [5].) 


similarly  defined  as  R  =  dS/dQ  [S«J~  ].   (The  radiometric  quantities  and  units  are 
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factors  and  their  interactions  are  combined  in  the  overall  flux  responsivity  R  . 

Another  consideration  that  must  not  be  overlooked  in  real  situations  is  that  of  lin- 
earity. A  radiation  detector  is  linear  if  its  flux  responsivity  is  the  same  for  any 
magnitude  of  flux  incident  on  the  detector  as  long  as  the  flux  distribution  in  position, 
direction,  and  wavelength  (and  time  and  polarization)  is  not  changed.   Fortunately,  many 
radiation  detectors  are  quite  linear  over  wide  ranges  and,  often,  they  are  used  only 
within  these  linear  ranges.   On  the  other  hand,  some  instruments  are  deliberately  designed 
to  have  non-linear  outputs.   For  example,  a  logarithmic  output  is  sometimes  used  in  appli- 
cations calling  for  relatively  rapid  measurements  of  widely  differing  amounts  of  flux.   In 
this  chapter,  we  assume  linearity,  i.e.,  that  the  responsivity  (the  spectral-directional- 
positional  flux  responsivity1)  is  not  a  function  of  the  level  of  input  radiation  nor, 
hence,  of  the  output-signal  magnitude.   This  simplifies  the  discussion  of  other  aspects  of 
the  measurement  situation,  and  it  will  still  be  adequate  for  a  great  many  applications. 

We  again  remind  the  reader  that,  for  the  present,  we  are  ignoring  dependence  on  time 
and  polarization.   We  also  call  attention  to  the  fact  that  responsivity  may  sometimes  show 
hysteresis,  where  its  value  is  a  function  of  time  and  depends  on  the  history  of  the  instru- 
ment's exposure  to  radiation.   This,  too,  we'll  ignore  for  the  present. 

EXAMPLES  of  the  MEASUREMENT  EQUATION.   We  have  now  developed  enough  of  the  fundamentals 
so  that  we  can  introduce  the  measurement  equation.   This  important  equation  relates  the 
output  signal  from  a  radiometric  instrument  or  measuring  device  to  the  distribution  of 
radiation  incident  on  it.   The  first  step  in  solving  this  equation  is  to  modify  it  so  that 
it  contains  the  radiometric  quantity  to  be  measured.   The  equation  and  its  solution  depend 
greatly  on  all  of  the  relevant  physical  factors  that  we'll  refer  to,  collectively,  as  the 
measurement  configuration  and  on  the  quantity  to  be  measured.   The  measurement  configuration 
always  includes  one  or  more  radiation  sources,  the  measuring  instrumentation,  and  the 
propagation  path(s)  for  the  radiation  beam(s)  between  them.  Accordingly,  rather  than  to 
attempt  a  general  discussion  that  will  cover  all  of  the  complex  possibilities,  it  will  be 


*Note  that  these  parameter  modifiers  denote  only  a  functional  dependence  in  the  case  of 
responsivity,  not  a  derivative  (concentration  or  distribution)  as  with  a  radiometric 
quantity.   For  example,  spectral  radiance  is 

-2   -1   -1 
L   =  dL/dA  [W*m   «sr   *nm  ] 

A 

-2   -1 
with  different  unit-dimensions  than  those  of  radiance  L  [W-m  *sr   ];  but  spectral  flux 

responsivity  is  not  a  derivative  with  respect  to  wavelength,  it  is  just  the  flux  responsiv- 
ity as  a  function  of  wavelength,  with  no  change  in  unit-dimensions: 

R$(A)   =  dS/d$|x  s      (dS/dA)/(d$/dA)|x  =     SX/$JA  [S-W-1]. 

The  effect  of  the  modifier  "spectral"  is  only  to  emphasize  the  functional  dependence  on 
the  spectral  parameter  (in  this  case,  wavelength). 
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simpler  to  start  by  setting  up  and  discussing  the  measurement  equations  for  three  simple 
illustrative  problems:   (1)  transferring  the  spectral-radiance  calibration  of  a  tungsten- 
ribbon  lamp,   (2)  determining  the  spectral  irradiance  near  a  large  source,  and   (3)  deter- 
mining irradiance  with  a  broad-band  (broad-spectral-band)  radiometer. 

Problem  1.   Transferring  the  spectral- radiance  calibration  of  a  tungsten-ribbon  lamp. 

This  problem  commonly  arises  when  one  has  obtained  a  spectral-radiance-standard  lamp 
(usually  a  ribbon- filament  tungsten  lamp)  from  a  national  or  commercial  standards  laboratory 
and  wants  to  calibrate  another  similar  lamp  against  this  standard  for  use  as  a  working 
standard,  to  conserve  the  limited  working  life  of  the  reference-standard  lamp.   Such  a 
comparison  calibration  measurement  is  usually  performed  by  observing  the  two  ribbon- 
filament  lamps  in  succession  with  a  spectroradiometer. 

A  spectroradiometer  typically  used  for  such  comparison  measurements  consists  of  (see 
figure  5.1) 

(1)  focusing  optics  (a  mirror  system  is  usually  used,  the  lens   (0)   is  shown 
for  simplicity  of  illustration)  that  image  the  desired  portion   (W)   of  the  source  emitting 
surface  onto  the  entrance  slit   (S..)   of 

(2)  a  monochromator   (M) ,   that  selects  and  transmits  a  narrow  wavelength  band 
of  the  received  source  radiation  to 

(3)  a  detector   (P) ,   that  produces  an  output,  usually  an  electrical  signal 
current  or  voltage,  proportional  to  the  incident  spectral  radiant  flux,  which  is  then  fed 
to 

(4)  an  electronic  amplifier   (E) ,  the  output  of  which  is  displayed  or  recorded 
in  the  desired  final  form  by 

(5)  an  appropriate  readout  device   (R) . 

Since  we're  concerned  primarily  with  the  basic  concepts  here  in  Part  I,  we'll  postpone  any 
discussion  of  the  details  of  these  devices  and  components  until  Part  II — Instrumentation 
and  Part  III — Applications,  and  will  confine  our  attention  now  to  just  their  main  functions. 

Our  first  concern  is  to  determine  just  how  the  instrumentation  selects  or  determines 
the  incident  beam  to  which  it  responds  in  terms  of  the  spectral-ray  radiance  distribution 
at  its  receiving  aperture.   The  three  parameters  with  which  we  are  concerned,  for  the 
moment,  are  position,  direction,  and  spectrum  (wavelength): 

(a)  position  (at  the  receiving  aperture)  —  Only  rays  through  points  within  the 
receiving  aperture  are  accepted  and  measured;  all  others  are  excluded. 

(b)  direction  (at  the  receiving  aperture)  —  At  each  point  within  the  receiving 
aperture,  only  the  converging  rays  in  the  solid  angle  subtended  at  that  point  by  a  desired 
target  area  on  the  source  emitting  surface  (the  ribbon  filament)  are  accepted  and  measured. 
Accordingly,  the  throughput  of  the  system  is  that  of  the  beam  defined  by  the  receiving 
aperture  and  the  target  area  in  combination.   Without  examining  in  detail  the  further 
paths  of  these  rays  through  the  monochromator  until  they  reach  the  detector  beyond  the 
exit  slit,  we  recognize  that  the  responsivity  may  vary  from  ray  to  ray  by  writing  it  as 
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R(x,y,6  ,<j>)  ,   as  in  the  preceding  section. 

(c)  spectrum  (wavelength,  at  the  receiving  aperture)  —  The  acceptance  and 
response  of  the  instrument  to  different  wavelengths  present  in  the  incident  radiation 
depends  on  (1)  the  extent  to  which  those  wavelengths  appear  in  the  beam  emerging  from  the 
exit  slit  of  the  monochromator  and   (2)  the  spectral  responsivity  of  the  detector  on  which 
that  beam  is  incident.   The  spectral  propagance  of  any  external  optics  and  optical  paths 
usually  has  a  relatively  minor  effect  although  filters  that  substantially  limit  the  spec- 
tral pass  band  are  sometimes  used  along  with  a  monochromator.   However,  the  overall  band 
width  AA   in  which  the  spectral  responsivity  of  the  instrument  as  a  whole  is  significantly 
non-zero  is  usually  determined  primarily  by  the  spectral  characteristics  of  the  monochroma- 
tor together  with  the  spectral  responsivity  of  the  detector.   If  AA   is  the  spectral 
pass  band  of  the  monochromator  (and  any  associated  optics  and  optical  paths) ,   AA  s  AA  . 
They  are  equal  if  the  detector  responds  to  all  wavelengths  passed  by  the  monochromator, 
etc.;  if  not,  the  overall  spectral  band  width  AA  may  be  less  than  the  monochromator  band 

width  AA  . 
m 

Without  examining  its  details,  it  is  sufficient  for  this  discussion  to  assume  that 

the  monochromator  is  an  instrument  designed  to  select  and  pass,  to  the  detector  at  its 

exit  slit,  a  desired  portion  of  the  spectrum  of  radiation  incident  on  its  entrance  slit. 

That  portion  is  characterized  by  a  wavelength  setting  A   and  a  spectral  bandwidth  or 

pass  band  AA  =  A„  -  A  .   All  wavelengths  A  reaching  the  detector,  lie  in  the  interval 

As  A  s  A_;   all  wavelengths  A  <  A.,   and  A  >  A„  are  blocked  and  do  not  appear  in  the 

emerging  beam  at  the  exit  slit.   The  spectral  bandwidth  of  the  monochromator  A A   is 

primarily  controlled  by  adjusting  the  widths  of  the  slits,  and  may  also  be  a  function  of 

the  setting  A  . 

°       o 

This  rather  sketchy  description  of  monochromator  operation  has  assumed  some  ideal 
conditions  and,  for  our  simple  example,  we  continue  to  assume  that  they  exist.   However, 
the  reader  is  cautioned  that  strong  incident  radiation  outside  any  spectral  bandwidth  AA 
may  not  be  reduced  completely  to  zero  (completely  blocked)  by  an  instrument  or  component 
designed  for  that  purpose.1  A  monochromator  is  no  exception.   Real  monochromators ,  because 
of  internal  scattering,  always  have  some  scattered  radiation  of  unwanted  wavelengths 
present  in  the  exit  beam,  although  with  good  designs  this  can  be  kept  extremely  low. 
However,  when  attempting  to  measure  weak  radiation  of  given  wavelengths  in  the  presence  of 
strong  radiation  of  other  wavelengths,  the  scattered  radiation  can  become  a  significant 
factor.   This  is  particularly  true  when  the  wavelength  separation  is  small.   A  good  example 
is  the  difficult  problem  of  measuring  the  very  weak  terrestrial  solar  ultraviolet  irradia- 
tion at  about  290  [nm] ,  that  is  passed  by  the  nearly  opaque  layer  of  atmospheric  ozone, 
when  this  is  done  in  the  presence  of  significantly  stronger  radiation  at  slightly  longer 
wavelengths  to  which  the  ozone  (and  the  rest  of  the  atmosphere)  is  highly  transparent. 


1This  caveat  applies  generally  to  filtering  and  blocking  with  respect  to  any  parameter,  not 
just  the  spectral  parameter. 
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Here  the  level  changes  by  about  six  orders  of  magnitude  within  a  wavelength  interval  of  30 
[nm]. 

Now  we're  ready  to  set  up  the  measurement  equation  for  the  incident  distribution  of 
radiation  at  the  receiving  aperture  of  our  instrument.   Figure  5.2  shows  a  single  ray  from 
a  point  within  the  target  area  at  the  source  to  the  area  element  dA  at  the  point  x,y 
of  the  receiving  aperture.   It  is  incident  there  within  the  solid-angle  element  dw  from 
the  direction  6,<)>.   From  eq.  (3.10  [5]),  the  element  of  radiant  flux  d$  incident  on  the 
receiving  aperture  along  this  ray  of  incident  spectral  radiance  L.  (x,y,6  ,<(>,A)   is  given 

A 

by 


d$(x,y,6,<|>,A)   =  L  (x,y,e,<t>,A)«dA«cos6-dwdA  [W] . 


(3.10) 


dA(x,y) 


W  —  target  area  (on  lamp  filament) 

A  —  surface  (imaginary)  area  across  receiving  aperture 

L^(x,y,6,<(>,A)  —  incident  spectral  radiance  at  area 
element  dA(x,y)   of  A  within  solid-angle 
element  dco  from  direction  9,<f>  and  in  wave- 
length element  dX  at  wavelength  A . 

Figure  5.2.    Incident  ray  from  lamp  filament  to  receiving  aperture, 


64 


We  can  think  of  such  an  element  of  flux  arriving  at  the  receiving  aperture  along  each  of 
the  rays  that  connect  every  point  of  the  target  area  with  all  points  across  the  receiving 
aperture.   Each  element  of  incident  flux  produces  a  corresponding  element  of  output  signal 

dS(x,y,e,(f.>X,Xo)   =  R$-d$ 

=     R.(x,y,e,<j),A,A    )«L.  (x,y,6  ,<|>  ,A)  •cos0*dwdA«dA    [S] .  (5.4) 

$  O         A 

where  R   is  the  flux  responsivity  of  the  spectroradiometer,  when  set  for  a  wavelength 
A  ,   for  this  particular   flux  element  (for  the  ray  from  the  direction  6,(|))   at  the  point 
x,y  of  the  receiving  aperture  and  at  wavelength  A.   The  total  output  signal  S  of  the 
spectroradiometer  is  then  the  sum  (integral)  of  all  such  signal  elements  dS  for  all  rays 
contained  in  the  optical  beam  accepted  by  the  measuring  instrument.   Included  are  all 
points  across  the  receiving  aperture  and,  at  each  point,  the  directions  of  all  rays  from 
every  point  of  the  target  area.   In  addition,  the  integration  must  be  carried  out  over  all 
wavelengths  for  which  the  responsivity  is  not  zero.   Thus  the  total  signal  is  given  by 

S(A,co,AA,Ao)   =  /   /  /  R$'Lx«cose«dwdA«dA  [S]  ,  (5.5) 

AJ  An 

where 

w  [sr]   is  the  target  solid  angle,  the  solid  angle  enclosed  by  the  extreme  rays 
from  the  target  area  that  converge  at  a  point  x,y  of  the  receiving  aperture 
[accordingly,   w  can  be,  and  usually  is,  a  function  of  the  position  x,y  and 
so  could  be  written  more  explicitly  as  w(x,y)], 

AA  [nm]   is  the  wavelength  interval  over  which  the  instrument  spectral  flux 
responsivity  is  significantly  non-zero,  i.e.,  where  R*(^  »X)  ^  0  when  the 
instrument  is  set  on  wavelength  A  . 

This  equation  [eq.  (5.5)]  is  the  measurement  equation  for  this  problem.   It  relates  the 
output  signal  of  the  measuring  instrument  to  the  spectral-radiance  distribution  incident 
at  the  receiving  aperture  of  the  instrument. 

This  measurement  equation  is  similar  to  eq.  (3.11)  [5],  which  gives  the  total  flux  in 
a  beam  at  its  intersection  with  a  reference  surface  as 

$  =  J  /  /  LA(x,y,6,(f),A)'dA'cose'dwdA  [W] .  (3.11) 

A  a)  AA 

The  difference  is  that,  in  the  measurement  equation,  each  element  of  flux  is  "weighted"  by 
the  corresponding  value  of  the  responsivity.   In  other  words,  we  can  regard  the  total 
output  signal  S  as  the  responsivity-weighted  flux  in  the  beam  received  by  the  instrument. 
For  the  measurement  equation,  the  reference  surface  is  the  instrument  receiving  aperture. 

The  first  step  in  solving  the  measurement  equation  for  a  particular  problem  is  to 
introduce  into  that  equation  the  radiometric  quantity  one  is  interested  in  measuring.   In 
this  problem,  the  quantity  to  be  measured  is  the  spectral  radiance,  but  it  is  the  spectral 
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£  1 

radiance  L   at  the  lamp1  rather  than  the  spectral  radiance  L   incident  at  the  receiving 
aperture  of  the  spectroradiometer  (see  figure  5.3).   If  we  assume  that  the  index  of  refrac- 
tion (that  of  air)  is  the  same  at  both  these  points,  eq.  (3.14)  [5]  gives  their  relationship 
as 


Lx(x,y,e,4>,A) 


L  (x,y,e,<J>,X)«T*(x,y,e,<|>,A)  [W-m~2«sr  1'nm  -1]  , 


(5.6) 


W  —  target  area  (on  ribbon  filament   F  of  lamp  Z) 
A  —  surface  (imaginary)  across  receiving  aperture 

£ 

L.  —  spectral  radiance  of  ray  just  outside  of  lamp  envelope 

L  —  incident  spectral  radiance  of  same  ray  at  receiving 

A 

aperture  after  traversing  path  of  propagance  t*. 


Figure  5.3.    Spectral  radiances  at  points  along  a  single  ray  from  target 
area  to  receiving  aperture. 


Because  of  the  way  in  which  radiance  standards  are  used,  the  spectral  radiance  is  specified 
at  the  outer  surface  of  the  lamp  envelope.   This  value  of  radiance  is  a  property  of  the 
lamp  (and  the  operating  conditions  —  voltage  or  current,  etc.)  regardless  of  the  propagance 
between  it  and  a  remote  instrument. 
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where 

T*(x,y,6,<j>,A)   is  the  spectral  propagance  over  the  intervening  path  along  the 

specified  ray  (between  the  lamp  envelope  and  the  receiving  aperture  of  the 

instrument) • 

I 
Thus,  the  equation,  as  finally  modified  to  contain  the  radiometric  quantity  L^  ,   that 

we're  actually  interested  in  measuring  in  this  problem,  is 

S(A,io,AA,A  )   =  /  /  /  R$*T*'LJJ'cose-dwdA-dA  [S].  (5.7) 

AX  A  a) 

Z  s 

Now  we  want  to  compare  L.   to  the  spectral  radiance  L.   of  a  standard  lamp.   The 

corresponding  equation  for  a  similar  observation  with  a  standard  lamp  is 

SS(A,uS,AX,Ao)   =  J  /  /  R$«T*S«L^'COse'dwdA«dA  [S],  (5.8) 

AA  A  w 

where,  by  the  use  of  the  superscript  s,  we  have  allowed  for  the  possibility  that  the 

g 

target  area  and  its  subtended  solid  angle,  the  target  solid  angle  u  ,   as  well  as  the 

g 
propagance  t*  ,   associated  with  the  measurement  of  the  standard  lamp  may  be  different. 

On  the  other  hand,  we  don't  need  the  superscript  s  on  the  remaining  quantities  since  the 

receiving-aperture  area  A,   and  the  wavelength  setting  A   and  pass  band  AX  and, 

hence,  the  spectral  responsivity  R*(^  >a)   of  the  spectroradiometer  are  all  kept  the  same 

for  both  measurements  [eqs.  (5.7)  and  (5.8)]. 

I 
We  can  easily  solve  eqs.  (5.7)  and  (5.8)  for  the  spectral  radiance  L.   of  the  uncali- 

g 
brated  lamp,  in  terms  of  the  spectral  radiance  L^   of  the  standard  lamp  and  of  the  two 

g 
observed  output  signals  S  and  S  ,   if,  for  each  lamp,  the  value  of  spectral  radiance  is 

the  same  everywhere  in  the  beam  and  at  all  wavelengths  in  the  interval  AX .  When  that  is 

a  s 

true,  we  can  bring  L,   and  L^   outside  the  integrals  in  each  case,  resulting  in 

A         A 

S  =  IyJ  /  J  Rd.-T*-cos0-dwdA-dA  [S]  (5.9) 


AA  A  a) 


and 


S   =  Iv/  /  /  gR$'T*S-cose'da)«dA'dA  [S].  (5.10) 

AA  A  a) 

By  positioning  the  standard  lamp  filament  in  exactly  the  same  location  as  that  of  the 
uncalibrated  lamp  when  it  was  measured  and  leaving  the  focusing  adjustment  and  position  of 

g 

the  spectroradiometer  unchanged,  the  corresponding  target  (solid)  angles  w  and  u   and 

g 
propagances  t*  and  x*   will  also  be,  respectively,  the  same.   In  other  words,  by 

making  the  measurements  under  identical  conditions,  the  integrals  in  eqs.  (5.9)  and  (5.10) 

are  identical  and  dividing  one  of  these  equations  by  the  other  gives 

L*  =   (S/SS)-L^  [W-m~2-sr~1.nm""1].  (5.11) 
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This  simple  relationship  is  often  used  in  making  a  comparison  calibration  of  two  such 
lamps.   By  deriving  it,  in  this  way,  from  the  full  measurement  equation  (at  least  in  terms 
of  the  parameters  of  position,  direction,  and  spectrum),  we  can  see  that  it  has  a  solid 
basis  and,  in  addition,  we  are  explicitly  aware  of  all  of  the  assumptions  or  conditions 
that  must  be  satisfied  for  it  to  be  valid. 

The  requirement,  in  deriving  eq.  (5.11),  for  constancy  of  spectral  radiance  throughout 
the  beam  (A  and  co)   and  at  all  wavelengths  in  the  interval  AX  for  these  ribbon-filament 
lamps  is  never  actually  realized.   There  are  temperature  and  emissivity  gradients  on  the 
tungsten  surface  that  produce  variations  in  spectral  radiance  relative  to  both  position  and 
direction.   In  addition,  the  spectral  radiance  for  tungsten  varies  significantly  with 
wavelength.   The  variations  relative  to  position  and  direction  may  be  minimized  by  making 
A  and  w  sufficiently  small  so  that  the  range  of  variation  is  about  one  per  cent  or  less. 
Then,  in  the  measurement  equation,  eq.  (5.7), 

S(A,w,AX,Ao)   =   /  £*•(/  /  R$«T*«cose«dwdA)«dX  [S]  (5.12) 

AX     A  go 

=  /  lJ(X,A,m)-bJ(X  ,XtA,«)-dX  [S],  (5.13) 

AX 

where 


J  J  L  «R  •T**cos6*dwdA 


L*(X,A,oj)   =  — M [Win  2-sr  1-nm  1]  (5.14) 


A 


/  /  Rft«T*'COS0«dwdA 
A  oj 


and 


•A 


R7(X  ,X,A,w)   =  J "/  R.'T*'COse«d(jj«dA  [S-W   «m  «sr]  .  (5.15) 

A  u 

-X,  I 

L  (X,A,o))   is  the  weighted  average  of  the  lamp  spectral  radiance  L..   over  A  and  u), 

I 
weighted  by  R  •t**cos0,   and  the  lamp-spectral-radiance  responsivity  R.  (X  ,X,A,uj)   is  the 

integral  over  A  and  w  of  this  weighting  function.1  A  similar  set  of  equations  exists 

-s  -2   -1   -1 

for  the  average  spectral  radiance  of  the  standard  lamp  L^(X,A,w)  [W«m  «sr   «nm  ]. 

The  variation  with  respect  to  wavelength  can  not   be  treated  in  exactly  the  same  way 
because  the  flux  responsivity  R   and,  therefore,  the  radiance  responsivity  R^  both  vary 
all  the  way  from  zero  to  their  respective  peak  values  in  the  wavelength  interval  AX. 


^ote  that  this  is  the  radiance  responsivity  to  the  total  beam;   R^  =  dS/dL  where  dS  is 
an  elementary  change  in  the  output  signal  and  dL  is  an  elementary  change  in  the  value  of 
radiance  of  an  incident  iso-radiance  beam  that  produces  the  signal  element  dS. 
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However,  lamps  such  as  these  are  usually  operated  so  that  their  relative  spectral  distribu- 
tions are  virtually  the  same.   In  this  case, 

L*(A,A,u>)   =  K  •JtjJ(A,A,a>)  (5.16) 

and 

L^(A,A,u))  =     K2-I®U,A,uj)  «  K2-I£(X.A,u),  (5.17) 

-£       -s 
where  £^   and  JL   are  the  relative  average  spectral  radiances  (averaged  over  A  and  w 

and  adjusted  to  have  the  same  maximum  value)  of  the  two  lamps  and  K..   and  K„  are  con- 
stants.  The  next  step  is  to  substitute  eq.  (5.16)  into  eq.  (5.13)  and  eq.  (5.17)  into  the 
standard-lamp  equivalent  of  eq.  (5.13)  and  take  the  quotient  of  the  resulting  expressions 
to  obtain 


r      -SI       I 
S(A,uj,AA,A  )  1  ;    A  X 


o'  AA  "  Kl 


SS(A,u,AA,Ao)  K2«J   £*-R£«dA       K2 


(5.18) 


But,  from  eqs.  (5.16)  and  (5.17), 

K±/K2     =  lZx    /  LSX    .  (5.19) 

Accordingly,  combining  eqs.  (5.18)  and  (5.19),  we  have 

L*  =   (S/SS)-L^  [W-m"2-sr"1-nm"1].  (5.20) 

A  A 

Equation  (5.20)  is  the  same  as  eq.  (5.11),  except  that  here  the  two  spectral  radiances, 
unknown  and  standard,  are  averages  over  A  and  u,   as  defined  in  eq.  (5.14). 

The  primary  objective  of  the  treatment  just  presented  is  to  introduce  the  measurement- 
equation  approach.   Discussion  of  further  details  of  actual  applications  of  eq.  (5.11)  or 
(5.20)  will  be  postponed  for  later  chapters.   However,  we  would  like  to  add  that,  by  limiting 
the  target-area  dimensions  used  in  such  comparison  calibrations  to  one  millimeter  or  less , 
aperture  angles  (subtended  by  the  instrument  receiving  aperture  at  the  lamp  source)  to 
about  0.01  steradian,  and  spectral  pass  band  to  about  2.5  nanometers,  the  uncertainties  in 
such  comparisons  have  been  limited  at  NBS  to  much  less  than  1%.   It  should  also  be  noted 
that,  when  average- spectral-radiance  values  are  used  in  a  maximum-accuracy  calibration, 
they  apply  only  to  the  particular  beam  from  the  particular  target  area   (W)   observed. 
Hence,  means  must  be  provided  for  accurately  identifying  a  desired  portion  of  the  lamp 
filament  and  a  specific  solid  angle  to  insure  identical  beam  geometry. 
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Problem  2.   Determining  the  spectral  irradiance  near  a  large  source. 

This  second  measurement  problem  is  illustrated  in  figure  5.4.  We  want  to  measure  the 
spectral  irradiance  produced  by  a  horizontal  fluorescent  lamp  at  a  point  x  ,y   on  a 
horizontal  plane  surface  beneath  it,  as  shown.   The  vertical  distance  z   from  the  lamp  to 
the  horizontal  plane  of  x  ,y   is  approximately  equal  to  the  length  of  the  fluorescent 

tube,  so  the  solid  angle  w   subtended  by  the  lamp  at  that  point  is  much  greater  than  the 

g 
solid  angle  oi   subtended  there  by  a  standard  of  spectral  irradiance  (e.g.,  a  small  cali- 
brated tungsten-halogen  lamp)  at  the  distance  for  which  it  is  calibrated  (usually  0.5  [m] ; 
see  figure  5.5).   As  we'll  see  presently,  we  don't  need  the  exact  values  of  these  solid 
angles  so  these  approximate  statements  will  suffice  for  the  description  and  analysis  of  the 
measurement  configuration. 


(xo,  yG) 


Figure  5.4.    Configuration  for  measuring  spectral  irradiance 

at  x  ,y   from  fluorescent  lamp  X. 
owo  r 


(x0,y0) 


Figure  5.5.  Configuration  for  spectral-irradiance 
calibration  measurement  with  standard 
lamp  £  . 
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An  apparatus  commonly  used  at  NBS  to  make  such  measurements  of  spectral  irradiance  is 
shown  in  figure  5.6.   It  consists  of  a  spectroradiometer  and  a  small  averaging  sphere.1 
The  spectroradiometer  is  just  like  the  one  in  figure  5.1,  used  to  compare  or  measure  spec- 
tral radiance  in  problem  1.   The  averaging  sphere  is  a  hollow  sphere,  typically  a  few 
centimeters  in  diameter,  with  uniform  diffusely  reflecting  inner  walls  and  two  ports,  as 
shown.   Radiation  incident  through  the  entrance  port  strikes  a  portion  of  the  inner  wall 
from  which  no  rays  can  be  reflected  directly   into  the  beam  to  the  spectroradiometer.   The 
spectroradiometer  is  focused,  by  its  external  optics,  with  its  target  area   (W)   entirely 
within  the  exit  port  of  the  sphere,  as  shown  in  figure  5.6.   The  function  of  the  averaging 
sphere  is  to  produce  isotropic  spectral- flux  responsivity  over  a  very  wide  solid  angle  of 
incidence.   The  multiple  internal  diffuse  reflections,  each  scattering  the  incident  radiation 
to  all  parts  of  the  sphere  wall,  from  which  it  is  again  so  redistributed,  produce  uniform 
irradiation  of  the  wall  opposite  the  exit  port,  regardless  of  the  direction  of  the  entering 
rays  through  the  entrance  port,  within  very  wide  limits.   Accordingly,  the  reflected  rays 
from  that  wall  toward  the  spectroradiometer  through  the  exit  port  all  have  the  same  value 
of  spectral  radiance,  and  that  value  is  proportional  to  the  spectral  flux  in  the  beam 
incident  through  the  entrance  port.  We  have  described  the  averaging  sphere  as  having  ideal 
properties.   With  a  real  sphere,  corrections  may  be  required  for  the  highest  accuracy,  but 
we  won't  be  concerned  with  them  in  this  chapter. 

The  measurement  is  made  with  the  entrance  port  of  the  averaging  sphere  in  a  horizontal 
plane  and  with  the  center  of  the  entrance  port  at  x  ,y  .   That  entrance  port  is  the 
receiving  aperture  for  the  entire  instrument,  limiting  the  incident  beam  to  be  measured. 
Again,  we  are  concerned  only  with  the  three  radiation  parameters  —  position,  direction, 
and  spectrum  (wavelength)  —  for  the  incident  rays  at  that  receiving  aperture: 

(a)  position  (at  the  receiving  aperture  of  the  entire  instrument)  —  By  its 
uniform  diffusion  of  the  incident  radiation,  the  averaging  sphere  insures  the  same  full, 
uniform  irradiation  of  both  the  target  area   (W)   and  the  receiving  aperture   (A)  of  the 
spectroradiometer   (see  figure  5.6)  by  all  rays  entering  the  entrance  port  of  the  sphere, 
regardless  of  position  there.   This  makes  the  overall  flux  responsivity  of  the  system  R. 
independent  of  position  x,y  within  the  entrance  port  (receiving  aperture  of  entire  instru- 
ment).  Also,  the  limits  of  the  beam  are  the  boundaries  of  the  entrance  port. 

(b)  direction  (at  the  receiving  aperture  of  the  entire  instrument)  —  In  the  same 
way,  the  overall  flux  responsivity  R   is  independent  of  incident-ray  direction  0,<(>  at 
the  receiving  aperture  (entrance  port  of  sphere)  over  wide  limits,  at  least  to  the  solid 
angle  subtended  there  by  the  fluorescent  lamp.   The  rays  of  the  desired  beam  from  the 
source  are  limited  by  the  dimensions  of  the  source,  not  by  the  measuring  instrument  which 
can  usually  accept  rays  over  a  wider  solid  angle.   Accordingly,  it  is  necessary  to  make 
sure  that  there  are  no  other  sources  present  from  which  a  significant  amount  of  additional 
incident  radiation  might  produce  an  erroneous  measurement.  When  this  is  not  possible,  as 
at  infrared  wavelengths  where  all  matter  at  room  temperature  is  emitting  thermal  radiation 


Sometimes  called  an  integrating  sphere. 
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(A) 

spectroradiometer 

^- 

(monochromator 
etc.) 


X  —  fluorescent  lamp 

S  —  averaging  sphere  (not  drawn  to  scale;  sphere  is  usually  only  a 
P 

few  centimeters  in  diameter) 

E  —  receiving  aperture  (entrance  port  of  averaging  sphere  S  ) 
n  p 

E  —  exit  port  of  averaging  sphere  S   [contains  target  area   (W) 
x  p 

for  spectroradiometer  (see  figure  5.1)] 
(A)  —  receiving  aperture  of  spectroradiometer  of  figure  5.1 


Figure  5.6.   Sketch  of  some  essential  parts  of  the  configuration  for  a  measurement 
of  spectral  irradiance. 
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in  significant  amounts,  special  techniques  for  cancelling  the  effects  of  unwanted  "back- 
ground" radiation  are  used  to  obtain  the  desired  measurement  results.  We  won't  go  into  the 
details  of  those  techniques  now,  however,   In  any  event,  the  incident  rays  of  the  desired 
beam  are  those  within  the  solid  angle  subtended  at  each  point  of  the  receiving  aperture 
(entrance  port  of  sphere)  by  the  source. 

(c)  spectrum  (wavelength,  at  the  receiving  aperture  of  the  entire  instrument)  — 
Ideally,  the  aver aging- sphere-wall  coating  is  selected  to  have  as  nearly  uniform  spectral 
reflectance  (as  a  function  of  wavelength)  as  possible.   The  wavelength  limits,  consequently, 
are  those  established  by  the  spectral  bandwidth  AA  at  the  wavelength  setting  A   of  the 
spectroradiometer.   Any  small  spectral  non-uniformity  of  reflectance  in  the  sphere  merely 
affects  the  overall  spectral  distribution  of  responsivity  which,  thus  may  be  slightly 
different  from  that  for  the  spectroradiometer  used  alone  as  in  problem  1. 

Again,  as  in  problem  1,  we  first  write  a  general  expression,  in  terms  of  the  flux 
responsivity  R. ,   for  the  signal  element  dS  produced  by  an  element  of  flux  d$  incident 
on  the  receiving  aperture  (entrance  port  —  see  figure  5.7) 

dS(x,y,6,<}>,A,Ao)   =   V<1$ 

=  R.(x,y,6,<J>,A,A  )«L..  (x,y,0  ,<|>,A)  «cose -dco-dA'dA  [S].    (5.4) 

$  O    A 

And,  again,  the  measurement  equation  is  obtained  by  integrating  eq.  (5.4)  to  obtain  the 
expression  for  the  total  output  signal 

S(A,oj,AA,A  )   =  /   /  /  R$*Lx«cos6'dwdA«dA  [S]  ,  (5.5) 

AA  A  to 

where,  now, 

u  [sr]   is  the  acceptance  solid  angle  that,  at  any  point  x,y  in  the 
receiving  aperture  (entrance  port  of  sphere) ,  includes  all  converging 
incident  rays  to  which  the  instrument  responds, 

2 
A  [m  ]   is  the  area  of  the  receiving  aperture  (entrance  port  of  sphere) , 

and 

AA  [nm]   is  the  wavelength  pass  band  of  the  entire  instrument,  including 

all  wavelengths  for  which  the  spectral  flux  responsivity  is  significantly 

non-zero;  i.e.,  where  R^(A  ,A)  ^   0. 

$  o 

The  first  step  in  solving  this  measurement  equation  is,  again,  to  introduce  the  radio- 
metric quantity  to  be  measured  (in  terms  of  the  incident  spectral-radiance  distribution). 
The  incident  spectral  irradiance  at  the  point  x,y  and  the  wavelength  A   [see  eq.  (4.4) 
and  the  last  section  of  Chapter  4  entitled  "Spectral  Radiometric  Quantities"]  is  given  by 

Ex(x,y,A)  =     d$x(x,y,A)/dA 

=  /  L.(x,y,8,<fr,A)'cos8'da>  [W*m  ].  (5.21) 


A 

01 
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Lx<x,y,0,*,X) 


X 


A  —  circular  area  of  receiving  aperture  (entrance  port 
of  averaging  sphere) 

x  ,y  —  center  point  of  A 
o'-'o  r 

x,y  —  any  point  within  A 

L.  (x,y,6  ,((i,A)  —  spectral  radiance  incident  at  x,y  from 

A 

direction   6,<|>   at  wavelength  A 


Figure  5.7.    Designation  of  spectral-radiance  distribution  incident 
at  receiving  aperture 
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In  eq.  (5.5),  if  R   were  brought  outside  the  inner   integral,  that  integral  would  be 
identical  to  that  of  eq.  (5.21)  [see  eq.  (5.22)].   But  we  have  already  seen  that  the  effect 
of  the  averaging  sphere  is  to  make  the  flux  responsivity  independent  of  the  position  and 
direction  of  the  incident  rays  within  A  and  w.   Accordingly,  we  can  write 
R  (x,y,6,(j),  A,  A  )  =  R4(A  ,  A)   so  that  eq.  (5.5)  becomes 

S(A,u),AA,A  )   =   /   /   /  LA(x,y,e,<t>,A)*cose-du)  -R$(Ao, A) «dA-dA  [S],       (5.22) 

AA  A  a) 

and  we  can  replace  the  inner  integral  (in  large  parentheses)  with  E  (x,y,A)   [eq.  (5.21)] 

A 

to  give 

S(A,AA,A  )   =   /   /  EA(x,y,A)-R  fA  ,A)'dA-dA  [S].  (5.23) 

O  .  ,  '    A  <V        O 

AA  A 

Actually,  when  it  is  clear  that  the  fluorescent  lamp  is  the  only  significant  source  irradi- 
ating the  receiving  aperture,  the  incident  spectral  radiance  L.   at  any  point   x,y   in 
that  receiving  aperture  is  zero  for  all  directions  outside  the  solid  angle   w   subtended 
at  that  point  by  the  lamp.   In  that  case  the  integration  in  eq.  (5.22)  need  be  carried  out 

a 

only  over  oo   rather  than  over  the  full  acceptance  solid  angle  w,   so  that 

E^(x,y,A)   =   /  L  (x,y,0,({i,A) 'cose-dca  [W'm~2-nm-1].  (5.24) 

Furthermore,  if  the  lamp  distance  is  large  enough  in  relation  to  the  size  of  the  receiving 

£ 
aperture,  neither  the  size  of  the  solid  angle  oj   nor  the  value  of  the  incident  spectral 

radiance  L   from  any  point  of  the  lamp  will  change  significantly  as  a  function  of  position 

A 

x,y  across  the  receiving  aperture.   In  other  words,  the  spectral  irradiance  will  be  con- 
stant, and  equal  to  its  value  at  the  center  x  ,y   across  the  entire  aperture 

E*(x,y,A)   =   E*  (xo,yQ,A)  [W-m"2 -nm"1] .  (5.25) 


If  we  also  assume,  for  the  moment,  that  the  spectral  irradiance  does  not  vary  significantly 

with  ^ 

that1 


with  wavelength  A  within  the  pass  band   AA   for  any  given  setting   A  ,   we  may  also  write 


E*(x,y,A)   =   E*(x  y   A  )  [W-m  2 -nm  1],  (5.26) 

A  A   O   o   O 

which  may  be  treated  as  a  constant  and  brought  outside  the  integrals  in  eq.  (5.23),  leaving 


tacitly  assuming  that   A   lies  within   AA.   This  may  not  be  true  for  a  nominal  instru- 
ment "setting"  but  will  hold  for  the  calibrated  or  "corrected"  value  corresponding  to  that 
"setting". 
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S  (A,AA,A  )   =  A-f  E;(x  ,y  ,A)-RA(A  ,A)-dA 
o        J  . ,  A  o  •'o     $  o 

AA 

■  ea(wV*a*'  WA),dA 

AA 

=  ea<wV'Vxo),a*ax  cs]-  (5-27) 

The  average  responsivity  RA(A  )  =" (1/AA) •/  R.(A  ,A)«dA  can  be  eliminated  along  with  A 

AA    " 
and  AA,   by  making  a  similar  measurement  of  the  spectral  irradiance  produced  in  the 

receiving  aperture  by  a  tungsten-halogen-lamp  standard  of  spectral  irradiance.   The  corre- 
sponding equation  for  the  standard  lamp  can  be  derived,  step  by  step,  in  the  same  manner  as 
for  the  fluorescent  lamp.   The  result,  corresponding  to  eq.  (5.27),  is 

SS(A,AA,A  )   =   E^(x  ,y  ,A  ) -r7(A  )-A-AA  [S].  (5.28) 

o       A  o  J  o     o   $  o 

By   combining  eqs.    (5.27)    and    (5.28),   we  obtain 

E^(xo,yo,Ao)      =      (S£/SS)-E^(xo,yo,Ao)    [W-m"2-nm_1],  (5.29) 

the  simplest  relation  for  measuring  the  spectral  irradiance  of  the  fluorescent  lamp  in 
terms  of  that  of  the  standard  lamp. 

Note  that,  since  we  are  concerned  with  irradiance  at  the  reference  surface  across  the 
entrance  port  of  the  sphere,  we  have  not  introduced  the  propagance   x*   as  we  did  in  prob- 
lem 1,  where  we  were  measuring  L.   at  the  lamp  source.   However,  if  the  propagance   x* 
over  the  relatively  short  path  in  the  laboratory  for  the  standard-lamp  measurement  should 
be  significantly  different  from  that  in  the  standards  laboratory  where  that  lamp  was  cali- 
brated, this  could  produce  an  error  in  our  results. 

The  most  critical  condition  for  this  measurement  is  the  constancy  of  directional 

responsivity  over  the  solid  angle  of  acceptance  to,   required  to  go  from  eq.  (5.5)  to  eq. 

(5.22).   Actually,  as  we  can  see  from  eq.  (5.24)  it  is  sufficient  if  the  constant  direc- 

l 
tional  responsivity  holds  only  over  the  smaller  solid  angle  w  ,   as  well  as  over  the  still 

g 
smaller  w   for  the  corresponding  treatment  of  the  standard  lamp,  providing  no  other 

sources  are  present.   Fortunately,  as  previously  mentioned,  averaging  spheres  that  perform 
well  in  this  respect  are  available.   Also,  it  is  not  difficult  to  verify  this  performance 
if  a  sufficiently  collimated  beam,  wide  enough  to  more  than  fill  the  entrance  port  of  the 
averaging  sphere,  is  available.   When  the  angle   6  between  the  collimated-beam  direction 
and  the  normal  to  the  plane  of  the  entrance  port  is  varied,  the  instrument  output   S 
should  vary  in  direct  proportion  to  cos6   (a  so-called  lambertian  responsivity  character- 
istic).  This  is  readily  apparent  in  eq.  (5.5)  if   oj  is  made  small  (a  small  distant 
source)  and  all  other  quantities  on  the  right-hand  side  of  the  equation,  particularly  R  , 
remain  constant  while   6   is  varied. 

H 
Equation  (5.25)  was  based  on  two  assumptions:   (1)  that  the  solid  angle  w   (subtended 

by  the  lamp)  does  not  change  significantly  as  a  function  of  position  x,y  across  the 
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receiving  aperture,  and   (2)  that  the  incident  spectral  radiance  L.   from  any  point  of  the 
lamp  also  does  not  change  significantly  with  x,y.   If  we  set  the  height  of  the  lamp  z 
above  x  ,y   equal  to  the  length  of  the  lamp,  and  the  receiving-aperture  diameter  approxi- 
mately equal  to  0.01*z  ,   it  is  easily  shown  that  the  angle  subtended  by  the  long  dimension 
of  the  lamp  will  change  by  much  less  than  a  tenth  of  one  per  cent  between  the  center  x  ,y 
and  the  edge  of  that  aperture.   The  change  is  obviously  much  less  for  the  shorter  dimension 
of  the  fluorescent  lamp  and  for  the  smaller  tungsten-halogen  lamp.   On  the  other  hand,  the 
possible  variation  of  lamp  spectral  radiance  L.   with  direction  from  each  point  of  the 
lamp  may  not  be  so  easy  to  measure  or  evaluate.   However,  the  constancy  of  irradiance 
across  the  receiving  aperture,  assumed  in  eq.  (5.25),  can  be  checked  directly  by  placing 
over  it  a  screen  with  the  smallest  possible  hole  that  allows  just  enough  flux  to  pass  to 
give  a  measurable  output  and  observing  the  constancy  of  that  output  as  the  hole  is  moved 
about  across  the  receiving  aperture.   If  necessary,  the  spectroradiometer  slits  can  be 
widened  to  increase  AX  and  pass  more  flux  to  produce  a  stronger  signal  output  for  this 
purely  geometrical  check,  since  there  is  no  reason  to  believe  that  the  directional  distri- 
bution  will  be  significantly  different  at  slightly  different  wavelengths. 

Finally,  we  have  assumed,  in  going  from  eq.  (5.25)  to  eq.  (5.26)  that  E.   is  constant 
with  respect  to  X  within  AX.   For  thermal-radiation  sources,  such  as  tungsten  lamps  and 
blackbody  simulators,  it  is  usually  possible  to  make  the  spectral  pass  band  AX  of  the 
spectroradiometer  narrow  enough  so  that  E.   does  not  change  by  a  significant  amount  within 
that  spectral  interval.   However,  when  there  are  spectral  lines,  as  in  the  radiation  from 
a  fluorescent  lamp,  or  where  thermal  radiation  falls  off  rapidly,  particularly  at  the 
short-wavelength  end  of  the  spectral  curve  of  thermal  emission,  there  may  be  large  varia- 
tions in  E.   within  the  interval  AX.   We'll  have  to  leave  much  of  the  treatment  of  that 
X 

problem  for  a  later  chapter,  after  the  concept  of  slit  function  has  been  developed.   For 
the  present,  we  recognize  that  when  AX  is  large  enough  so  that  there  are  significant 
variations  in  the  spectral  irradiance  E.(X)   over  the  spectral  pass  band  AX,   the  situa- 
tion is  essentially  no  different  than  that  of  a  broad-band  radiometer  with  a  non-uniform 
spectral  responsivity,  which  we  take  up  next  in  problem  3.   Accordingly,  all  of  the  dis- 
cussion in  problem  3  (following)  is  equally  applicable  here.   As  will  be  clear  from  that 
discussion,  there  are  always  problems  when  both   the  incident  spectral  distribution  and  the 
spectral  responsivity  vary  significantly  across  the  pass  band.  With  a  monochromator  there 
is  always  large  Variation  of  the  overall  instrument  spectral  responsivity,  about  which  we 
will  learn  more  in  a  later  chapter. 

Problem  3.   Determining  irradiance  with  a  broad-band  radiometer. 

A  practical  example  of  this  problem  is  determining  the  irradiance  (incident  watts  per 
square  meter)  between  400  and  500   [nm]  produced  at  the  treatment  distance  by  a  bank  of 
blue  lamps  used  in  the  phototherapy  of  infant  jaundice.   A  schematic  diagram  of  such  an 
experimental  setup  is  shown  in  figure  5.8.   The  measurement  configuration  is  similar  to 
that  of  problem  2  (figure  5.4)  except  that  a  filter  radiometer  that  responds  to  incident 
radiation  in  a  broader  spectral  pass  band  is  used  in  place  of  the  spectroradiometer.  We 
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RADIOMETER    RECEIVING 
APERTURE   (with  diffuser) 


BROAD-BAND 
RADIOMETER 


Figure  5.8.   Sketch  of  parts  of  configuration  for  measurement  of  irradiance 
with  broad-band  (broad-spectral-band)  radiometer. 


78 


still  need  spatially  uniform  responsivity,  so  some  form  of  diffuser,  such  as  an  averaging 
sphere,  is  again  used  to  achieve  uniform  flux  responsivity  to  all  rays  of  any  given  wave- 
length incident  on  the  receiving  aperture  from  a  wide  range  of  directions.   The  wider 
spectral  pass  band  is  now  determined  primarily  by  a  spectral  filter  in  conjunction  with  the 
spectral  responsivity  of  the  detector  or  transducer.   The  spectral  characteristics  of  other 
portions  of  the  internal  optical  path,  including  the  diffuser,  may  also  contribute  to  the 
varying  overall  spectral  responsivity. 

The  measurement  equation  for  this  problem  is  obtained  in  the  same  manner  as  that  for 
problem  2,  and  is  identical  to  it  [eq.  (5.5)]  except  that  the  spectral  pass  band  AA  is 
larger,  now  of  the  order  of  100  [nm] ,   and  there  is  no  wavelength  "setting"  A  : 

S(A,w,AA)   =  /   /  /  R  •Lx-cos6«dwdA«dA  [S]  ,  (5.30) 

AA  A  a) 

where  now 

to  [sr]   is  the  acceptance  solid  angle  of  the  instrument  at  the  point 
x,y  within  its  receiving  aperture, 

2 
A  [m  ]   is  the  area  of  the  receiving  aperture,  and 

AA  [nm]   is  the  spectral  pass  band  of  the  instrument,  including  all 
wavelengths  A  for  which  the  spectral  flux  responsivity  is  signif- 
icantly non-zero  —  R^^)  "*  0.   [This  is  determined  by  the  combined 
spectral  characteristics  of  the  detector  and  all  internal  path(s) 
and  optical  elements (s),  including  the  diffuser  and,  especially 
any  spectral  filter(s).] 

As  in  problem  2,  it  is  essential  that  the  flux  responsivity  R   be  constant  every- 
where within  the  acceptance  solid  angle  w,   or  at  least  as  much  of  that  solid  angle  as  is 
filled  by  rays  of  sufficient  radiance  to  make  a  significant  contribution  to  the  output 
signal  of  the  instrument.   If  there  are  no  other  sources  present,  it  is  the  portion  of  tn 
filled  by  rays  from  the  bank  of  lamps  being  measured,  including  reflections  and  scattering 
from  the  fixture  in  which  they  are  mounted.   Then,  as  before,  the  responsivity  can  be  taken 
outside  of  the  inner  integral,  the  integral  with  respect  to  solid  angle,  so  that  that 
integral  gives  us  the  spectral  irradiance.   Then  the  equation  becomes  [eq.  (5.23] 

S(A,AA)   =  /   /  E  (x,y,A)-R  (A)-dA'dA  [S].  (5.31) 

AA  A 

Again,  as  before,  by  a  suitable  choice  of  the  size  of  the  receiving  aperture  in  relation  to 
the  working  (treatment)  distance  plus  adequate  experimental  verification,  we  may  treat  E, 
as  well  as  R.   as  constant  for  all  points  x,y  over  the  aperture  area  A,   so  that 

S(AA)   =  A-J   R$(A)-EA(A)-dA  [S] .  (5.32) 

AA 
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Although  eq.  (5.32)  appears  to  be  a  relatively  simple  equation,  with  only  one  variable 
parameter,  the  wavelength   A,   it  is  a  remarkably  troublesome  one  and  it  represents  a 
situation  that  arises  again  and  again  in  radiometry.   So  we'll  take  a  particularly  careful 
look  at  it.   Note  that  the  quantity  that  we  want  to  measure  in  this  case  is  the  irradiance 

500 
E  =   /   E  (A)-dA  [W«m  ],  (5.33) 

400 

the  integral  from  400   to  500  [ran]   of  the  spectral  irradiance  E  (A).   Unfortunately, 

A 

however,  eq.  (5.31)  does  not  give  us  a  measurement  of  this  desired  quantity  E   [of  eq. 
(5.33)]  because  the  integral  in  eq.  (5.32)  is  a  weighted  integral  of  the  spectral  irradiance, 
with  the  spectral  flux  responsivity  R(fi(A)   as  tne  weighting  function,  and  AA  can  seldom 
be  made  to  correspond  exactly  to  the  desired  wavelength  interval,  in  this  case  400-500 
[nm]. 

If  we  were  so  fortunate  as  to  have  an  instrument,  with  uniform  spectral  responsivity1 
R(A)  =  R  over  its  entire  pass  band   AA  =  A  -  A..   and  zero  responsivity  R(A)  =  0  outside 
of  that  band  [for   A  <  A..   and   A  >  A„],   and  with  the  limits  exactly  matching  the  desired 
interval   (A  =  400  [nm]   and   A  =  500  [nm]),   as  depicted  by  the  dashed  lines  in  figure 
5.9,  there  would,  of  course,  be  no  problem.   Then  eq.  (5.32)  simplifies  to 

500 
S(AA)   =  A-R  -J   E  (A)-dA  =  A-R  -E  [S],  (5.34) 

400 

and  we  need  only  make  a  similar  calibration  measurement,  with  a  known  standard  source 
substituted  for  the  unknown,  to  also  obtain 

500 
SS(AA)   =   A-R  •/   Ef(A)-dA   =  A-R  -ES  [S].  (5.35) 

$  400  A  * 

We  can  then  combine  eqs.  (5.34)  and  (5.35)  to  obtain 

500 
E  =   (S/SS)-ES   =   (S/SS)-J   E^(A)-dA  [W-m  Z],  (5.36) 

400  A 

g 
which  can  be  evaluated  from  the  values  of  E  (A)   furnished  with  the  spectral-irradiance 

A 

standard.   But  such  an  ideal  instrument  response  is  very  difficult,  if  not  impossible,  to 
approximate  adequately. 

Now  let's  consider  instruments  with  non-uniform  spectral  responsivity.   In  figure  5.9, 
in  addition  to  the  dashed  lines  showing  the  ideal  "rectangular"  responsivity  just  discussed, 
the  solid  curve  shows  a  more  realistic  spectral-responsivity  function  for  a  broad-band 


•^Here  R  may  be  flux  responsivity   R  ,   irradiance  responsivity  R^  =  A-R  ,   [see  below, 
eq.  (5.49)],  etc.   It  is  R,(A)  =  R   that  we  use  in  eq.  (5.34). 
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X  [nm] 


Figure  5.9.   Spectral  responsivity  R(A)   of  broad-band  radiometer. 

The  dashed  lines  represent  the  ideal  spectral 
responsivity  R(A)   desired  for  problem  3. 


radiometer.   The  responsivity  is  not  constant  over  all  of  the  pass  band  and  the  pass  band 
does  not  have  sharp  limits  outside  of  which  the  responsivity  vanishes  completely.   In  such 
cases,  we  can  still  obtain  the  desired  value  of  measured  irradiance  E,   given  by  eq. 
(5.33),  if  we  know  both   the  relative   responsivity  r(A)  [dimensionless]  and   the  relative 
distribution  of  incident  spectral  irradiance  e.  (A)  [nm  ],  defined,  respectively,  by 


R$(A)   5  K^rCA)  [S-W  1] , 


(5.37) 


and 


EX(A)  =     K2-ex(A)  [W-m  2-nm  1] , 


(5.38) 


,-lr 


-2. 


where  K-  [S«W  ]   and  K~  [W-m  ]   are  both  constants  with  respect  to  wavelength.1   Then, 
when  we  substitute  from  eq.  (5.37)  and  (5.38)  into  eq.  (5.32),  we  have 


1We  are  not,  at  this  point,  addressing  the  question  of  how  to  determine  either  r(A)   or 

e..  (A)   (which  can  be  a  very  difficult  problem).  All  we're  concerned  about  now  is  the  con- 

A 

cepts  involved  when  we  do   have  at  least  this  much  information. 
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X2 
S(AX)   =  A'K^K^-J  r(A)-eA(X)-dX  [S],  (5.39) 

Xl 

where,  now,  the  limits  of  AX  (5X.  -  X.)   extend  below  400  [ran]   and  above  500  [ran]   to 
include  all  wavelengths  where  R  (X)  ^  0.   In  figure  5.9  it  is  evident  that,  at  both  limits, 
the  responsivity  approaches  zero  gradually  so  that  exact  values  for  X1   and  X„  are  not 
clearly  established.   The  values  assigned  to  them  will  depend,  in  most  cases  on  the  noise 
level  of  the  measuring  equipment.   In  addition,  as  suggested  earlier,  this  choice  also 
depends  on  the  incident  radiation.   The  integrated  contribution  of  very  strong  incident 
radiation  at  wavelenghts  where  R^C^)   is  supposed  to  be  zero  but  is  actually  only  very 
small  can  be  significant  while  ordinarily,  with  normal  levels  of  spectral  irradiance  at 
these  wavelengths,  there  is  no  observable  contribution.   Accordingly,   X..   must  be  chosen 
small  enough  and  X„  large  enough  to  include  all  significantly  non-zero  values  of  R*(^) 
in  the  given  situation. 

For  a  corresponding  measurement  of  a  standard  of  spectral  irradiance,  we  similarly 
substitute  only  from  eq.  (5.37)  into  eq.  (5.32),  giving 

A2 
SS(AX)   =  A-K  •/   r(X)-E^(X)-dX  [S].  (5.40) 

Al 
Next,  we  divide  eq.  (5.39)  by  eq.  (5.40)  to  obtain 

g  /  r(X).eA(X)-dX 

=  K0-  — - ,  (5.41) 


,s  r   /,s  ^s 

JX 


/   r(X)-E^(X)«dX 

AX 


from  which  we  can  evaluate  K„   as 


S 

/     r(X)-E^(X)-dX 

AX 

[W«m  2]. 

2 

ss 

/     r(X)-e   (X)-dX 

AX 

(5.42) 


Finally,  combining  eq.  (5.33),  (5.38),  and  (5.42),  the  result  is 


500 
E  =  K  •/   e  (X)-dX 
400 


/   r(X).E^(X).dX    50Q 
AX  r     /,n  ,,  rT,  -2- 


/   e  (X)-dX  [Win   ],  (5.43) 

SS     /   r(X)-e  (X)-dX    400 

AX 
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where  all  of  the  quantities  on  the  right-hand  side  of  the  equation  are  known. 

In  practice,  accurate  measurements  of  irradiance  with  a  broad-band  radiometer  require 
the  solution  of  eq.  (5.43)  because  instruments  with  ideal  "rectangular"  spectral  responsiv- 
ities  are  just  not  available.   If  the  relative  spectral  responsivity  r(A)   of  the  radiom- 
eter and   the  relative  spectral  distribution  ei(^)   of  the  irradiance  being  measured,  which 
are  both  required  for  a  solution  of  eq.  (5.43)  are  not  known,  one  must  resort  to  narrow- 
band measurements  using  a  suitable  spectroradiometer.   In  this  case  the  spectral  irradiance 
is  measured  as  a  function  of  wavelength  and  the  desired  irradiance  is  obtained  by  integra- 
tion over  the  wavelength  interval  of  interest,  as  in  eq.  (5.33). 

NORMALIZATION  of  BROAD-BAND  MEASUREMENT  RESULTS.    A  quantity  that  is  occasionally  useful 
and  that  can  be  obtained  with  a  broad-band  radiometer,  when  only  the  relative  spectral  flux 
responsivity  r(A)   of  the  radiometer  is  available,  is  a  lower  bound  for  the  value  of 

E(AA)   =   /   E  (A)'dA  [W«m~2],  (5.44) 

AA 

the  irradiance  over  the  spectral  band  AA  =  A  -  A.,   of  the  radiometer.1   This  can  be  done 
by  setting  the  peak  relative  responsivity  r(A  )  =  1,  where  A   is  the  wavelength  at 
which  the  spectral  responsivity  has  its  maximum  or  peak  value.   This  is  equivalent  to 
setting  the  constant  IL   in  eq.  (5.37)  equal  to  that  maximum  responsivity:  K-  =  R*(A  ). 
At  all  other  wavelengths  r(A)  £  1,   so  we  can  write 

/  r(A)-E  (A)-dA   <;  /  E  (A)-dA.  (5.45) 

AA  AA 

This  is  often  referred  to  as  peak  normalization  and  the  quantity  on  the  left-hand  side  of 
eq.  (5.45)  is  called  the  peak-normalized  value  of  incident  irradiance  E  .  Accordingly, 
eq.  (5.45)  can  be  combined  with  eq.  (5.44)  and  rewritten  as 

E   (AA)  <;   E(AA)  [W«m~2],  (5.45a) 

np 

where 

En  (AA)   E   /   r(A)-Ex(A)-dA,    for   r(A)  <;  r(A  )  =  1, 

AA 

=  K2-/  r(A)-eA(A)-dA.  (5.46) 

AA 

With  this  last  expression,  eq.  (5.46),  it  clearly  follows  from  eq.  (5.41)  that 


*E(AA)   given  by  eq.  (5.44)  is  equal  to  E  of  eq.  (5.33)  only  when  \^  =   400  [nm]   and 
A„  =  500  [nm];  i.e.,  when  the  instrument  has  exactly  the  desired  pass  band. 
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E   (AX)   =  -2—  •/   r(X)-E^(X)-dX  [W-m  2],  (5.47) 

np  SS    AX      X 

which  can  be  evaluated  because  all  factors  on  the  right-hand  side  of  the  equation  are 

known.   Thus,  our  measurement  gives  a  value  for  the  peak-normalized  irradiance  E   (AX) 

np 

which,  by  eq.  (5.45a)  is  a  lower  bound  for  the  irradiance  E(AX) .   In  addition,  to  the 
extent  that  the  spectral  responsivity  Ra(^)   approximates  the  rectangular  dashed  "curve" 
in  figure  5.9,  with  limits  at  400  and  500  [nm] ,   E   (AX)   is  an  approximation  to  the 
lower  bound  for  the  desired  irradiance  E  of  eq.  (5.33). 

There  are  other  "normalization"  practices  where,  for  various  reasons,  different  choices 
are  made  for  the  value  of  K.,   which,  as  we  have  seen,  can  be  assigned  arbitrarily.   We 
won't  go  into  all  the  details  of  this  rather  complex  subject  here.   There  is  a  very  general 
treatment  [21]  for  anyone  who  wishes  to  dig  more  deeply. 

When  actual  values  can't  be  obtained,  normalized  values  usually  are  reported  in  order 
to  try  to  obtain  results  of  some  value  from  such  broad-band  measurements  of  incident  radia- 
tion of  unknown  spectral  distribution,  made  with  instruments  of  non-uniform  spectral  respon- 
sivity.  Unfortunately,  the  normalization  methods  used  are  not  always  clearly  defined  or 
specified,  which  defeats  their  purpose.   Also,  repeated  use  of  terms  applied  to  some  of 
them  seems  to  build  up  confidence,  based  mostly  on  familiarity,  that  blinds  many  people  to 
their  limitations.   There  just  isn't  any  way  to  get  around  the  fact  that,  without  the 
incident  relative  spectral  distribution,  quantities  measured  in  lumens,  roentgens,  the 
"effective  watts"  used  by  many  in  the  military- infrared  community,  etc.,  and  other  similar 
quantities,  can't  be  converted  directly  to  the  actual  quantities  in  absolute  units.  Nor- 
malized quantities  (spectral-responsivity-weighted-integral  quantities)  can  be  directly 
related,  in  general,  only  to  measurement  results  from  instruments  having  the  same  relative 
spectral  responsivity  r(X).   When  they  are  combined  with  approximate  information  about  the 
incident  spectral  distribution,  they  can  provide  correspondingly  approximate  values  or 
estimates  of  the  true  irradiance  (or  other  radiometric  quantity) .   The  strong  feelings  of 
confidence  that  are  often  attached  to  such  normalized  units  and  quantities  probably  stem 
from  the  development  of  a  faculty  for  making  such  estimates,  based  on  estimates  or  assump- 
tions about  the  incident  spectral  distribution  of  which  the  individual  may  not  even  be 
consciously  aware. 

We  note,  also,  in  passing,  that  similar  considerations  arise  in  connection  with  non- 
uniform responsivity  with  respect  to  any   of  the  radiation  parameters,  not  just  the  spectral 
parameter.   For  example,  many  instruments  do  not  have  a  uniform  directional  responsivity, 
particularly  near  the  edges  of  the  field  of  view  which  is  often  not  very  sharply  defined. 
When  that  is  the  case,  the  output  signal  is  again  a  responsivity-weighted  integral  of  the 
incident  radiance  distribution  and  similar  normalization  considerations  are  applicable, 
since  the  mathematical  relationships  are  exactly  parallel  to  those  in  the  spectral  case 
just  analyzed.   This  is  all  covered  in  the  general  treatment  [21]. 
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RESPONSIVITY  CALIBRATIONS.    There  are  many  situations  where  the  user  of  a  radiometer  or 
radiometric  instrument  may  not  make  a  calibration  measurement  or  measurements  each  time  he 
makes  a  measurement  of  an  unknown  radiometric  quantity.   This  may  be  for  a  number  of 
reasons  that  we  won't  attempt  to  analyze  except  to  point  out  that,  when  substantial  periods 
of  time  elapse  between  a  measurement  and  the  related  calibration,  care  should  be  taken  to 
make  sure  that  the  instrument  stability  is  adequate  to  meet  the  needed  measurement  accuracy. 
In  any  case,  when  this  is  done,  say,  for  a  broad-band  measurement  of  irradiance,  as  in 
problem  3,  the  calibration  measurement  [eq.  (5.35)]  is  then  treated  as  a  measurement  to 
establish  the  flux  responsivity  of  the  instrument  as 

R$  =   SS/(A-E3)  [S-W-1]  (5.48) 

or,  even  more  directly,  the  irradiance  responsivity  as 

RE  =   SS/ES  =  A-R$  [S-W_1-m2].  (5.49) 

This  value  of  responsivity  may  then  be  used  to  evaluate  the  measured  values  of  irradiance 
from  unknown  sources  in  subsequent  measurements  as 

E  =   S/R^  =   S/(A«R  )  [W-m-2]  •  (5.50) 

as  long  as  the  instrument  stability  is  adequate  to  meet  the  desired  accuracies. 

GENERAL  DISCUSSION  of  the  MEASUREMENT  EQUATION  and  ITS  USE.   We  have  introduced  the  meas- 
urement equation  by  deriving  it  and  obtaining  its  solution  for  three  rather  common  radio- 
metric problems.   For  each  problem,  the  equation  turned  out  to  have  the  same  form,  namely 

S  =  /  /  /  R  »L  •cose-dwdA'dX  [S].  (5.30a)1 

AA  A  a) 

Actually,  this  is  what  one  would  expect,  because  the  equation  is  simply  a  summing  up  (inte- 
grating) of  all  the  elements  of  radiometer  output  signal 

dS  =  R  -L  -cose-dwdA-dA  (5.4a) 

$   A 

produced  by  all  the  elements  of  radiant  flux 

d$  =  L  •cosO'doi'dA'dA  [W]  (3.10a) 

A 

that  are  incident  on  the  receiving  aperture  of  the  radiometer.  These  statements  are  quite 
general;  accordingly,  eq.  (5.30a)  is  applicable  to  any  optical  radiation  measurement  prob- 
lem involving  incoherent  radiation. 

Unfortunately,  there  is  no  unique  general  solution  to  this  measurement  equation.   Even 


^his  is  identical  to  eq.  (5.30)  except  that  functional  dependence  of   S   is  not  shown 
explicitly. 
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if  the  spectral  responsivity  R  (x,y,6,<j>,  A)   is  completely  specified  as  a  function  of 
position  x,y   (on  A),  direction  6,<f>   (in  w) ,   and  wavelength  A   (in  AX),   there  are 
an  unlimited  number  of  spectral  radiance  functions  L  (x,y,6  ,<f>,A)   that  would  produce  the 
same  observed  output  signal  S.   Moreover,   L^   in  eq.  (5.30a)  can't  be  replaced  by  a 

desired  radiometric  quantity  that  is  an  integral  of  L..   (e.g.,   E.  ,   E,   or  F  )   unless 

A  At 

the  equation  can  be  modified  to  .contain  that  integral.   This  requires  that   R   and  possibly 
other  quantities  in  the  measurement  equation  be  constant  with  respect  to  one  or  more  radia- 
tion parameters.   These  difficulties  are  inherent  in  the  multidimensionality  of  optical 
radiation.   In  principle,  it  is  always  possible   to  use  a  spectroradiometer  with  small 
enough  AA,  A,   and  w  so  that  L   is  sufficiently  constant  over  each  of  those  integra- 
tion intervals  that  it  can  be  removed  from  all  the  integrals.   Then  measurement  of  an 
unknown,  and  of  a  known  standard,  allow  one  to  obtain  L^ .   If  this  is  done  for  all  possible 
combinations  of  the  radiation  parameters,  any  radiometric  quantity  desired  can  then  be 
calculated.   This,  however,  is  usually  an  extremely  difficult  and  tedious  task  and  is 
generally  impractical. 

Accordingly,  the  practical  solution  is  usually  to  try  to  select  a  measuring  instrument 
and  a  measurement  configuration  to  satisfy  certain  conditions,  at  least  within  a  desired 
degree  of  approximation,  that  will  make  it  possible  to  modify  the  measurement  equation  to 
include  the  desired  radiometric  quantity  and  to  obtain  a  unique  solution.   The  kinds  of 
conditions  we  are  talking  about  are  illustrated  in  the  three  sample  problems.   They  include 
the  constancy  or  approximate  constancy  of  a  radiation  quantity  or  of  the  responsivity 
relative  to  one  or  more  radiation  parameters  so  that  these  functions  can  be  brought  outside 
one  or  more  integrals,  the  use  of  an  average  to  replace  one  or  more  of  the  integrals,  and 
the  use  of  the  relative  spectral  distribution,  if  known,  of  the  otherwise  unknown  radio- 
metric quantity  being  measured.   The  major  advantages  of  using  the  measurement  equation  are 
that  it  makes  clear  that  such  conditions  must  be  sought  and  provides  insight  and  a  system- 
atic approach  towards  finding  them.   In  fact,  without  such  an  approach,  we  believe  it  is 
highly  unlikely  that  one  can  make  state-of-the-art  measurements,  or  even  less  accurate 
measurements,  with  a  meaningful  estimate  of  the  uncertainty. 

To  assist  the  reader  in  making  use  of  the  measurement  equation,  we  have  prepared  a 
list  of  steps  to  follow.   This  is  not  a  comprehensive  list,  however,  because  there  are  two 
more  radiation  parameters  not  yet  treated  as  well  as  other  concepts  still  to  be  developed 
(see  next  section).   In  later  chapters,  the  kinds  of  problems  that  can  be  addressed  will  be 
extended  and  alternatives  in  approach  will  be  provided.   In  fact,  the  entire  Self-Study 
Manual  can  be  thought  of  as  a  treatment  in  which  concepts  and  techniques  are  developed  for 
use  in  the  measurement  equation  and  expertise  in  its  use  is  developed  through  discussion 
and  analysis  of  various  applications.   Nevertheless,  the  steps  that  are  listed  below  are  a 
useful  beginning  and  we  recommend  that  they  be  followed  routinely  in  applying  the  measure- 
ment equation. 
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1.  Prepare  a  diagram,  approximately  to  scale,  of  the  entire  measurement 
configuration,  showing  all  pertinent  radiation  sources,  both  the 
source (s)  of  the  radiation  beam  to  be  measured  and  any  source (s)  of 
unwanted  or  background  radiation,  including  objects  from  which 
radiation  might  be  scattered  or  reflected  into  the  measuring  instru- 
mentation.  Also  show  the  measuring  instrumentation  in  sufficient 
detail  to  establish  the  features  required  in  the  following  steps. 

2.  Place  the  calibration  standard  that  is  to  be  used  on  this  same  diagram. 

3.  Determine  the  approximate  size  and  position  of  the  receiving  aperture 
A  and  the  acceptance  (field)  solid  angle  u)  of  the  instrumentation 
and  show  them  on  the  diagram.   Show  a  few  typical  rays  within  them  as 
well  as  the  extreme  rays  that  indicate  their  limits.   These  values 

of  A  and  w,   along  with  AA   from  the  next  step  and  an  estimate 
of  R,   can  be  used  in  the  measurement  equation  to  obtain  approxima- 
tions of  the  expected  output  signals. 

4.  Prepare  a  graph  of  the  relative  spectral  responsivity  of  the  measuring 
instrumentation  as  a  function  of  wavelength,  if  this  information  is 
available,1  showing  the  limits  A   and  A„  of  the  spectral  pass  band 
AA  =  A„  -  A1  .   On  the  same  graph,  plot  also,  if  available,  the  approx- 
imate relative  spectral  radiance  distribution  in  the  incident  beam  that 
is  to  be  measured. 

Taken  together,  such  a  diagram  and  graph(s)  provide  a  rough  picture  that  is  very  helpful  in 
visualizing  the  important  interrelationships  between  the  different  parts  of  the  measurement 
configuration  when  they  are  more  carefully  analyzed  in  later  steps.   It  also  provides  a 
basis  for  preliminary  estimates  of  the  adequacy  of  the  instrumentation  to  make  the  desired 
measurement,  the  appropriateness  of  the  spectral  responsivity  relative  to  the  wavelength 
range  present  in  the  beam  to  be  measured,  the  size  and  configuration  of  the  instrument 
throughput  relative  to  the  geometry  of  the  beam  to  be  measured,  etc.   If  everything  appears 
to  be  in  order,  indicating  that  the  proposed  measurement  is  feasible,  we  can  continue  with 
the  modification  and  simplification  of  the  equation. 

5.  Write  down  the  mathematical  relationship  (from  previous  chapters, 
particularly  Chapter  4)  between  the  radiometric  quantity  to  be 
measured  and  the  spectral  radiance  incident  on  the  receiving 
aperture  of  the  radiometer.   Determine  how  the  measurement 
equation  must  be  modified,  and  whether  such  a  modification  can 
be  carried  out  with  sufficient  accuracy,  so  that  this  relation- 
ship can  be  used  to  replace  the  incident  spectral  radiance  L 


■^Chapter  7  addresses  the  problem  of  determining  the  relative  spectral  responsivity  r(A  ,A) 
and  the  pass  band  AA   of  a  spectroradiometer  containing  a  monochromator . 
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in  the  measurement  equation  by  the  desired  radiometric  quantity. 
An  example  of  this  is  given  in  problem  2  where  it  was  necessary 
to  remove  the  flux  responsivity  R   from  within  the  integral 
over  solid  angle.   This  integral  then  became  the  spectral  irradi- 
ance  E  ,   the  desired  radiometric  quantity.   The  substitution 
could  be  made  only  if  R   had  the  same  value  throughout  w  at 
every  point  of  A.   A  constancy  such  as  this  would  have  to  be 
investigated  and  confirmed  experimentally  in  an  auxiliary  measure- 
ment to  an  accuracy  close  to  that  finally  desired.   If  this 
required  condition  were  not  found  to  exist,  the  instrument  would 
have  to  be  replaced  by  one  in  which  it  did  exist. 

6.  Following  the  introduction  of  the  desired  quantity,  the  resulting 
equation  is  simplified  so  that  a  unique  solution  is  possible. 

As  illustrated  in  the  examples,  this  usually  requires 
that  the  unknown  radiometric  quantity  be  constant 
over  one  or  more  of  the  remaining  integration 
intervals  so  that  it  can  be  removed  from  the 
integrals, 

resorting  to  averages,  or 

utilizing  a  relative  spectral  responsivity  function 
and  a  relative  spectral  distribution  (or  distribution 
with  respect  to  other  radiation  parameters)  of  the 
unknown  radiometric  quantity. 

7.  Repeat  steps  A,  5,  and  6  with  respect  to  the  calibration  standard 
to  be  used. 

8.  Combine  the  simplified  equations  for  the  unknown  and  the  standard 
so  as  to  eliminate  the  responsivity  and  associated  geometrical 
factors  common  to  both  equations.   This  results  in  a  solution  for 
the  unknown  radiometric  quantity  in  terms  of  the  two  measured 
output  signals  (unknown  and  standard)  and  the  known  value  of  the 
calibrated  standard.   Alternatively,  some  prefer  the  mathemati- 
cally equivalent  steps  of  first  solving  the  equation  for  the 
calibration  measurement  of  the  standard,  to  evaluate  the  respon- 
sivity for  the  desired  radiometric  quantity  (combining  the  flux 
responsivity  and  associated  geometrical  factors) ,  and  then  using 
that  responsivity  to  evaluate  the  unknown  quantities  represented 
by  the  output  signals  from  measurements  of  those  quantities.   This 
alternative  approach  is  used  primarily  when  less  frequent  calibra- 
tion measurements  are  made  due  to  less  stringent  accuracy  require- 
ments. 
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An  experienced  radiometrist  will  carry  out  as  much  as  possible  of  this  entire  procedure 
first  as  a  "thought  experiment"  while  planning  the  measurement.   In  this  way  he  may  be  able 
to  determine  whether  or  not  the  measurement  is  feasible  with  the  available  instrumentation, 
or  to  determine  what  instruments  will  be  needed.   It  enables  him  to  plan  ahead  for  the 
calibration  measurements  and  for  any  auxiliary  experiments  that  may  be  required.   Thus,  the 
measurement  equation  provides  a  systematic,  comprehensive  approach  for  planning  how  to  make 
a  particular  measurement  as  well  as  for  carrying  it  out  and  finally  estimating  its  uncer- 
tainty. 

Finally,  when  a  unique  solution  is  not  possible  because  not  even  the  relative  incident 
distribution  is  known,  we  have  seen  that  results  of  limited  usefulness  and  significance  can 
be  obtained  in  the  form  of  normalized  quantities.   Normalization  is  discussed  immediately 
following  the  treatment  of  problem  3 . 

LIMITATIONS  of  the  MEASUREMENT  EQUATION  DEVELOPED  in  THIS  CHAPTER. 

Most  of  the  limitations  have  already  been  pointed  out.   However,  because  of  their  impor- 
tance, they  are  summarized  here.   The  measurement  equation  developed  in  this  chapter  is 
incomplete.   It  does  not   include 

the  radiation  parameters  of  time  and  polarization, 

environmental  and  instrumental  parameters,  such  as  ambient  temperature, 
humidity,  or  magnetic  fields, 

corrections  for  diffraction  effects  [as  a  perturbation  or  departure 
from  pure  geometrical  (ray)  optics],  or 

cases  where  the  responsivity  is  non-linear. 

Moreover,  the  sample  problems  are  stated  as  simply  as  possible  to  illustrate  the  important 
concepts  involved,  and  likely  sources  of  error  in  carrying  out  actual  measurements  are  not 
pointed  out  in  most  cases.   We  intend  to  incorporate  these  additional  parameters  and  cor- 
rections and  to  go  into  greater  detail  in  our  treatment  of  the  measurement-equation  ap- 
proach to  radiometry  in  subsequent  chapters. 

SUMMARY  of  CHAPTER  5.   An  essential  factor  in  the  measurement  equation  is  the  flux  respon- 
sivity, given  here  as  a  function  of  the  spatial  and  spectral  radiation  parameters  only, 
since  we  have  not  yet  taken  up  the  temporal  and  polarization  parameters  nor  the  environmen- 
tal and  instrumental  parameters: 

R$(x,y,e,<f>,A)  B     dS(x,y,e,4>,A)/d$(x,y,e,<t>,A)  tS-w"1].  (5.1a) 

In  terms  of  this  flux  responsivity,  the  output-signal  element  produced  by  an  incident-flux 

element  d$(x,y,6,<(>,  A)  [W] ,   associated  with  an  incident  ray  of  spectral  radiance 

-2   -1   -1 
L^  (x,y,6,<(>,A)   [W*m  *sr   «nm  ],   is 
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dS(x,y,6,cM)   =  R$(x,y,e,<|),A)'d$(x,y,e,<j>,A) 

=  R.  (x,y,e,<(),A)«L,  (x,y,6,(f>,A)  "cose'dwdA'dX  [S]  .  (5.2a) 

When  "reducing"  measurement  data,  some  people  prefer,  instead,  to  use 

d$  =  dS/R$  =  K$-dS  [W],  (5.3) 

where  K  =   -'-/R-*  [W«S   ]   is  sometimes  called  the  calibration  "constant"  although,  like  the 
responsivity,  it  is  obviously  also  a  function  of  all  of  the  radiation  parameters  and  may  be 
quite  variable.   More  generally,  the  responsivity  of  a  detector  element,  or  of  an  entire 
radiometric  instrumentation  system,  is  defined  as  the  output  signal  per  unit  input  of  any 
incident  radiometric  quantity  [17]. 

The  measurement  equation,  relating  the  output  signal  to  the  incident  distribution  of 
spectral  radiance  at  the  receiving  aperture  of  a  radiometric  instrument,  is  introduced  by 
deriving  a  simple  equation,  involving  only  spatial  and  spectral  radiation  parameters,  for 
three  illustrative  problems:   (1)  transferring  the  spectral-radiance  calibration  of  a 
tungsten-ribbon  lamp,   (2)  determining  the  spectral  irradiance  near  a  large  source,  and 
(3)  determining  irradiance  with  a  broad-band  (broad-spectral-band)  radiometer.   It  is 
evident  from  these  examples  that  the  measurement  equation,  in  terms  of  just  the  three 
radiation  parameters  of  position,  direction,  and  spectrum  (wavelength) ,  has  basically  the 
following  form  for  all  radiometric  measurements: 

S(A,w,AA)   =   /   /  /  R$(x,y,0,<J),A)'LA(x,y,e,(j),A)-cose-dca-dA'dA  [S].       (5.30b)1 

AA  A  to 

R  (x,y,0  ,<j>,  A)  [S*W   ]   is  the  instrument  spectral-ray  flux  responsivity  for  the  flux  ele- 
ment 

d$(x,y,6,<j>,A)   =  L,(x,y,6,<l>,A)«cose«dwdA«dA  W  (3.10) 

A 

-2   -1   -1 
associated  with  the  ray,  of  spectral  radiance  L^  (x,y,6  ,<j>,  A)  [W«m   «sr   *nm   ],   incident 

2 
on  the  receiving-aperture- area  element   dA   [m  ]   at  the  point  x,y  within  the  solid-angle 

element  dw  [sr]   from  the  direction   8,<f>   and  within  the  spectral-wavelength  element 

dA  [nm]   at  the  wavelength  A.   The  integration  limits  or  intervals   AA,   A,   and  u 

[more  explicitly,   co(x,y)]   are  selected,  based  on  the  measurement  configuration,  to  include 

all  of  the  incident  radiant  flux  that  might  possibly  contribute  to  the  instrument  response 

(output  signal).   Caution  is  urged  in  excluding  or  ignoring  parameter  values  for  which 

either  R  or  L  ,   alone,  appears  to  be  zero,  particularly  if  it  approaches  zero  gradually 

$       A 

or  asymptotically  as  a  function  of  that  parameter.   Even  though  L.   is  very  small,  the 
integrated  contribution  over  a  wide  wavelength  band  in  which  R   is  substantial  may  be 


^his  is  identical  to  eq.  (5.30)  except  that  the  functional  dependence  of  R   and  L.   on 
the  radiation  parameters  is  explicitly  shown  here. 
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surprisingly  large,  and  the  same  may  be  true  if  R   is  very  small  and  L.   is  substantial 
over  such  a  parameter  level. 

Note  that,  in  addition,  both  S  and  R   may  be  functions  of  a  wavelength  setting 
A  ,   as  in  eqs.  (5.4)  and  (5.5),  when  the  instrument  includes  a  monochromator.   However,  we 
have  given  eq.  (5.30),  in  the  more  explicit  form  of  eq.  (5.30b)  above,  as  the  general 
relationship  from  the  standpoint  that  each  new  setting  of  the  monochromator,  like  the 
insertion  of  a  new  spectral  filter  in  a  filter  radiometer,  can  be  considered  as  defining  a 
new  instrument  with  a  new  spectral  responsivity.   It  is  the  spectral  responsivity  of  such 
an  instrument,  with  a  fixed  wavelength  setting  or  fixed  filter,  that  we  discuss  in  general 
terms  in  this  summary. 

Unfortunately,  there  is  no  unique  general  solution  to  the  measurement  equation.  Even 
with  R  (x,y,6  ,<|>,A)  completely  specified  with  respect  to  all  parameters,  there  are  unlim- 
ited possibilities  for  L.  (x,y,6  ,<j>,  A)   that  will  satisfy  eq.  (5.30b).   Furthermore,  we 

A 

cannot  replace  L.   in  eq.  (5.30b)  by  the  radiometric  quantity  we  want  to  measure  when  it 

A 

is  an  integral  of  just  L   without  R   unless  eq.  (5.30b)  can  be  modified  to  contain  that 
integral,  which  requires  that  R   be  a  constant  with  respect  to  one  or  more  radiation 

parameters.   In  principle,    if  the  parameter  intervals  are  made  small  enough  so  that  L^ 

A 

varies  negligibly  within  the  intervals,   L   can  be  evaluated.   If  this  is  done  for  enough 

A 

different  parameter  values  to  specify  the  function  L   over  some  range,  the  integrals  over 

A 

that  range  can  be  evaluated.   However,  this  is  usually  extremely  difficult  and  tedious  and 
is  generally  impractical. 

The  practical  approach  to  solving  the  measurement  equation  is  to  select  and  adjust  the 
instrumentation  and  the  measurement  configuration  to  make  possible   (a)  substitution  of  the 
desired  radiometric  quantity  and   (b)  simplification  to  obtain  a  unique  solution.   A  com- 
plete and  exhaustive  procedure  or  algorithm  covering  every  conceivable  case  is  not  possible. 
In  fact,  development  of  this  subject  will  be  our  main  concern  throughout  most  of  the  rest 
of  this  Manual.   For  now,  we  suggest  a  set  of  steps  for  solving  the  measurement  equation, 
eq.  (5.30b),  as  a  good  way  to  start  dealing  with  almost  every  case.   These  steps  cover  a 
systematic  review  of  the  pertinent  features  of  the  measurement  configuration.   They  estab- 
lish approximate  values  of  the  terms  in  eq.  (5.30b)  for  estimating  the  expected  output 
signal  to  assess  the  feasibility  of  a  proposed  measurement  and  the  adequacy  of  proposed 
instrumentation.   They  review  the  constancy  of  R   required  for  transformation  to  the 
desired  radiometric  quantity.   They  deal  with  the  further  simplification,  which  usually 
requires 

that  the  unknown  radiometric  quantity  be  constant  over  one  or  more  of 
the  remaining  integration  intervals  so  that  it  can  be  removed 
from  the  integrals, 

resorting  to  averages,  or 

utilizing  a  relative  spectral  responsivity  function  and  a  relative 
spectral  distribution  (or  distribution  with  respect  to  other 
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radiation  parameters)  of  the  unknown  radiometric  quantity. 

When  adequate  simplification  has  been  achieved,  the  simplified  equations  for  the  unknown, 
and  for  a  calibration  measurement  of  a  standard,  can  be  combined  so  as  to  eliminate  the 
responsivity  and  associated  geometrical  factors,  leaving  the  solution  for  the  unknown 
radiometric  quantity  in  terms  of  just  the  two  measured  output  signals  (unknown  and  standard) 
and  the  known  value  of  the  standard.   Alternatively,  the  equation  for  the  standard  can  be 
solved  to  evaluate  the  responsivity;  then  this  calibrated  responsivity  can  be  used  to 
evaluate  unknown  quantities  from  the  signals  produced  when  they  are  measured  in  the  same 
way.   Finally,  if  a  unique  solution  is  not  possible  because  not  even  the  relative   incident 
distribution  is  known,  results  of  limited  usefulness  and  significance  can  be  obtained  in 
the  form  of  normalized  quantities.   Normalization  is  briefly  discussed  following  problem  3. 

Readers  are  reminded  that  the  measurement  equation  developed  in  this  chapter  is  still 
not  complete.   It  does  not  include  the  radiation  parameters  of  time  and  polarization, 
environmental  and  instrumental  parameters,  corrections  for  diffraction  effects,  or  cases  of 
non-linear  responsivity. 

Regular  use  of  the  measurement-equation  approach  to  radiometry  is  useful  in  planning 
for  and  estimating  the  feasibility  of  a  proposed  measurement,  including  calibration  measure- 
ments and  needed  auxiliary  experiments.   It  provides  a  systematic,  comprehensive  approach 
for  carrying  out  a  measurement  and  for  estimating  its  uncertainty. 
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Appendix  3.    PROJECTED  SOLID  ANGLES,  THROUGHPUTS,  and  CONFIGURATION  FACTORS 
by  Fred  E.  Nicodemus 

When  the  spatial  distribution  of  (radiant)  flux  is  given,  as  we  have  presented  it,  in 
terms  of  ray-radiance,  the  beam-geometry  configurations  and  measures  are  naturally  given  in 
terms  of  beam- intersecting  projected  areas  and  solid  angles  (or  intersecting  areas  and 
projected  solid  angles),  and  throughputs  (see  Appendix  2  [5]).   However,  most  of  the  heat- 
transfer  and  illumination  engineering  literature  deals  with  source  exitance,  irradiance  on 
a  receiver,  and  what  we'll  call  "configuration  factors."  Even  though  the  configuration 
factors  are  usually  defined  in  terms  of  radiant  energy  exchange  (with  simplifying  assump- 
tions that  aren't  always  emphasized  sufficiently),  they  are  actually  purely  geometrical 
quantities  that  can  be  directly  related  to  projected  solid  angles  and  throughputs.   We 
develop  those  relations  in  this  appendix  because  it  is  frequently  useful  to  transform  a 
projected  solid  angle  or  throughput  to  the  equivalent  configuration  factor  so  that  it  can 
be  evaluated  by  means  of  the  fairly  extensive  published  tables  of  configuration  factors, 
thereby  avoiding  otherwise  long  and  difficult  computations. 

The  quantities  involved  in  the  engineering  treatment  are  called  by  a  wide  variety  of 
names:   angle  factor,  angle  ratio,  configuration  factor,  geometrical  factor,  geometrical 
configuration  factor,  interchange  factor,  shape  factor,  or  view  factor,  sometimes  with 
"coefficient"  in  place  of  "factor".   We  prefer  the  term  "geometrical  configuration  factor", 
or  just  "configuration  factor",  used  by  Siegel  and  Howell  [22],  but  our  notation  and 
defining  equations  follow  more  nearly  those  for  the  corresponding  "angle  factor"  of  Sparrow 
and  Cess  [23].   As  we  use  the  term  "configuration  factor",  then,  it  includes  or  is  closely 
related  to  the  CIE-IEC  [6]  terms  and  concepts  of  "configuration  factor"   c,   "form  factor" 
f,   and  "(mutual)  exchange  coefficient"   g.   Although  we  do  not  use  or  define  here  these 
CIE-IEC  quantities,  their  relationships  to  each  quantity  that  we  do  define  is  given  for  the 
benefit  of  those  who  may  want  to  use  them  or  material  involving  them. 

It  is  important  to  understand  the  configuration  factors  not  only  because  they  are  so 
widely  used  in  the  engineering  literature  but  also,  as  already  mentioned,  because  of  the 
existence  of  extensive  tables  of  formulas  [22,23]  that  can  be  used  in  evaluating  them  for  a 
very  wide  variety  of  beam  configurations.   By  knowing  the  exact  relationships  between  those 
factors  and  the  corresponding  projected  solid  angles  and  throughputs,  the  latter  can  also 
be  evaluated  from  the  tables  with  very  substantial  savings  in  time  and  effort.   Further- 
more, the  formulas  in  the  tables  can  be  extended  to  still  other  configurations,  in  many 
instances,  by  the  use  of  the  so-called  "configuration-factor  algebra"  [22]  or  "angle-factor 
algebra"  [23]. 

We  denote  the  configuration  factor  (angle  factor)  for  the  beam  consisting  of  all  rays 

between  points  on  a  reference  surface  of  area  A   and  a  distant  surface  of  area  A.   as 

F   .   or,  more  explicitly,   F     .   Either,  or  both,  of  these  surface  areas  A  may  be 
i— J  A^-Aj 

just  an  element  of  area  dA.   The  reference  surface  is  ordinarily  the  locus  of  the  vertices 
of  the  solid  angles  or  projected  solid  angles  that  combine  with  its  area  to  make  up  the 
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throughput  between  it  and  the  distant  surface  area.   We  use  these  designations,  "reference" 
and  "distant"  in  place  of  the  usual  ones  of  "source"  and  "receiver"  surfaces  to  emphasize 
the  fact  that  we  are  concerned  only  with  purely  geometrical  relationships  between  the  two 
surfaces.   Also,  the  reference  surface  is  not  necessarily,  or  always,  the  source;  sometimes 
it  is  the  receiver  and  the  distant  surface  is  the  source.   Furthermore,  as  indicated  earlier, 
the  configuration  factor  F.  .   is  usually  defined  in  terms  of  the  interchange  of  radiant 
flux.   That  definition  is  as  follows:   when  the  reference  surface  is  perfectly  uniform  and 
diffuse  (lambertian) ,  so  that  all  rays  through  every  point  of  the  reference  surface  have 
the  same  radiance,   F.  .   is  the  fraction  of  the  total  flux  into  (or  from)  the  full  hemis- 
phere  above  the  reference  surface  area  that  intercepts  the  distant  surface  area.   (If  the 
reference  surface  is  an  extended  curved  (non-planar)  surface,  this  refers  to  the  hemisphere 
above  each  surface  element,  taken  element  by  element,  that  makes  up  the  reference  surface 
area,  as  will  be  clear  in  the  defining  relations  given  below.) 

The  purely  mathematical  (geometrical)  defining  relations  for  the  four  possible 

configurations  shown  in  figures  A3-1,  -2,  -3,  and  -4  are  listed,  respectively,  in  eqs. 

(A3-1) ,  (A3-2) ,  (A3-3) ,  and  (A3-4) .   They  are  given,  each  time,  on  the  first  line  as  they 

usually  appear  in  the  engineering  literature.   On  the  second  and  succeeding  lines, 

in  each  case,  are  given  the  relationships  to  the  corresponding  projected  solid  angle 

fi  e  Jcos6'd(j)  =  J  /cosS'sine'de  •<!$     and  to  the  corresponding  throughput   0  =  JJdA'dfi,   as 

well  as  to  the  CIE-IEC  [6]  "configuration  factor"   c,   "form  factor"   f,   and   "(mutual) 

exchange  coefficient"   g.    D   is  always  the  slant  distance  between  the  elements   dA.   and 

dA.. 
3 

(1)   The  configuration  factor  for  the  beam  between  two  infinitesimal  surface  elements 
dA.  (reference)   and  dA.  (distant)   (see  figure  A3-1,  p.  96)  is 

? 
dF,,   JA    =  cosG.  .  'cos0.  •dA./CiT'D  ) 
dA  -dA.         ij     ji   j     s 


i   j 


=   (l/ir)'COs6   •dw..   =   (l/tr)«dn..   =   dc .  . 
ij   ij  ij       ij 


=   [l/(Tr.dA±)]'d0  =   (1/dA.Wg  =   df±.  (A3-1) 

(2)   The  configuration  factor  for  the  beam  between  an  infinitesimal  surface  element 
dA.  (reference)   and  an  extended  surface  area  A.  (distant)   (see  figure  A3-2 ,  p.  97)  is 


dF,A   .    =  /   [cos6.  . -cose.  ./(Tr«D2)]«dA. 
i  j     A     1J     J1     s     J 


=   (1/tO'J   cosG   du  =   (l/ir)'J   d^.      =   (1/ir)  -Q± 


ij  -«         '  'w_  ■  =   °ij 


[l/(iT-dA.)]-de  =   (l/dA±)-dg  =   f±.  (A3-2) 
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(3)   The  configuration  factor  for  the  beam  between  an  extended  surface  area  A. 
(reference)   and  an  infinitesimal  surface  element   dA.  (distant)   (see  figure  A3-3,  p.  98) 
is 

dFA   JA    =   (1/A.)-/   [cose..  -cosG.  .-dA./(ir-D2)]  -dA. 
A±-dA  1   A      1J     J1   J     s     x 

J  i 

=   [l/(ir«A.)]'f   (cosO.  .  -dcj..)-dA.   =   [1/ (tt-A.)  ]  •/  dfl. .  «dA. 
1    JA      ij    ij     l  x    JA        IJ    1 

i  i 


=   (1/tt)  *dfi.  .   =   dc.   [both  averaged  over  the  area  A.] 
ij       ij  b  iJ 

=   [l/(ir«A  )]«d0  =   (l/A.)-dg  =  df . .  (A3-3) 

(4)   The  configuration  factor  for  the  beam  between  two  extended  surface  areas  A. 
(reference)   and  A.   (distant)   (see  figure  A3-4,  p.  98)  is 

FA   .    =   (l/A.)-f   J   [cosO..  «cose.  ./(ir'D2)]«dA.«dA. 
A. -A.         iJAJA      ii     li     s     j    l 
l  J  A.  A.       J      J  J 

i  J 

=   [l/(ir«A  )]  •/   /   cos9.. -doj.  . -dA.   =   [1/(tt-A.)  ]  •/  Q    . -dA. 
1    A.  (a..     1J    1J    X  1    A.   J    1 

=   (1/tt)«Q..   =   c..   [both  averaged  over  the  area  A.] 

=   [1/(1^).]  -0  =   (l/A±)-g  =   f±.  (A3-4) 

It  is  also  easily  shown  that  the  following  reciprocity  relations  hold  with  respect  to 
all  of  the  configuration  factors: 

dA.-dF,.   JA   =   dA.'dF,.   ,.  (A3-la) 

l   dA.-dA.        i    dA.-dA. 
i   J  J    i 

dAi'dFdA-A.   "   VdFA.-dA.  (A3"2a) 

i  J  1        i 

VdFA.-dA.   "   ^j'^dA.-A.  (A3"3a) 

i   J  J   i 

A.«FA   A    =  A.-FA   ,  (A3-4a) 

l   A. -A.       i   A. -A. 
i  J  J   i 

Note  that  eqs.  (A3-2a)  and  (A3-3a)  are  redundant.   Both  are  included  here  only  to  maintain 
the  one-to-one  relationship  between  this  set  of  equations  and  the  preceding  set  of  defining 
equations  [eqs.  (A3-1) ,  (A3-2) ,  (A3-3) ,  and  (A3-4) ] . 
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In  the  engineering  literature,  as  previously  mentioned,  there  are  extensive  tables 
[22,23]  of  configuration  factors,  or  equivalent  quantities  by  other  names  (most  of  which 
are  listed  in  the  second  paragraph  of  this  appendix) .   These  can  also  be  readily  combined 
to  obtain  the  configuration  factors  for  still  other  combinations  of  surfaces  by  application 
of  the  so-called  "configuration-factor  algebra"  [22]  or  "angle-factor  algebra"  [23]  that 
has  been  developed.   Accordingly,  by  expressing  projected  solid  angles  and  throughputs  in 
terms  of  configuration  factors,  in  other  words,  by  reversing  the  relations  of  eq.  (A3-1) , 
(A3-2) ,  (A3-3) ,  and  (A3-4) ,  the  tables  of  configuration  factors  can  be  used  to  evaluate 
many  projected  solid  angles  and  throughputs  without  carrying  out  the  integration  for  a 
projected  solid  angle  (eq.  (2.31)  [5] 


n 


/  cos0*dw  =   / /cosQ'sine'dO'dcf)  [sr] 


(A3-5) 


or  a  throughput  (eq.  (2.27)  [5] 


0   =   J/dfi'dA  =   ////cose -sine -d6 «d(fi'dx'dy  [m  -sr] 


(A3-6) 


which,  while  they  are  conceptually  straightforward,  often  entail  complex  and  difficult 
computations.   The  required  inverse  relations,  needed  for  using  the  configuration-factor 
tables  to  evaluate  a  projected  solid  angle  or  throughput,  are  given  in  eqs.  (A3-7)  to 
(A3-12) ,  inclusive: 


, —  dA 


dA:  —J 


Figure  A3-1.    Configuration  for  interchange  between 
two  infinitesimal  surface  elements. 

(a)   The  element  of  projected  solid  angle  subtended  at  a  surface  element   dA.   by  a 
surface  element  dA,   (see  figure  A3-1)  is 
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dfl..   =   TT-dF,A       =   Tr-(dA./dA.)-dF,.      [sr] 
11         dA.-dA.  i    1    dA.-dA. 

i   J  J    i 


(A3-7) 


The  corresponding  element  of  throughput  between  dA.   and  dA.   is 

? 
d0  =   dA.-dfi..   =   ir'dA.'dF,.   ,.    =   TT«dA.«dF,A   , .   [m  -sr] 
l    it  l   dA.-dA.  i    dA.-dA. 

i   3  J    i 


(A3-8) 


i— dA 


Figure  A3-2.    Configuration  for  interchange  between  an  infini- 
tesimal surface  element  and  an  extended  surface. 

(b)   The  projected  solid  angle  subtended  at  a  surface  element   dA.   by  an  extended 
surface  area  A.   (see  figure  A3-2)  is 


a--      ■  v^ah      a   ■  TT'(A./dA.)'dF.   ,.   [sr] 
xj        dA.-A.  ii    A.-dA. 

i  J  3i 


(A3-9) 


The  corresponding  element  of  throughput  between  dA.   and  A.   is 

dQ  ~     dA.'fi..   =   ir-dA .'FJA   A   =   tt-A.  «dF.   ,.   [m  -sr] . 
l   ii         i  dA.-A.        i   A.-dA. 
l  j  j    l 


(A3-10) 
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r-dA 


Figure  A3-3.    Configuration  for  interchange 

between  an  extended  surface  and  an 
infinitesimal  surface  element. 


(c)   The  element  of  throughput  between  an  extended  surface  area  A.   and  a  surface 
element   dA.   (see  figure  A3-3)  is 


d0  =   *'A  -dF        =   ir-dAj-F      [■  -sr] 


(A3-11) 


r-dA 


Figure  A3-4.    Configuration  for  interchange 
between  two  extended  surfaces. 


98 


(d)   The  throughput  between  two  extended  surface  areas  A.   and  A.   (see  figure  A3-4) 


0  =   tt«A.*F.   .    =   Tr«A.'F.   .   [m  «sr] 
i  A. -A.        i   A. -A. 
i  J  J   i 


(A3-12) 


Figure  A3-5.    Two  parallel,  coaxial  circular  discs. 


2     A        2 
A.  =  ttt  . ;   A.  =  ir«r . ; 

i      i     J      J 


D   is  the  perpendicular  distance  between 
centers  along  the  common  axis. 


We  will  illustrate  the  use  of  configuration-factor  tables  by  obtaining  the  formula  for 
the  throughput  of  a  very  common  beam  configuration,  the  beam  between  two  parallel  coaxial 
circles.   Most  optical  systems  have  circular  elements.   In  particular,  they  employ  circular 
stops,  so  the  images  of  those  stops,  the  pupils  and  windows,1  are  also  circular.   Accord- 
ingly, it  is  quite  common  for  a  beam  to  be  defined  by  two  parallel,  coaxial,  circular  areas 
A.   and  A.   separated  by  a  distance  D  along  the  axis  between  their  centers,  as  shown  in 

figure  A3-5.   For  example,  in  figure  4.10,   A.  =  AA   and  A.  =  AA   or,  again,  in  figure 

x     s        j     r 

4.13,  we  might  designate  the  beam-defining  apertures  as  A  =  A^  and  A.  =  A   or  else  as 
A  =  A   and  A.  =  A„.   In  any  case,  returning  to  figure  A3-5,  we  designate  the  radius  of 


See  any  standard  textbook  on  geometrical  optics  [13,14]. 
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2  2 

A.   as   r.   and  that  of  A.   as   r.   so  that  A.  =  ttt.   and  A.  =  ttt..   We  will  use  the 

11  JJ  iiJJ 

first  part   of   eq.    (A3-12) 

G     =      vA-I     _A      [m2-sr],  (A3-12a) 

i     J 

so  what  we  need  from  the  tables  is  the  formula  or  equation  for  evaluating  the  configuration 

factor  F       for  the  configuration  of  figure  A3-5. 

i~  3 
In  the  "Angle  Factor  Catalogue"  in  Appendix  A  of  [23],  we  find  the  desired  geometrical 

configuration  illustrated  as  "configuration  3"  with  the  following  entry  (adjusting  notation 

to  correspond  with  figure  A3-5) : 

X  =  r./D,   Y  =  D/r.,   Z  =  1  +  (1  +  X2)Y2, 

FA  _A   =   (1/2)  [Z  -  (Z2  -  bxh2)1"2].  (A3-13) 

i  3 

Alternatively,  in  the  "Catalog  of  Selected  Configuration  Factors"  in  Appendix  C  of  [22],  we 
also  find  the  desired  geometrical  configuration  illustrated  in  item  18,  with  the  following 
entry  (again  we  make  straightforward  substitutions  of  our  notation) : 

R.  =  r./D,    R.  =  r./D,   X  =  1  +  (1  +  R2)/R2, 
i    i       J    J  J    i 

FA  _A   =   (1/2) {X  -  [X2  -  4(R  /R±)2]h}.  (A3-14) 

i  j 

Starting  with  either  formula,  eq.  (A3-13)  or  (A3-14) ,  for  F  _   ,   we  can  obtain,  after 

i  J 
some  algebraic  manipulation,  the  following  expression  for  the  desired  throughput: 

0  =   (w2/2)-{(r2  +  r2  +  D2)  -  [ (r2  +  r2  +  D2)2  -  4r?-r*]'S}  [m2-sr].  (A3-15) 
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Appendix  4.    Some  Nomenclature  Considerations  —  Signal   S,   Responsivity  R   [or  ft], 
Reflectance  Factor  R.   —  by  Fred  E.  Nicodemus 

Dr.  G.  Bauer,  Chairman  of  the  Vocabulary  Subcommittee  of  the  CIE  Technical  Committee 
TC-1.2,   Photometry  and  Radiometry,   has  pointed  out  some  instances  where  our  nomenclature 
departs  from,  and  may  conflict  with,  the  CIE-IEC  International  Lighting  Vocabulary  [6]  and 
some  revisions  to  that  document  that  are  being  proposed  by  his  Committee.   First,  although 
our  use  of  the  term  "signal"  and  the  symbol  S   (or  s   for  the  relative  quantity)  to 
designate,  in  general,  the  output  of  a  radiometer  or  radiometric  instrument  is  based  on  the 
IEC  [20],  the  symbol   S,   at  least,  is  also  used  by  the  CIE-IEC  [6]  to  denote  the  "relative 
spectral  energy  [power]  distribution  (of  radiation)".   To  date,  we  have  not  had  occasion  to 
make  much  use  of  that  general  designation,  since  we  are  usually  concerned  with  the  spectral 

distribution  of  a  specific  radiometric  quantity,  e.g.,  spectral  radiance  L.(A)   or  spec- 

A 

tral  radiant  flux   $  (A).   Furthermore,  when  we  do  want  to  refer  to  such  quantities,  in 

A 

general,  without  designating  any  one  of  them  specifically,  we  prefer  to  let  the  symbol  X 
designate  any  such  radiometric  quantity,   $,  I,  E,  L,  etc.   as  in  [21].   Then  the  spectral 
distribution  is  given  as  X.(A),   where  the  subscript   A   is  a  clear  reminder  that  this 
quantity  is  a  derivative  with  respect  to  the  spectral  variable,  .something  that  is  not 
brought  out  by  the  CIE-IEC  symbol   S(A)   for  a  spectral  distribution. 

In  choosing  to  use  the  IEC  term  "signal"  and  symbol   S   for  the  output  of  a  radiometer, 
we  are  strongly  motivated  by  the  fact  that  this  is  the  "language"  of  communications  engi- 
neering.  We  feel  that  it  is  well  to  recognize  that  such  an  output  is  indeed  a  "signal" 
that  is  used  to  convey  information  about  the  input  of  incident  radiation.   In  fact,  that  is 
the  most  significant  thing  about  it,  regardless  of  whether  it  is  a  voltage,  a  current,  a 
photographic  film  density,  or  a  displacement  of  a  meter  needle  or  of  a  line  on  a  strip- 
chart  record.   Furthermore,  the  use  of  this  nomenclature  and  point  of  view  may  facilitate 
access  to  powerful  analysis  tools  in  the  literature  of  communications  engineering  that  may 
usefully  be  brought  to  bear  on  some  radiometric  measurement  problems. 

Dr.  Bauer  also  notes  that  they  wish  to  recommend  the  use  of  the  symbol   s   for  respon- 
sivity (sensitivity)  and  that  our  use  of  R  and   r   for  responsivity  may  conflict  with  the 
use  of  R  for  "reflectance  factor"  by  the  CIE-IEC  [6].   We  strongly  support  the  position 
of  R.  Clark  Jones  [24],  that  the  use  of  the  term  sensitivity  is  to  be  deprecated  because  of 
its  ambiguity  and  that  the  more  explicit  and  unambiguous  terms  "responsivity"  [output  per 
unit  input]  and  "detectivity"  [reciprocal  of  noise-equivalent  input]  should  always  be  used 
in  its  place.   Accordingly,  we  prefer  to  avoid  even  the  symbol   s   for  responsivity,  because 
of  its  association  with  the  undesirable  term  "sensitivity."  Furthermore,  we  do  not  antici- 
pate serious  difficulty  with  the  possible  confusion  between  R  for  responsivity  (and   r 
for  relative  responsivity)  and  for  reflectance  factor.   The  context  will  usually  make  it 
quite  clear  which  use  is  intended.   Also,  while  both  may  have  angular  or  directional 
dependence,  the  directional  dependence  of  the  reflectance  factor  always  involves  two  sets 
of  directional  designations,  those  for  both  the  incident  and  reflected  ray  or  beam  geometry, 
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e.g.,   R(0.,<|).;6  ,  <j>  )   or,  most  generally,   R(u.;co  )  [26].   The  directional  dependence  of 
responsivity,  on  the  other  hand,  involves  only  the  direction  of  the  incident  rays  R(6,<j>). 
Neither,  of  course,  is  a  derivative;  both  are  weighting  functions  or  weighting  factors  with 
respect  to  the  radiometric  quantities  and  related  signals.   In  those  instances  where  it  may 
be  necessary  to  use  symbols  for  both  responsivity  and  reflectance  factor  in  the  same  dis- 
cussion and  the  same  equations,  we  would  probably  shift  to  a  script   R,  for  responsivity, 
in  order  to  keep  them  separate,  but  we  do  not  believe  that  this  will  occur  often  enough  to 
make  it  worthwhile  to  use  the,  otherwise  less  convenient,  script   R,  for  responsivity  all 
of  the  time. 

We  are  grateful  to  Dr.  Bauer  for  his  comments  and  we  understand  that  there  may  be 
quite  different  and,  for  them,  valid  reasons  for  the  positions  taken  by  him  and  his  com- 
mittee.  However,  for  the  reasons  given  above,  we  believe  that  we  need  to  take  the  fore- 
going positions  for  this  tutorial  presentation.   We  urge  that  these  reasons  also  be  given 
consideration  by  those  currently  engaged  in  revising  the  CIE-IEC  International  Lighting 
Vocabulary  [6].   It  should  be  emphasized  that  these  are  the  views  of  the  author,  only,  and 
not  necessarily  those  of  the  National  Bureau  of  Standards  or  of  any  other  group  with  which 
he  may  be  associated. 
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NBS  TECHNICAL  NOTE  910-1  —  ERRATA 

Reverse  of  Title  Page,  last  line  of  footnote: 
price  should  be  given  as  $2.10. 

p.  8,  next  to  last  line: 

"directed"  should  read  "detected". 

p.  11 ,  Figure  2.1 : 

ribbon  filament  should  be  extended  downward  so  that  point  1  is  at  the  middle  of 
the  filament. 

p.  21,  last  line  of  footnote: 

first  word  should  be  "wish"  (not  "with") . 

p.  25,  2nd  line  of  eq.  (2.15): 

2 
insert   6   so  that  it  reads:   "=  L  -cosG  «dA  -cose  *dA  /D   [W],". 

p.  32,  eq.  (2.24): 

insert  multiplication  dot  between  dw  and  dA. 

p.  32,  2nd  line  after  eq.  (2.24): 
insert  "e"  in  "includes". 

p.  52,  1st  line  of  eq.  (3.10): 

change  "sinO"   to   "cos6", 

p.  52,  1st  line  after  eq.  (3.10): 

change  "eqs.  (2.23)  and  (2.28)"  to  "eqs.  (2.24)  and  (2.29)." 

p.  53,  1st  and  3rd  lines: 

change  "eq.  (2.24)"  to  "eq.  (2.25)". 

p.  53,  line  9: 

change  "eqs.  (2.25)  through  (2.33)"  to  "eqs.  (2.24)  through  (2.29)". 

p.  69,  next  to  last  line  of  paragraph  following  eq.  (A2-8): 
change  "f  "   to  read  "F  ". 
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